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Energy Accounts of Earthquakes in and near Japan. 


By 


Chuji Tsusor 


Geophysical Institute, Faculty of Science, Tokyo University, Tokyo. 


Abstract 


By using the newly deduced formula: 


M=1.73 log 4+log A—0.83, 


where M is the magnitude of an earthquake, and A is the maxium ground amplitude 


measured in micron at the epicentral distance 4 measured in km., 


the magnitudes of 1025 


earthquakes which took place from 1931 through 1955 in and near Japan have been determined. 
The mean annual number N of earthquakes according to the magnitude has been 


found to be expressed by 


log N= —1.08+0.72 (8—M), 


‘for 4M=0.1. 


The mean rate of energy release by earthquakes in and near Japan has been estimated 
to be 1.510% ergs/cm? year or 0.1 x 10-® calory/cm? sec. 
As a whole, Japanese area is contributing well over 10% to the total seismicity of the 


world. 


The total amount of energy which is sent 
out in the form of earthquake waves from a 
certain region during a certain span of time 
can be estimated if the energy of each in- 
dividual earthquake concerned is known from 
seismogram studies. But’ this logically 
straightforward procedure is well known in- 
evitably to introduce extremely laborious 
calculations rather than desirable simplifica- 
tions. For statistical investigations, it is far 
better and far more practical first to deter- 
mine the magnitude of each individual earth- 
quake and then to convert the magnitude into 
energy on the basis of some appropriate 
formula. This is indeed the procedure which 
has been followed by many authors in the 
studies of this kind. 

The magnitude of an earthquake was ori- 
ginally defined by C. F. Ricuter (1935) 
relation to the maximum trace amplitude on 
the seismogram obtained by a standard tor- 
sion seismometer having prescribed instru- 
mental constants. For this very useful con- 
cept of magnitude to find wider applications, 


it is desirable to have a fofmula deduced in 


which the magnitude of an earthquake is ex- 


pressed in terms of the maximum ground 
amplitude observed at a point, instead of the 
maximum trace amplitude recorded by. the 
“standard” seismometer, and the epicentral . 
distance of that point. This is especially 
disirable for studying Japanese earthquakes 
because there are great many of them and 
because they are observed ae PASEO AO 
other than the “ standard” one. 

In one of his previous papers (Tsusot, 1952), 
the present author deduced formulae for this 
purpose. To begin with, the form 


M=a4+8 log A+r 


was adopted, with several grounds for justi- 
fication, where 4 is the epicentral distance, A 
the maximum ground amplitude, and a, 8 
and 7 are numerical constants. In order to 
determine the numerical values of the three 
constants in the formula, the following pro- 
cedure was adopted. From the magnitude 
list of earthquakes given in GUTENBERG 
and Ricuter’s book “Seismicity of the 
Earth ”, 78 Japanese shocks were chosen and 
the values ‘of the three constants’ in’ the 


-formula were determined by the method-of 
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least squares in such a way that the formula than 500km. The two formulae obtained are 
will, with 4 and A observed at Japanese as follows: 
seismological stations, give as close the M=0.204 +0.62 log A+3.80 (4<5), 


magnitudes as possible to those given in the = 

book. It was found that two different for- M=0.0334-+0.60 log A+5.005 (4>5) 
mulae are needed for the purpose according 
as the epicentral distance is greater or less micron. 


where 4 is measured in 100 km. and A in 
The first formula is appicable for 


Table I. Yearly Number of Earthquakes 
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earthquakes of which the epicentral distances 
are less than 500 km. and of which the maxi- 
mum ground amplitude is presumably due to S 
waves. The second formula, on the other hand, 
is applicable for earthquakes of which the dis- 
tances are greater than 500 km. and of which 
the maximum ground amplitude is presumably 


in and near Japan according to Magnitude 


due to surface waves. By means of these 
formulae, the magnitudes of 735 shallow 
earthquakes were determined which had taken 
place during 20 years from 1931 through 1950 
in and near Japan, the results having been 
published by the Central Meteorological Ob- 
servatory (1950). On the basis of the deter- 
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minations, some aspects of the magnitude- 
frequency relation of earthquakes in and near 
Japan were studied by the author (Tsusot, 
1952). : 

“In 1953, R. Yosuryama (1953) rightly point- 
ed out that (Mce—Mr) shows systematic in- 
crease with increasing Mz, where Mz is the 
magnitude of an earthquake determined by 
GuTenserRG and Ricuter and Mr, that of the 


same earthquake determined by the present — 


author. The reason for this systematic 
tendency has been searched and finally it was 
found (Tsuso1, 1954) that the form of the 
formula assumed at first was not very satis- 
factory and that the form 


M=a log 4+log A+c 


is better and more reasonable. There are 
two numerical constants to be determined in 
the new formula. By means of the procedure 
similar to what was followed in the previous 
case, the values of the two constants have 
been determined with the result 


M=1.73 log 4+log A—0.83. 


In this, 4 is measured in km. and A in mic- 
ron. There remain no systematic discre- 
pancies between the magnitudes determined 
by Gurenserc and Ricuter and those deter- 
mined by means of the above formula. Also 
no need was felt any more to use two dif- 
ferent formulae according to epicentral dis- 
tances. 

On the basis of this new formula, the 
magnitudes of 1025 earthquakes have been 
determined. Earthquakes taken were those 
which occurred in and near Japan from 1931 
through 1955 and which were felt beyond 
200 km. from respective epicentres. In order 
to determine the magnitudes of individual 
earthquakes, observational materials obtained 
at the six seismological stations, Sapporo, 
Sendai, Tokyo, Nagoya, Kobe and Fukuoka, 
have been used. Making use of observational 
materials obtained at the six stations, six in- 
dependent estimates of the magnitude of one 
andthe same earthquake .are obtained which 
do not generally agree exactly with one ano- 


ther. The mean value of the magnitudes has' 
been taken to be the final magnitude of that 
earthquake. “ *t 
Table I gives the yearly numbers of earth- 
quakes according to their magnitudes as 
determined by the above procedure, while 
Fig. 1 shows the mean annual numbers of 
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Fig. 1. Mean Annual Number of Earthquakes 
in and near Japan according to Magnitude. 
(1931~1955). 


earthquakes according to their magnitudes. 
The decrease of the number toward smaller 
magnitudes than M=6 is of course merely an 
apparent tendency caused by our taking those 
earthquakes only which were felt beyond 
200 km. from the origins and is therefore re- 
presenting no real situation. The portion of 
the curve in Fig. 1 that shows the decrease 
of the number toward larger magnitudes than 
M=6 can be trusted upon. Least squares cal- ( 
culations have shown that this portion of the 
curve can be well expressed by 


log N=—1.08+0.72 (8—M), 


the class interval for M being taken as 0.1. 
It is interesting to compare this formula with 
the corresponding one deduced by B. GuTEN- 
BERG and C. F. RicuTer (1954) for the whole 
world: 


log N=—0.48+0.90 (8—M). 


The mean annual numbers of earthquakes 
which have larger magnitudes than M in the 
whole world and in Japanese area and which 
are expected from the above formulae have 
been computed and are compared in Table II. 
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Fig. 2. Cumulative Sum of Released Energy in and near’ Japan. 


From the values given in the table, it may 


be said that earthquake frequency in and 
near Japan is well over 10% of that in the 


whole world taken altogether. 
Table II. Mean Annual Number of Earthquakes 
Having Magnitudes larger than M. 
M Whole world ‘| Japan 
6.0 a+110: B+14.0 
6.5 at 37.8 B+ 5.9 
TeLD a+ 12.3 B+ 2.3 
7.5 a+ 3/3 BY'O.7 
A ae a | B 
Now that the magnitude of each of the 
earthquakes in and near Japan has been 


determined, the computation of the cumula- 


tive sum of released energy is a matter of 


simple arithmetic. In doing this, however, 
the final result will decidedly be affected by 
the choice of formula which gives EF asa 
function of M. Here the 1956 formula of B. 
GUTENBERG and others (1956): 


log H=11.8+1.5 M 


will be used. Fig. 2 shows the cumulative 
sum of released energy in and near Japan. 
The increase in energy with time is notably 
uniform. The final value reaches about 
1.2 10° ergs in 25 years, so that the mean 
annual release of energy is about 0.5 x10 
ergs/year. If the total area concerned is taken 
to be 3.5x10!°cm?, the mean rate of release 
per unit area will be 1.5 x 10% ergs/cm? year or 
lx OCP cal cm? sec. 

The annual release of earthquake energy in 
the whole world has been estimated by B. 
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GuTenBerG (1956) to be 10% ergs/year. As 
compared with this value, the contribution 
from Japanese area is only 1/20 of it which 
appears to be too low. This is likely to have 
been caused by the difference of the formulae 
connecting M and £ as used by GUTENBERG 
and by the present author so that the above two 
values might not have been to be compared. 
To deduce results which are truly comparable 
with each other, computations must of course 
be carried out on the same basis. It has been 
stated that the mean annual numbers of 
earthquakes within the magnitude range 
M~ (M=+0.05) in the whole world and in 
Japanese area are expressed by 


log N= —0.48+0.90 (8—M) 


=6.72—0.90 M Whole world, 
log N=—1.08+0.72 (8—M) 

=4,68—0.72 M Japan. 
Table III. Mean Annual Energy Released by 


Earthquakes Having Magnitudes larger than M. 
(in 1022 ergs) 


M Whole world Japan 
6.0 362 +a 93+ 
1 361+-a@ 93+, 
2 359 +a 92+ 8 
3 358 + @ 92+ 
4 356 + a 92+ 
5 353+ a@ 92+. 
6 351 +a 91+ 8 
if 348 +a 91+, 
8 344+ a 90+ 8 
9 340 +a 90+ 8 
30 336 + a@ 89+ p 
1 330+ a@ 88+-B 
2 324+ a 87+ 8 
3 317 +a 86+-B 
4 309 + a 84+ 
5 300+ a@ 83+, 
6 290 + @ 80+ 8 
@ 278 +a 78+ 8 
8 264-+a 75+ 8 
9 248 +a 71+ 
8.0 230 +a 67+-B 
1 209 +a 62+, 
y, 185+ a 55+ 8 
3 158+ a 48+ 8 
4 126+a 39+ 8 
5 90+ a 28+ 
6 48 +a 15+ 


If these formulae are combined with 
log E=11.84+1.5 M, 
the expressions 


Whole wnat 
Japan 


log NE=18.52+0.60 M 
log NE=16.48+0.78 M 


are obtained. By means of these formulae, 
the mean annual energy released by earth- 
quakes which have magnitude larger than M 
has been calculated. The results are given 
in Table III. 


From the table, it may be said that the con- 
tribution from Japanese area is roughly 1/4 
of the total amount for the whole world. 
This conclusion agrees with that drawn 
regarding the comparative earthquake fre- 
quencies. 


It will hardly be necessary to add that 
while the conclusions arrived at as to the 
seismicity in and near Japan relative to that 
in the whole world are reliable, the absolute 
values of energy and related quantities deduc- 
ed in this article are uncertain, possibly by a 
factor 10. 


References 


Central Meteorological Observatory: 

1950 “The Magnitude Catalogue of Major Earth- 
quakes which Occurred in Vicinity of 
Japan.” Seismological Bulletin of Central 
Meteorological Observatory for the Year 
1950. 

GUTENBERG, B.: 

1956 “The Energy of Earthquakes.” Quarterly 

Journ. Geol. Soc. London, 112, 1-14. 
GUTENBERG, B. and RICHTER, C. F.: 

1954 “Seismicity of the Earth.” 

versity Press. 
GUTENBERG, B. and others: 

1956 “Progress Report, Seismological Labora- 
tory, California Institute of Technology, 
1955.” Trans. Amer. Geophys. 
Simeeoue 

RICHTER, C. F.: 

1935 “An Instrumental Earthquake Magnitude 

Scale.” Bull. Seism, Soc, Am., 25, 1-32, 


Princeton Uni- 


Union, 


Energy Accounts of Earthquakes in and near Japan if 


TsuBol, C.: 


1952 


“Determination of the RICHTER-GUTEN- 
BERG’s instrumental magnitudes of earth- 
quakes occurring in and near Japan.” 
Geophys. Notes, Geophys. Inst. Tokyo 
Univ., 4, No. 5. 

“ Magnitude-Frequency Relation for Earth- 
quakes in and near Japan.” Journ. Phys. 
Earth, 1, 47-54. 


1954 


“Determination of the GUTENBERG- 
RICHTER’s Magnitudes of Earthquakes oc- 
curring in and near Japan.” Zishin (Journal 
Seism. Soc. Japan) 7, 185-193. 


YOSHIY AMA, R.: 


1953. At the Annual Meeting of the Seismolo- 


gical Society of Japan. 


. 7 
t senegal, “gt tb ni entaged 


watevatao 342 Ta mittienicavatoed nike 
wr adnan d Sustain 2 RSTHONA 
pega nifi, © fe wen piso i 


webi eatad oith ‘a yaiesold Taunt: oft 4 CEPT 
sue to wii Savi, 


by the Inter 
in 2024 


Vo tnewtone Ogee 
“e 


7 ~ 
ws 


oe 
rv 6) Lory 
ee fee Ay = 


ye, Lip fag “¢ 


=m 


tne ike ON mend gd 
https://archive.org/details/journal-of-physics-of-the-earth 1987-05 5 sh tes 


« 


Se 


is 


JOURNAL OF PHYSICS OF THE EARTH, VOL. 5 No. L957 9 


Observations of Tidal Strains of the Earth 


By 


Izuo OzAwa 


Disaster Prevention Research Institute, Kyoto University 


Summary 


Observations of tidal strain components of the earth’s surface by means of extensometers 
were initiated by K. SASSA and the value of constant J, usually called SHIDA’s number, was 
obtained seperately from other constants h and k. 

Observations of tidal strain components in three horizontal, one vertical and two dia- 
gonal directions are made at Osakayama Observatory. Using the observed values, horizontal 
areal strain, vertical strain and cubical dilatation are obtained, and the numerical values 
of h—3l, v and C=v+2h—6l are caluculated therefrom. 


Introduction 


Observations of strain of the earth’s surface 
caused by the tide-generating force by means 
of extensometer were initiated by Kenzo Sassa, 
and the value of a constant 7, called Sxtpa’s 
number was obtained thereby. Eiichi NisHimura 
(1950) obtained value of D=1+k—h=0.66 from 
the tiltmetric observation at Barin and value of 
L=1+k—l=1.20 by analyzing tidal latitude 
variations at six stations. Kenzo Sassa, Izuo 
Ozawa and Soji Yosutkawa (1952) made obser- 
vations of the earth tide by means of extenso- 
meters devised by K. Sassa at Ikuno and 
) Osakayama, and obtained the value of the 
constant 7. I. Ozawa examined the results of 
the observations at Osakayama and Makimine 
and studied the influence of the oceanic-tide. 
Michio Taxapa observed the tidal variation of 
vertical strain at Ide. Hitoshi Takeucai cal- 
culated tidal strains of the earth with internal 
constitutions inferred from recent studies of 
seismology combined with K. E. BUuLLen’s 
studies of density distribution within the earth. 
Kyoto university is conducting this kind of 
observation and research at Osakayama, Ide, 
Kishu, Suhara and several other places. 


Fundamental theory 
The components of radial, meridional and 
longitudinal displacement #,, me and wg due 


to the potential of tide-generating force W, 
are given as follows. 
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where F(r7) and G(r) are functions of 7 only, 7 
is the distance from the center of the earth, 
6, colatitude and ¢, longitude. The disturbing 
potential W, is taken to be proportional to 
S.7?, where S; is a surface spherical harmonic 


of order 2 and simple harmonic function of 


time. 
Making use of (1) and the partial differential 
equation 
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we find that the conponents of linear and longi- 
tudinal strains due to the disturbing potential 
are written as follows 
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where @9+éy¢ is horizontal areal strain and 4 
cubical dilatation. Extension in any direction 
E; will be given by 
Ei = Cree? + 9082+ Cer? + €rea@B + Coph 
+Erg7A , (4) 
where a, 8 andy denote the direction cosines 
of the direction 7 with respect to 7, @ and ¢ 
directions respectively. We can obtain the 
components of the strain by observing £&; 
and solving equation (4). 
Introducing the notations 
h=H(a) 
1) 
h_ 2F(a)+G(aj@ 
gi a 
hey dia) 
se F 
where h is Love’s number, J, Shida’s number, 
a, radius of the earth and g, gravitational ac- 
celeration at the earth’s surface, we have at 
the earth’s surface 


Table 1. 
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where v and C are constants. 
Observation 
Osakayama observatory is situated at 


34°59’.6N and 135°51’.5E. The observing 
room is located 150 meters under the earth’s 
surface and in the middle of subway of 700 
meters in length. No artifical disturbance 
prevents us from our observation. The neigh- 
bouring ground is the clay-slate belonging to 
the Chichibu palaeozoic system. The ampli- 
tude of temperature variation in the observing 
room is about 0.4°C in annual variation and 
smaller than 0.01°C in daily variation, and the 
amplitude of linear strain is the order of 
10-%c.g.s. in daily variation. At this obser- 
vatory, extensometers are placed in six direc- 
tions—three horizontal directions, one vertical 
and two diagonal directions. Two additional 
extensometers are set up in order to make 
comparative observations in particular direc- 
tions. The types and sensitivities of the 
instruments which are related to the present 
paper are shown in Table l. 

The method of analysis is as follows. We 
first read the curve on the photographic re- 
cording paper every hour and by taking over- 
running mean of every twenty-five hour, 
extracted variations in one day, 1/2 day...... 
from which the components My, S:, Ki, O1 
etc. were derived by harmonic analysis of one 


Sensitivity of Instruments. 


nn 


Direction of 


Observatin Type of Instrument 


Azimuth N 38° E Horizontal Sassa-type 
Azimuth N 29° W Horizontal Sassa-type 
Azimuth N 61° E Horizontal Spring-type 
Vertical Spring-type 

Azimuth S 52° W 

Fievationsla® Rod-pendulum-ty pe 
Azimuth N 38° W ‘ 

Blevakion 43° Wire-roller-type 


Length of Sensitivity 
Observation in 10-%c.g.s./mm 
20.0 m. | 0.41—0.63 
4.2m. 3.0 —9.8 

9.6 1.82 
4.4 1.85—2,24 
By 1.83 
6.4 2.46 


a 
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Table 2. Analysed Values for Component Tides. 
—E leo ine) | ww eo a) 
iS Direction of esta | so'as ca Sa | Sa es 
n | = Q, 2) a, a, jo) 
2 Epoch Analysed | 2%) ¢ Uren So wie I = 
an Observation | 22. | Sot | Sot | Seon | Sot | Sot 
A o | OeadS | OaUS |OURUGS | OBUS | OBGCO 
OA | GES. | LEST | GEST | Geass | Less 
| ] g | ad S| ad Ss | Mas |e s | od? 
é | 1947 Oct. 24. |0.41- adits era 
Be) aNeg8" E 1949 Feb. 24. 0.63|9-33 48 
| O66 254 0b 8s 2a 
° 1953 Feb. 19.- ° | 
2 N 29° W cee 8 | 108828 | 
O22 ool OOmmre 
° 1953 Feb. 12.- ° 8 
3 | N6IlE Mar. 12. | 3-1 {0-50 9°) 
. 1952 Dec. 25.— 5 
4 | Vertical 1953 Jan. 23, 2-24 |0.44 215 
| | 
| 4 1953 Feb. 12.- | ° | 
5 | Vertical Mar tis. | 2°22 0.63 201 | | 
: 1953 Mar. 12.- o | 
6 | Vertical Apr. 10. | 2.22 |0.64 193 | 
1954 May 2.- * | s : - 
Th | N 29° W May 331. fe-? LESS O Onin S4 Oral O24) ee ioye || Apt e)— oH! 
8 | N61°E 1954 ad 35° | 1-82 |0.56 14°/0.13 30° 0.04 30°/0.18 359°/0.36 29° 
9 Vertical W954 Oct. 7 | 1.85 0.74 174° 0.65 158°| 0.16 158°|0.73 187°] 0.61 232° 
| 
10 Vertical 1954 Nov. 5 | 1.85 [0.84 192°| 0.39 142°|0.11 142°] 0.62 193° 
. . | | 
| | 
Azimuth S 52° E 1955 Jan. 7.- o| 
1 Elevation 45° Jans 21% 1.34 |0.12 56 | 
Azimuth S 52° E | 1955 Sep. 20,- ° 
12 | Elevation 45° | Octal 9: 1.34 | 0.06 271 | 
| Azimuth N 38° E | 1956 Jul. 13.- | o | 
13 | Elevation 48° | Aug. 21. 2.46 | 0.26 331 | 
Table 3. Ratios of Component Tides. month period. Further adjustment necessiated 
Division | Direction) S»/M,z ky /Mz O;/My by some differences in time-length of the 
l x ] n= a TS running-mean and in the period of component- 
3 . ae Mis 7 ai tides is made at the final step. The amplitude 
7 N 29° W| 0.56 0.62 0.86 and phase of component-tide and amplitude 
8 |N6I°E | 0.23 0.32 0.64 ; ; Rees 
9 | Vertical | 0188 | 0.99 | 0.82 Tatio S,/M;, K,/M;, and O./M, thus o 
10 Vertical 0.92 0.74 are shown in Table 2 and Table 3. 
Table 4. 
nnn ee 
§ €& | Horizontal Areal Strain) Cubical Dilatation | 
@ 5 Amplitude Phase Amplitude Phase v h-3l CG 
3 o in 10-8 lag in 10-8 lag 
2-6 | Ms 2.47 20° 1.91 20° ~0.22 0.47 0.73 
| 
2-6 | Ss 0.80 36° | | 0.59 
7-10 M, 1.94 10° spel 149 —0.30 0.37 0.45 
7-10 So 0.87 Slay” 0.42 iW —0.75 0.64 0.62 
7-10 Ky 1.04 20° 0.39 38° —0.38 0.29 0.22 
(2%) O7 Pe Do) Bos 0.98 22° —0.40 | 0.51 0.64 ie 
Weighted —0.28+0.04 — 0.45+0.03 0.60+0.04 


| 


| 
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Results 

Using observational values, the vertical 
strain, horizontal areal strain and cubical 
dilatation were obtained, the constants for 
which v, h—3/ and C=v+2h—6l were derived 
and shown in Table 4. 

Further, using observations of six com- 
ponents of extensions, the six strain components 
of M,-tide is obtained as follows. 

€rr =0.64 x 10-8 cos(2#—193°.0) 
€99=0.96 x 10-8 cos(2t—29° .6) 
€ye= 1.32 x 10-8 cos(2#—5°.5) 
€re=0.46 x 10-8 cos(2t—226° 4) 
€rg=0.54 x 10-8 cos(2#—183°.3) 
€oy=1.65 x 10-® cos(2#—74°.8) 

In Fig. 1. (a), (b) and (c), are shown the 
relation between the azimuth and amplitude 


of M,-tide and the relations between the 
amplitude and elevation along the sections 


Earth's center 


Bigvedk. 


Izuo OZAWA 


se 
x08 
Ww De Sepals 
S 


wosxio? 


ra ees 


Earth center 


1.0 05 xio® 


iN) 


Earth's center 


Rigen 


of north-south and east-west directions res- 
pectively. In Fig. 2, is shown the relation 
between amplitude of M,-component of exten- 
sion and azimuth or elevations, calculated by 
using theoretical values of F(E), G(é), dF(E)/dé, 
and dG(&)/d—, where — is v/a, obtained by H. 
TaxEucH!. According to Fig. 1-(a), it seems 
that there is a maximum in amplitude in the 
direction of south-east to north-west. As 
there is not only the direct influence of a 
heavenly body but also the influence of oceanic 
tides, so when we want to calculate the value 
of 7 by using @ey or €ea—€yy, it is necessary 
to exclude thoroughly the influence of oceanic 
tides. But by using the Boussines@ solution, 
it can be shown that the influence of oceanic- 
tide upon the horizontal areal strain is very 
small. Thus we can obtain the direct influence 
of heavenly bodies from the horizontal areal 
strain. Hence, 


h—31=0.45+0.03 


as was obtained in Table 4. 

Thus for h=0.60, we have 7=0.05 and for 
?=0.04, which was obtained by K. Sassa~ and 
others, we obtain h=0.57. The mean of 


( 
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amplitude ratios S./M,, K,/M, and O,/M, are 
0.26, 0.68 and 0.57 respectively in the present 
observation against the ratios 0.47, 0.63 and 
0.58 in tide-generating force. 

Conclusion: We have observed the tidal 
variation of six components of extension at 
Osakayama, and have obtained their vertical, 
horizontal areal strain and cubical dilatation, 
from which we got v=—0.28-+0.04, h—3/ 
=0.45+0.03 and C=0.60+0.04. The in- 
fluence of oceanic tide in a place like Osaka- 
yama which is about 65km. away from the 
nearest ocean has not yet ascertained. But 
as the horizontal areal strain is scarcely in- 
) fluenced by the oceanic tides the observed 
constant h—3/=0.45+0.03 is considered to give 
an accurate value for the primary-effect of 
the earth tide. This observed value is in good 
agreement with the value derived theoreti- 
cally by Dr. H. Taxerucut. 

In conclusion, the author wishes to express 
his sincere thanks to Prof. Kenzo Sassa and 
Prof. Eiichi Nisutmura of Kyoto University 
for their valuable instructions throughout this 
research. 
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Time characteristics of magneto-hydrodynamical dynamos 


By 


Shoshichiro Kato 
Geophysical Institute, Faculty of Science, Tokyo University: 


Summary 


In connection with the self-exciting dynamo problem, we study here time characteristics 
of magneto-hydrodynamical dynamos driven by fluid motions of magnitudes larger or 


smaller than those for stationary dynamos. 
numbers of dynamos of different time characteristics. 


For a given fluid velocity, there are infinite 


The fundamental dynamo (P, solution 


in section 4) is developing or decaying according as the fluid velocity is larger or smaller 


than that for the stationary dynamo. 
time characteristics. 
variable stars. 


$1. 


By many studies on the dynamo theory of 
the earth’s magnetic fields using magneto- 
hydrodynamics (Ensasser: 1946, 1947; Butuarp: 
1949; BuLtLarp and GELLMAN: 1954; TaKEucui 
and Suimazu: 1952, 1953; Taxerucut: 1956), 
we have learned. roughly characteristics of 
the dynamo actions in the earth’s core. But 
almost all studies up to the present time are 
on the stationary magnetic fields and little 
have been made on non-stationary cases. Thus 
the questions on the secular variations or the 
reversals of the earth’s main magnetic field 
have not yet been settled. In the present 
paper, we shall make a study on the time 
characteristics of the dynamo actions, or in 
other words, we shall make a study how the 
magnetic fields coupled with fluid motions 
change with time. 

Before entering theoretical calculations, we 
have to consider some observational facts. 
There are some interesting facts in the earth’s 
magnetism, such as the westward drift of 
magnetic fields, slow variations in the dipole 
field and the reversals of the dipole field in 
geological times. But we cannot directly see 
hydromagnetic phenomena in the earth’s core 
which are producing the above observed facts. 
Therefore, if we want to study the dynamo 
actions more directly, we must turn our at- 
tentions to some stars having magnetic fields. 
Many magnetic variable stars haye recently 


There are also dynamos having quasi-periodical 
This periodical dynamo may be useful as a model for magnetic 


been found. Most of them are varying ir- 
regularly with time but some are found to 
have periodically varying magnetic fields with 
constant amplitudes. Little is known on fluid 
motions in these stars. The sun seems to 
have a very weak general magnetic field, 
order of one gauss, of dipole type. But be- 
cause of observational difficulties, it is not 
certain whether the general magnetic field is 
changing with time or not. There are some 
indications that the sun possesses a very strong 
magnetic field of toroidal type (T.° type) os- 
cillating periodically with 23 year cycle. Thus 
the sunspots appear mostly as a pair of two 
spots or two spot groups, namely a bipolar 
group, and the line connecting the two com- 
ponents of a bipolar group is parallel toa 
circle of heliographic latitude. The observed 
law of magnetic polarity is as follows: the 
western spot or spot group has a fixed polarity 
for all bipolar groups appearing during one 
114 year sunspot cycle, and the eastern group 
has the opposite polarity. This polarity is 
reversed from one spot cycle to the next. 
The law of polarity in the southern hemisphere 
is the exact reverse to that in the north hemi- 
sphere. It is quite certain that the dipole field 
and the toroidal field in the sun are connected 
by fluid motions of toroidal types manifested 
such as the so-called equatorial accelerations. 
Although local fluid motions in the sun are 
turbulent and far from stationary, the above 
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mean motion of equatorial acceleration type 
is remarkably stationary. There is also some 
indication for the existence of stationary fluid 
motion of poloidal type (S.° type) in the sun. 

Turning now to the theoretical sides of the 
problem, we shall at first pay our attention 
to the ScuwarzscHILp’s theory on magnetic 
variable stars (1949). In his theory, he as- 
sumes a strong external magnetic field and 
gets changing magnetic fields and fluid motions 
as some perturbations from this initial state. 
Thus the fluid motions in his theory are 
changing their directions periodically. It is, 
however, difficult to see how such motions 
can subsist in a rapidly rotating system with 
a strong magnetic field. From what was said 
in the last paragraph, it seems better to as- 
sume the existence of some stationary fluid 
motion at first and study time characteristics 
of magnetic fields thus caused. This will be 
our main point of view in the present paper. 
This is also a natural extention of the results 
obtained by Taxkeucut and Suimazu (1952a, b; 
1953) and Butuarp and GELLMAN (1954) on the 
stationary dynamo. 


§2. 


According to the program in section 1, we 
must solve, under certain boundary conditions, 
the following Maxwell’s equations 


OH 


curl H+ at =()), (2.1) 
curl H=47rJ=4rne(E+VxA), (2.2) 
bh TEA, (2.3) 
div H=4z pc*e-1 (2.4) 


for given velocity field V. Customary notations 
are used in the above equations. From (2.1), 
(2.2) and (2.3), we get 


4 4 pat Ane curl (Vx A)=0. 


In spherical coordinates, 
and (2.5) is given by 
H={SH3+SH7r}e"' , (2.6) 
where the (7,9,¢) components of S-type 
(poloidal type) and T-type (toroidal type) 


magnetic fields are given respectively by 
S-type magnetic field 


(2.9) 


HT satisfying (2.3) 


—n(n+1)Si(nr23¥n 
. OYn » 
[_ [r ds” tease ae 


Hy; = dr (250) 
ds” (es Wo OX 
|. E ee * s(nt1se sin 000 
T-type magnetic field 
- 
tien ah OS ances 
Hips) okeia Odum Meee 
OX a 
Tir ag 
in which 
Yn = Fa (cos 9) © cin (2.9) 


On the other hand, V is assumed to be ex- 
pressed as an arbitrary linear combination of 
the following two fields: 

S-type velocity field 


—n(n+1)Vsnlr)r2Yn , 


” Ws Ya 
_j-[ve —+(n-+1) Vm pr nee , 


dies ayn 
rie pas 1 Vs. t—1__ 
-|7 oe alr sin 006 


(2.10) 


with m=2 and 22=2c [2c means cos 2¢ in the 
expression (2.9)]. 
T-type velocity field 
0, 
= nm (r)r” OY = 
Vr= ig sin 00¢ ’ 
OYe 


Vin(r)r” 36. 


(2.11) 


with m=1 and m=0. 

The expressions for the magneto-hydro- 
dynamical couplings of general (poloidal and 
toroidal) magnetic fields by general (poloidal 
and toroidal) velocity fields are given as fol- 
lows. 

With the suffixes 


@:; primary (inducing) magnetic field 
8: velocity field 
7: secondary (induced) magnetic field, 


the orthogonal conditions for Maxwell’s 
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equation (2.5) are expressed by 


afr (2u ‘ 
\ \ \ D| —Areb le + pHa b4re curl (Vex Ha) |- Hy sin Odrd6d¢=0 . (2.12) 


0oJOJ0 @ 


Carrying out the integrations with respect to @ and ¢, we get 


| E SeSy +E Tab, |rdr= 0 
# = = , (C463) 
\ | = SeTy+> TD, |rdr=0 
o| @ @ 
where 
a a (2r 
S8,=| | | —ArebH se +7*Hse +e curl (= Vex Hse) | Hs sin 0d0d¢, 
0 Jo 
; (2.14) 


TaSy= | | | —ArabHre +P?*AHrea+4r« curl (= Vex Hira) | Hy sin 6dédé¢ , 


o Jo B 
ete: 
We shall call Ee, Tags ; hye and Fil, the coupling coefficients. By putting p=0 in (2.14), 
we have the coupling coefficients already discussed in the papers above refered to. Introducing 
the term into (2.14) will not bring any difficulty, so we shall omit to write down the 


detailed forms of the coefficients. The following change of the independent and dependent 
variables will be convenient for future discussions. 


r=ae (a: radius of the earth’s core or a star) 


( n (7) )- 1 ee ) (eee )= ply (Ge ) (2.15) 
Vs.nlr) a Vs, n,mV's.n(E) Vil) a” Vr,n,mVr, n(&) 

where Vs,n,m and Vr,n,m are some constants designating the magnitudes of the corresponding 
velocity fields. In considering boundary conditions, we shall assume an insulator free of 
fluid motions outside r=a. Introducing the time-dependent term into the fundamental equations 
will bring no change in the boundary conditions. This is understood as follows. At first, 


putting «=0 in (2.5) will elliminate the time-dependent term from the equation. Thus, we 
have in the insulator the same magnetic fields 


n(E)= AE ntn , Tn(E)=0, 


where A is some constant, (2.16) 


as in the time-independent case. Thus the boundary conditions will give 


wn py pagan Gon ag “ Ls (2.17 
mE)+ on 14 de 20),., “Te)s0° at €=1. 2.17) 


Assuming a rigid boundary at €=1, we have the following conditions for the velocity fields. 
Vr(E)= Vs(E)=—_, =): Ate ee?, (2.18) 
Considering (2.18), we shall take rather tentatively 


1 
Vs (7) — a Va, 20 V9 (E) , 
(2.19) 


1 
Vo, (n=, Vr oViE) ) 
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2 = 0<&<1 
—&) it E) re lsc (2.20) 


2e(E)= k V,! A 
Vi%(6)= {3 ae 4 


where V»s,.- and V;,9 are constants designating the magnitudes of the corresponding velocity 
fields. 


§ 3. 

As in the case of the self-exciting dynamo, we are considering a group of magnetic fields 
which are coupled by 7° and S,”” motions and in which we can find a dipole (S,°) type 
magnetic field. As is easily seen, there are infinite numbers of magnetic fields in the system 
we are considering. It is almost prohibitive to take all of these magnetic fields into con- 
sideration, so we shall confine ourselves to the case when there are magnetic fields of 
n,m<2 only. In this case, we have 


1 coe cae as ee 
| Es = Pie's. = ViPpeSs o ar TS» |prae =0 , 
0 


ye (3.1) 
1 aa “aa Taga ZN 
| | Ser 4 FT O4 FTO + Tete lpg ath: ete 
0 
where 
a 16/dS° 4 dS g*1a @S, 4 dS dS,° 
S,9S)o=- ee \sy 10 rues le 21 pee” “\K 1 °) 
wing geet gery Mghe ae gee aR Mees 
4 dS 
(OR DE ) 
3 [ ane (ae +25,)'], 
T,S=0 , 
ae 48x24 /0 48x12 /0 ld dS,° 
T7285 9= —- 2V7,2¢ T2085 0 ma 2 a 47, 2¢T 20,8 1 
1 P PNG 27° T7)*S, - es Ye gee Ye T; (Fae + 2S,° 
SPT) =—"P aS ae ET,o , 
ey Te tad WARE W Be a eae 
Ti0T Oat 2 2 \¢s po_** peep. oy2 
PT rere Each Gnas ATR eee 


“90, 87 288 x6 / 1° d d 
T2287 9 = — y4 . 0.8) a : - 
35 \O 41 GE -(E4V,% 77° ) + gg EVs" JET? ery, 


So, 3.1- 98 Ee v4 Co {servis : (aes) ae ESV/,2)S,0 ]enee. 


a OF 288x6 /1% 
DEAT fe I 1V,2°7,0 
: 35 0 DY ms )+ dé 


0,87 20 exe pUT 8" | dT: 
1h AS (TA en 82a 
Es. 


ee VETS ET, 


\0 
20, 8 eee 48 x6 pay ey! 1T’,2¢ 
ty T,8= 5 Gi we got oe "ET" — tke Sab TeV OT 
0\. 6x48 


— 47, 2¢,8 728 
7 as ae aa ae (3.2) 
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Gb ed ne Ss 
In (3.2), Sap in the expression for 7.2%%S,°, for example, means 0 and 1 for T,2° and 


T,"*, respectively, and z, y and P are non-dimensional parameters given by 
x=4rKaV,,9, y=4rKaVo,. and P=4rrKap . (333) 


From (3.2) and (3.3), we have the following simultaneous differential equations for the radial 
functions S,°, T,°, T.?° and T,?*. 


@s,° 4 Sv 9 x 24 
—PES,° ee 1 
ES +E— dé bas 5 
Oh Aes if hel 12 dV1°¢ 4 72 Bia 
: i 


yee V.% f ’ 


ic eke an SEF Va) a + ge EViMDETS | 


be arte Pde Moe de 
TT,” ,dT:* 6 
==) 2) Iho = a 6 2 SSG A GE 47,2¢ i) (E38 |7. 26 0|—97V,9 28 
ET MEL de? = re Ae V2"T> nO Us ers | 2eV PET» 
he eT. | dT 2 d(dSp , 2 d 
— PET, + +6 = 22 VET, eae 3 V2 Me ibe 37,20 
ae GR ey eet ela He a rae +aelae as |. 


(3.4) 
It is needless to say that by putting P=0 in (3.2) and (3.4), we have the results already 
obtained by Taxerucui and Suimazu. As is seen from (2.17) and (3.4), we have here an eigen 
value problem. Thus for given x and y, we get the eigen values P which determine the 
time characteristics of the magnetic fields concerned. 


§ 4. 

In solving the eigen-value problem defined in the last section, we shall use the approximate 
method devised by Taxkrucui and Suimazu in solving the self-exciting dynamo problem. Thus, 
taking the boundary conditions (2.17) into consideration, we shall put 


Si)= A 14 b) + ARIF) ts + Acee(1— Ege, (4.1) 
T(E) By ( B, By, 
| rac }-0-»] ake Eee pon] (4.2) 
‘ T.5(E) Dy D, Dn, 
where Ay, Ai,--:: are constants to be determined. Each radial function under each constant 


satisfies the boundary condition separately. Inserting (4.1) and (4.2) into (3.1), we get the 
equations of the following types: 

195 8 24 72 
ay ee SAS! ae yDy— ee tc. : 
( 6 P)As+( ; =P )A— Se yDy— Le yDrt ++ =0, ete (4.3) 


In Table I, are shown the coefficients of the equations when we assign two undetermined 
constants to each radial function. 
In order that these equations are compatible, the determinant formed by the coefficients 
must be equal to zero. When we assign one undetermined constant to each radial function, 
we can reduce the determinant equation to 
P!+-P3 x 1.17473684 x 102+.P?(4.91115789 x 103 +-0.362? +4.82153483y") 

+ P(8.34896842 x 10+ 16.370526312? + 3.55674005 x 1074?) 

+ (4.42004210 x 10°+1.22778947 x 102x? —0.18602870827y? +6.55555504 x 10%y? 

—1.65251476y*)=0 (4.4) 
Similarly, when we assingn two undetermined constants to each radial function, we have 
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P®+ P’-4.47142823 x 10?-+.P®(8.23885278-+ 1.05050503 x 10-42?+-0.297129360 x 10-2y2) x 10! 
+P*(8.11219451 +3.39124570 x 10-*x? +1.01348934 x 10-74?) x 106 
+P"(4.64648145 +-4.36486426 x 10-42? +-5.75305086 x 10-822y2-+0.132231402 x 10-824 

+1.34904708 x 10-y?+0.372336199 x 10-8y*) x 108 
+P9(1.58551541 + 2.82720382 x 10-4a?-+ 0.129065369 x 10-8a2y2+4-0.241164891 x 10-824 

+0.890744496 x 10-?y?+-0.978920402 x 10-8y*) x 101° 
+P(0.315821042 + 0.942943824 x 10-4a? +9.31028155 x 10-8x2y?—0.961831261 x LO™ 2G 

+0.136209793 x 10-82 —0.759558608 x 10-1aty? 

+0.306994192 x 10-2y?+-0.885702498 x 10-%y4 —2.87105340 x 10-124) x 1012 
+P(3.33911212 x 10-?+-0.146822109 x 10-42? +2.27749165 x 10-82?y?—2.01435413 x 10-1224! 


+:2.85018778 x 10-2 —3.31001411 x 10-Maty? 


+5.23901081 x 10-*y?+0.319302975 x 10-®y*—6.31631807 x 10-12y) x 1014 

+ (0.137213997 x 10-?+-0.750857133 x 10-%x?-+0.133079886 x 10-8a2y?— 1.04327783 x 10-!2x2y4 
+0.166611568 x 10-12*—2.30403529 x 10-Maty? 
+0.345266759 x 10-*y?+4.01066177 x 10-8yt—3.47397522 x 10-12y®) x 1018 


=(). 


For given y and 


xz 

ia te (4.6) 
we solve these equations with respect to P. 
The results obtained from (4.4) and (4.5) are 
shown in Table II and III respectively. 
From (4.4), we get two real roots P; and Pr 
and one pair of complex roots Py. Similarly, 
from (4.5), we get five groups of roots Py, 
Pua, Pus, Pma and Pip. Pr root in 
Table II or in the first approximation seems 
to separate into Pya and Pyg in Table III 
or in the second approximation. Similarly, Pr 
and Py in the first approximation seem to 
correspond to P; and (Pma and Pip) in the 
second approximation. It is needless to say 
that we have one independent group of radial 
fOMetIONS 9)", 222, Le° and Z.7? for each P 
and thus for given w=4zKaV;,. and y 
=4rKaV2,., we have many independent so- 
lutions of the problem. 

We shall now consider time characteristics 
of the independent solutions corresponding to 
IP, Jex ehatel Jeanne 

Py,-solution: The value of P in this series 
is real, so the corresponding magnetic fields 
are developing, stationary and decaying ac- 
acording as P20. As was shown by TakeE- 


10 40 


100 40 


30 | 


1000 | 40 
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(Gila) 
Table II. 
P, Pr Pr 

= 93 %43 —46.84 | —23.61+4 44.85% 
LH ww —52.56  —23.59+ 69.06% 
1.99 —58.31  —23.59+4 92.90% 
— 6.236  -64.07 | —23.58+116.6 i 


— 0.4761 —69.83 | —23.58+140.2 


— 6.303 


42.64  -34,274128.2 


— 0.5315 | —48.47 | —34.244192.6 


5.798 


12,3894 


19.05 


— 3.215 


2.543 

8.885 
15.56 
22.39 


3.124 
2.580 
8.933 
15.58 
22.39 


—54.78 | —34.234256.9 
—61.38 | —34.224321.2 
—68.08 | —34.22+4385.5 


~42.30 —35.98-+1200 
—48.06 —35.98-+1801 
—54.40  —35.94+2402 
—61.05 —35.99+3021 
—67.86 | —35.98+3603 


~ 42.30 —36.00-+ 12000 
~48.06 —36.00-+18000 
—54.41 —36.00+24000 
—~61.05 —36.00+30000 
— 67.86  —36.00-4 36000 


a 
a 
a 
a 


oa a7 en er] 


°. 


a 
a 
a 
D 
a 


a 


ucHI and Suimazu, for each a=a/y in this 
series, there is one value of ¥, Ya,o, Say, for 
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Table III. 
a a 
a y Pi Pra Pir Pia Pure 
W209) 223.8 ~37.03 {= 96-72) | 42.414 18.25) | 52.684 97.891 
| 40} —18.59 | — 46.22 va | 
A 60 | 249.3710 — 54.06 ‘- 99/98} | ~39-61+ 72.41i —53.874321.4 7 
80 | — 6.106 — 62.86 ao | 
100 | — 0.1295 | ~75.03 lran's } —39.62+123.3 i —53.944539.2 7 
| 20| — 6.838 | —45.09+ 5.613) | — 98.82 | -66.97+ 85.68% —58.69+209.3 i 
40) 2.032 | = 98.26 - | 
10 60 11.06 =—55 6525 .15.4malere Olde —65.474244.6 i —59.284650.3 t 
80 | 19.68 = 06.42 oa 
100 | 28 20 —65.43436.79 i ~ 95.14 | —65.344406.2 4 —59.3341087 7 
| = 
OO. 3.295 |) — 40. 2112.21 4 — 98.87 —71.704766.2 4 —60.58+1904 4 
40 6.148 | — 98.29 
100 60 15.28 —50.254+28.90i — 97.39 —71.69+2296 i —60.594+5714 1 
| 
Peso} 923-41 4 — 96.16 
100) 32.18 | —59.88441.22i | — 95.00 —71.6943827 i ~60.59+9524 7 
ie 3231 ~40.05+12.66 i | ~ 98.76 | —71.77+7648 « | ~60.61-+19020 % 
40 6.207 | — 98.29 | 
1000 60 15.34 —50.17+28.97 4 | — 97.38 | -—71.77422940 i —60.61457060 i 
80 23.99 | ~ 96.23 
100 32.22 —59.83441.24 4 | — 95.00 | —71.77438240 i — 60.61 +95090. 7 


which P=0, that is, the corresponding magnetic 
fields are stationary. Refering to the above 
Ya,o, We can say that the magnetic fields 
in this series are developing, stationary and 
decaying according as ya = Ya,o. One physical 
interpretation of this result is as follows. The 
magnetic fields are stationary when the Joule 
dissipation of the fields are compensated by 
the work done by the fluid motions. y is a 
meassure of the work against the Joule dis- 
sipation. We may expect that when the work 
done is larger or smaller than just to com- 
pensate the Joule dissipation, the correspond- 
ing fields are developing or decaying. Thus 
we get the above results. 

Py, -solution: The value of P in this series 
is real and negative in the first approximation. 
In the second approximation, we have two 
roots, one is real and negative and the other 
is complex. The real parts of these Py are 
always negative, so the corresponding fields 
are decaying even when we have developing 


fields corresponding to Py. For complex roots, 
we have some oscillating fields. But as the 
damping of the fields are so large that we 
can hardly see their oscillating characteristics. 

Pyy-solution: The real part of P in this 
series is always negative, so the correspond- 
ing fields are decaying. For large a, the real 
part is almost constant whereas the imaginary 
part of P becomes larger for the larger ay. 
Thus we have 


real part of P<imaginary part of P 


when ay or «>1. (4.7) 


As the imaginary part of Pya,p is pro- 
portional to the frequency of the oscillating 
fields, we have here quasi-periodical fields, 
whose damping during one period is very 
small. The frequency of Pyyp-solution is 
about twice that of Pyra-solution. Py may 
be an over-tone of Pyra. Anyway, the Pyy- 
solution we have here may be useful as a 
model for magnetic variable stars. 
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In the quasi-periodical case, there are some 
phase lags among the four radial functions S,°, 
T,°, T,?* and T,?5. Since we have the terms 
as —PES,° (P is almost pure imaginary) in 
(3.4), we cannot satisfy (3.4) by using only 
the radial functions of cos p¢ (p is pure real) 
types. We must introduce the radial functions 
of sin pt types to satisfy the equations, thus 
we have the phase lags above refered to. 
Refering to ParKER’s work (1955) on dynamo 
waves, ELSASsER (1956) says that there may be 
a model for magnetic variable stars in which 
energy exchange is made between two mutual- 
ly orthogonal types of modes of the magnetic 
field, poloidal and toroidal, rather than between 
magnetic field and fluid velocity field as in 
the ScHwarzscHILp’s model. The result obta- 
ined here corresponds to his remarks. 


§5 

In this section, we shall apply the results 
obtained in the last section to the interpre- 
tation of magnetic phenomena in the earth 
and the sun. Refering to (2.6) and (3.3), we 
define the decay time rt and period T of the 
fields as follows. 


t=1/real part of p=4zxa?/real part of P, 
ket) 
T=2z/imaginary part of p 
=87?xa?/imaginary part of P. (5.2) 


For fixed P, both the decay time and period 
are proportional to «a?. For the earth’s core, 
we take a=3x10%cm and «=3x10-° e.m.u. 
a@ is the ratio of the magnitude of the toroidal 
fluid velocity to that of the poloidal velocity 
and is considered to be large (about 100) in 
the earth’s core. Thus from the results in 
Table II and III, we have the decay time and 
period of several hundred to several thousand 
years. It is interesting that these time scales 
are of the same order of magnitude as those 
in the westward drift of the earth’s magnetic 
fields, etc. 

The radius of the sun is 7x 10®°km=7 x10" 
cm and according to T.G. Cowuine (1953), 
electrical conductivity in the sun is between 
10-8 e.m.u. (at the surface) and 10~* e.m.u. 
(at the center). 

Corresponding to k=10-° aiatal D4 paar 


we get from Table II and III decay times of 
the order of 10° and 10° years. These time 
scales are much larger than those in the sun. 
We shall now assume that the (quasi-) 
period in the Pyy,q-solution is corresponding 
to the 23 year period in solar activity. In 
order that this may be so, we must have 


a(=4rKaV,,9) ~ O(10" ~ 108), 
when kK=10-* e.m.u. (73) 


The time ¢, say, during which the toroidal 
flow makes a complete round around the 
solar axis is given by 

2na _ 8x*Ka? 


h= Vi: ; ive (sec.) . (5.4) 


> 


Putting a=/xl0"cm, £=10- emu sand 
(5.3) into (5.4), we get 


600 day 107 
ty & = 
; | Sine ‘ioe OD 


If we take «=10-* e.m.u., we have ¢ which 
is 10* times that in (5.5). In the so-called 
equatorial acceleration in the sun, fluid parti- 
cles near the equator make one complete 
round relative to those near the poles in about 
130 days, which agrees well with (5.5). Since 
the depth of the convective layer in the sun 
is small (about 1/10) compared to the radius, 
and 7) in (5.4) is proportional to @?, we may 
be able to adjust ¢) so as to agree the observ- 
ed value even for a larger x. 
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Summary 


Wave generations from line sources within the ground are studied. 


For a compres- 


sional origin of step function type in time, we calculate displacements at several points 


within the ground. 


In contrast to the surface displacement, we have, deep in the ground, 


three distinct pulses corresponding to the direct P, reflected P and S phases and no dis- 
placement corresponding to the RAYLEIGH wave. 
Next, we calculate profiles of surface displacements for an impulsive compressional 


origin. 


§1. In a previous paper (H. Takeucut and N. 
KoBAYASHI: 1955) titled wave generations from 
line sources within the ground, we calculated 
ground displacements to be obtained at the 
free sourface bounding the ground. The 
source was assumed to be of step function 
type in time. The noted results thus obtain- 
ed are as follows. At the epicenter, the dis- 
placement is a simple pulse followed by a 
gradual decrease to a permanent displace- 
ment. At a certain distance from the epi- 
center, RAYLEIGH waves appear. Thus the 
displacement at a distant point is composed 
of a P pulse, a S pulse (for S type origin 
only) and a movement of relatively low fre- 
quency corresponding to the RAYLEIGH wave 
which tends asymptotically to a permanent 
displacement. As pointed out by C. H. Dix 
(1955), to study the generation mechanisms 
of this long wave is of theoretical and 
practical importance. In the present paper, 
in order to make this point clear, we shall at 
first calculate displacements at several points 
within the ground. Secondly, we shall calcu- 
late profiles of surface displacements at certain 
fixed times. In the first calculation, we shall 
take the source of step function type time 
yariation as in the previous paper, whereas 


The importance of surface reflections in generating long waves is shown. 


in the second calculation we shall take the 
source of impulsive type. 

§2. Refering to the rectangular coordinates 
in Fig. 1, we take a compressional line source 


Ox ; 


ge se e'? HR) 


Uo 


= Wes cosh $ ) Jp, 
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where 
wars trfyten?, OD 


at a depth f below the free surface. The 
corresponding experessions for % (2, y, 2) and 
v (x, y, t) will be 

+ co 


E Ee e-#(F+4) 
Pays 


u=telrt 
| : 


aca (f-¥) | 
tes Si =e a 
- (as 4 Zi 


_9 PEOET— ED 01-64 |e dé , 


v=o" | ve ees a e-uS+0) 
ae 2 


at 2 eer) o-0f-By } elt dE , 


where 
F(€)=(2&?—k?)? 4a RE? , 
F (E)= (2&2—k*)?+4aBe , 
a’= &—h?, B? = F2— k?, 


p 

k=py/ a a 
As is seen below, the first, second and third 
terms in w and v in (2.2)lrepresent respective- 
ly the movements beginning at the arrival 
times of the direct P, reflected P and reflect- 
ed S. In order to get the expressions of pro- 
gressive waves at woo, we must have 


(2.2) 


the real parts of a and B>0. (2.3) 


We can satisfy these conditions by interpret- 
ing the integral expression for wz in (2.2), for 
example, as the limit of the contour integral 


—7pipt a, 1 iq = 
yi € FO) -Ve=F + 
2VC—h F(f) 


VOC —K) 
F() 


-VG—Taf - roe Y 


-/@-FA(f-y) 


Xue ets dl, 


F)=00= BY AOV OB/OB, 


FO)=20?—k)? 40° C/O, 
(2.4) 


when semi-circles around six singular points 
+h, +k and -:« are made vanishingly small 
and /C?—h? and 1/€?—R? at the right end 
of LZ are assumed to be positive real. In 


Pigs 2. 


Fig. 2, +:« are the poles of the integrands in 
(2.4) for which 


F(+«)=0. (2.5) 


The first and second intergrals in (2.4) are of 
the same forms as those we had in the pre- 
vious paper and may be calculated as there. 
The third integral which is of the form 
—-V¥C—-hf-VC —ky+ i¢x+ipt dt, 

(2.6) 
where <€> is some algebraic function of €, 
is the one we have for the first time. In 
order to calculate 7 by the method of contour 
integral, we shall integrate the integrand in 
(2.6) on the path L+(B-F—D->A-E->G-C), 
where BF and GC are parts of a large circle 
with the center at the orgin O and D>A->E 
is the path on which 


=-YO—} f—V CR y+ ita (2.7) 


in (2.6) is pure imaginary, and on FDor EG, 
the imaginary part of Z is kept the same as 
that at Dor E. The points D and E are 
made far away in the following analysis. 
Since’ there is no singular point inside the 


1=| <e>e 
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path and since the integrand vanish on BF 
and GC, we may put J equal to the integral 
on G>E>A-—D-F. As in the previous 
paper, there is some reason in omitting the 
contributions from the paths G>E and D->F. 
Thus we have 


1=| <a eS eAtrint dt z 


§ 3. In arder to transform (2.8) into more con- 
venient form for numerical calculations, we 
shall make some mathematical considerations 
in the-present section. 

(a) € at A, &o, say, is real and negative. 
C, Ve?—h? and 1 C2—k*? on “AD are of the 
forms —(-+)+2(+), where (+)’s are some 
positive real numbers. € on AE is the con- 
jugate complex of that on AD, and /C?—h? 
and '1//€?—k? on AE are —(conjugate com- 
plexes) of those on AD, respectively. 

(b) Zon AD and AE can be put 


(2.8) 


Z=—tlyX, UaXK— oo, (3.1) 
where 


hoa VevaPs pan Sra) 


It must be remembered that 1/€)?—h? and 
VC,2—k? are positive imaginary and —€, is 
positive real and thus € is (3.2) is positive 
real. 


(c) Taking the differential of (3.1), we have 
f<C>dl=6, dX, (3.3) 
where 
ant BAGS ee Se 
VOB VOB ff: 
: (3.4) 


Inserting the above results into (2.8), we have 


LB 
} AD A | i ae 
From what was said in (a), it can easily be seen 
that <€> and <C> on AE are conjugate 
complexes of those on AD, Thus, we haye 


i(pt—¢oX) é, dX. 


(3.5) 


nase i(pt- GX) 


? 
where V(€) is the imaginary part of <¢€> 


l<t>. 
(d) Putting 


I | VOQe b°dx, (3.6) 


(337) 

| 8 
and assuming iS (¢) instead of e’”* in (2.1) for 
the time variations of displacement and stress 
near the origin, we have 


+, 2a 
if 
By VEIT py VET yer 


I [-viens@—absx) éy dX, 


(3.8) 

S’ means the differentiation with respect to 
the argument in the bracket. 

(e) The radial displacement wp and radial 

stress Terr in the immediate neighbourhood of 
the wave origin are given by 


eo) 


2 
aa Tre=——_ 9). 3 
R? RE= — pe) (3.9) 
In the following numerical calculations, we 


snall assume 


S()=0, —o<?t#<g; S(é)=const, 720 
or S’‘()=0, t0; SW=o, S’Qdti=s, 
isayauen— 0: (3.10) 


(f) For the time variation in (3.10), we 
have 


2S 
T=—-ViG), 
IAGO) 


where €’ is given by 


igen ae ee 
eo our ae 


iV 7m athe . (3.12) 


(3.11) 


§ 4. Practical numerical calculations are made 


as follows, At first, giving sy x/f and 
p 
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S- 


Fig. 3. 


Cala 1073 


y/f, we solve (3.12) and get €’, the real and 
imaginary parts of which are negative and 
positive, respectively. We next calculate 
V €72—] and VY 2’2—m?, the real and imaginary 
parts of which are negative and positive, 
respectively. The selections of the signs 
above made are according to what was said 
in (a) in §3. After getting €’, Ye72—] and 
V €’/2—m?, we calculate 


V(€’)=imaginary part of <¢’>/<¢’> 


where <¢’> and <€’> in the present case 
are by (3.4), (2.8), (2.4) and (3.7) 


aie seas 


, (4.1) 


OOS 00m OV CV Cm 
(4.2) 

, C’ i pen ee 
<= ps tiveicm pt 


Other integrals in (2.2) are calculated in the 
similar way. All of these integrals diverge 
as f>oo or |€’|->0o. But it can be shown 
that the three terms in the expression of u 
and v in (2.2) adjust themselves to give finite 
wu and v. CaaniaArp has pointed out this 
divergence effect (CaaniaRp, 1939, p 128-133, 
see also Dix above cited, p 98) and said that 
this is an indication that the potentials are 


not physically significant quantities. In our 
numerical examples, we assumed 
eng Me aay 3 AT (4A) 
a Vs 


and calculated the displacements at #=0, 


y=f/2, 2f and 10f; «=10f, y=/f/2 and 
10f. wu and v thus calculated are shown in 
Fig. 3. The results to be noted are as follows: 

(a) As it should be, there is no displace- 
ment before the time t=7/V», which is the 
time of transit of P wave across the distance 
rag? + f2. 

(b) Horizontal displacement w on the verti- 
cal line =0 is zero all through the time. 
This will be understood by the consideration of 
symmetry as shown in Fig. 4. As was shown 


Fig. 4. 


in the previous paper, vertical displacement v 
at («=0, y=0) was a sudden displacement 
corresponding to the direct P phase followed 
by a gradual decrease to a permant displace- 
ment. Vertical displacement at («=0, y<0) 
is composed of a direct P and a reflected P 
pulse followed by a gradual decrease to a per- 
manent displacement. As it should be, the 
direct P pulse is down or upward according 
as y2f. Vertical displacements on the line 
2=( have recently been calculated by Drx, but 
in his figures only the correction terms to the 
perfect image point reflection are shown, 
whereas in our Fig. 3, the total displacements 
are shown. 

(c) As was obtained in the preceeding paper, 
each component of the displacement at (7=10/, 
y=0) is composed of a P pulse and a move- 
ment of relatively low frequency corresponding 
to the RayLEIGH wave. Displacements at (v= 
10f, y= f/2) are similar to those at («=10f, 
y=0), but they have two additional pulses 
corresponding to the reflected P and S phases 
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in the neighbourhood of the initial motion. 
At (x=10f, y=10/S), these reflected pulses are 
well. separated from the initial pulses. . Thus 
we have here three distinct pulses and no 
displacement corresponding to the RAYLEIGH 
wave. 


§5. In order to understand generation mecha- 
nisms of RAYLEIGH waves, we shall calculate 
profiles of surface displacements at certain 
fixed times for the origin of impulsive type, 
for which S(#) in (3.9) is 


Se VSUy SQ=0, S(d)di=S, 
say at ¢=0. 


i): 
(5.1) 


The displacements for the origin of impulsive 
type are obtained by differentiating those for 
the origin of step function type with respect 
to time ¢. Profiles of the surface displace- 
ments thus obtained for ¢ /- L=2, 5 and 10 
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Fig. ‘5. (c) 


are shown in Fig. 5. At ¢#=2f/V», which 
is the time when the wave front is at 
the distance twice the depth. of the origin, 
horizontal displacement # is toward the 
epicenter and vetical displacement v is 
downward everywhere on the free surface 
except in the immediate neighbourhood of 
x=V 3 7, where the wave front cuts the free 
surface. At #=5f/V», there appear a region 
of. upward. vertical displacement and. two 
regions of large positive (away from the 
epicenter) and small negative (toward the 
epicenter) horizontal displacement nearer to 
the epicenter. At ¢=10f/V», the above cited 
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Fig: 6. 
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regions expand in their dimensions. Thus we 
may say that the characteristic profiles of 
surface displacements are formed between 
t=2f/V» and 5f/V», probably nearer to the 
former. The time distance curve for the point 
of zero vertical displacement bounding the 
regions of positive and negative vertical dis- 
placements is shown by (a) in Fig. 6. The 
curve gives the velocity corresponding to the 
RAYLEIGH wave. In Fig. 6, are shown also 
the time distance curves (b) and (c) for the 
points of zero horizontal displacements bound- 
ing the region of large positive horizontal 
displacement. (d) is the curve for the direct 
Pwave. As time goes on, the distance between 
the two points corresponding to (b) and (c) 
becomes larger at first to approach a final 
distance of about 1.57. We may safely con- 
sider that the final distance corresponds to 
a half wave length of the RAYLEIGH wave. 
Thus the wave length of the RAYLEIGH wave 
formed is about three times the depth of the 
wave origin. In practical cases, the wave 
length will be determined by the type of time 
variation of the displacement etc. at the origin 
rather than by the depth of the origin. The 
above result is, however, useful in considering 
the effect of the origin depth on the wave 
length of the long wave formed. 
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Fig. 7. 


The (vertical) displacement at the epicenter 
is shown in Fig. 7. The initial motion is, of 


cource, upward, but its duration time is very 
short. It is omitted in Fig. 7. Next comes the 
time interval of downward displacement fol- 
lowed by the time of upward displacement. 
The time interval of downward displacement 
ends at a time a little later than 2//V». Thus we 


Fig. 8. 


see that the time interval of downward dis- 
placement at the epicenter corresponds to that 
of surface profile as shown for =2//V> in Fig. 
5. This is the reason why we put the time 
when the characteristic surface profile such 
as for t=5f/V>» in Fig. 5 is formed at a time 
a little later than 2f/Vp. Anyway, we have 
the schematic representations of surface pro- 
files as shown in Fig. 8. It is interesting 
that the above results agree well with those 
obtained by Drx in his study on the mechanism 
of generation of long waves from explosions. 
In his paper (p 88), he says 

(a) The spreading spherical pulses from the 
buried source is converted into non-spherical 
secondary disturbances on striking the free 
Surlace: si.. 

(b) An appreciable amount of this energy 
goes back toward the source or toward a 
vertical axis through the source. Thus energy 
is held near the source and near the free 
surface for a time long compared with the 
time it takes for energy of the signal to be 
transmitted away from the source neighbour- 
hood. 
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(c) This long held energy spreads outward Dix, C. H.: 
in the available modes, one of which is the 1955 “The mechanism of generation of long 
surface wave mode of RaAyIgIGH. waves from explosions” Geophysics, 20, 


It is not difficult to see that the time inter- 87. 


vals (a), (b) and (c) above correspond to our T#euchi, H. and Kobayashi, N. > 
1956 “Wave generations from line sources 


a), (b i ig. 8. 
fee PARRA ea within the ground” J. Phys. Earth, 3, 7. 
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Rayleigh and Sezawa Waves Generated by Explosions. 


By 


Akira KuspoTera 


Geophysical Institute, Faculty of Science, Kyoto University, Kyoto. 


Summary 


Characteristics of ground motions caused by small dynamite charges were studied. 
Seismometers were set up at various depths within bore-holes, and varying the depth of 


shot-hole, the explosive waves were recorded. 


In this experiment, surface waves which could be considered as RAYLEIGH and SEZAWA 


waves were observed. 


The characteristics of these waves are as follows. 


RAYLEIGH wave: 
The wave velocity is 60 m/sec. 


The period of motion is about 0.22 sec. 
The amplitude decreases exponentially with depth. 


SEZAWA wave: 
The wave velocity is 500 m/sec. 


The period of motion is about 0.1~0.08 sec. 
As to the vibration mode, there is one node in the upper layer. 
These results are in good agreement with those predicted by the theories of RAYLEIGH 


and SEZAWA waves. 


Introduction 


Putting seismometers at various depths 
within bore-holes, and varying the depth of 
shot-hole, explosive waves were recorded. 

The most interesting result obtained in these 
experiments was that the RayneicH and Ss- 
ZAWA waves could be produced by explosions, 
and their amplitude distribution with depth 
could be observed in detail. 

The experiments of the same kind have 
been carried out by the writer (1955). At that 
time, interface waves, which could be con- 
sidered as STONELEY waves, were observed in 
the interface between sand and clay. 

There is a detailed theoretical investigation 
by K. Sezawa and K. Kanat (1935 a) on disper- 
sive Rayleigh-type waves in a superficial layer 
on a semi-infinite medium. According to the 
study, there are two kinds of dispersion curves, 
one corresponding to the RayneicH (M,-) 
wave and the other to the Srezawa (M.-) 
wave. The particle orbits and the amplitude 
distributions with depth of these waves have 
been calculated. 


On the other hand, it is known experiment- 
ally that there exist M,-waves in the explosion 
waves. Studies have been made on them by 
B. F. Howe tt (1949), M.B. Dosrin (1951) and 
A. Kusotera (1951). 

On M,-waves, K. Tazime (1956) and H. 
Suma (1956) have recently made experimental 
studies as the members of Seismic Explora- 
tion Group of Japan. G.V. Keer (1955) 
reported his observations of M,-waves too. 
In all these experiments, they have drawn 
their inferences from their observations at the 
ground surface. In the present experiment, 
two surface waves were very cleary observed. 
In addition, observing their amplitude distribu- 
tions within the ground, we can reasonably 
consider these waves as corresponding to M,- 
and M,-waves theoretically obtained. 


$1. 


Instruments 


Experimental Procedures 


The instruments employed in our experl- 
ments are those used in usual seismic prospect- 


ings. 
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The seismometer is of the electromagnetic 
type (Electro-Technical Laboratories’ Seismic 
Detector, E. V.S.-2A) with a vertical com- 
poent, the natural frequency being 7 cycle/sec., 
and the damping being critical. 24 seismo- 
meters were used; 12 of them were con- 
nected with attenuators to the galvanometers 
with the natural frequency of 30 cycle/sec., 
the remaining 12 being connected with ampli- 
fiers (E.T.L. P.R.T.-type) to the galvanometers 
with the natural frequency of 100 cycle/sec.. 


Location 


The observations were carried out by the 
‘‘Seismic Exploration Group of Japan’”’ in 
September 1956, in the basin of the River 
Mogami to the east of Sakata City. Accord- 
ing to the bore-hole test to the depth of about 
60m, the layer consists of clay and sand. A 
seismic prospecting was carried out at the 
same time. The time-distance curve obtained 
is shown in Fig. 1. In a superficial layer, the 
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Fig. 1 Time distance curve by the refraction test. 
(1) Shot depth; 5m. 
(2) ” 30/10 
(3) Y : 3m, 


longitudinal wave velocity is 430 m/sec., while 
in the lower layer it is 1550 m/sec.. The 
velocity distribution down to the depth of 4m 
was complicated. This velocity distribution is 
shown in Fig. 2. As is observed from the 
refraction test, the field seems to be approxi- 


Ge) 
“ 102m /sec 


10° 12) 14 sie 18720 


Depth G 
8 

10 

(m) 42 


24 


Fig. 2 Velocity distributions with depth. 
(after K. TAZIME) 


mated by a superficial layer of thickness of 
about 100m on a semi-infinite medium. The 
velocity contrast between the two layers is 
fairly sharp. 


Shooting Procedures 


As is shown in Fig. 3, two shot holes named 
A and B were made 175m apart, their depth 
being 60m. At four points between these 
two shot-holes, 5 bore-holes were dug to the 
depth of 5m, 10m, 20m, 30m and 60m, to 
set seismometers. The holes were 2m apart 
and perpendicular to the measure line (A-B 
line). The casing pipe was set up 2~5m 
shallower than the bottom of the seismometer’s 


Table 1. 
No. Shot Depth | Shot Hol Charge 

1 lm A 225 x ; gr 
2 | 60 m A 225 gr 

3 lm B ” 

4 5m B w 

5} 10m B w 

64 20 m B la 

7 30m B Wd 

8 60 m B ” 
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Fig. 3 A map of disposition of seismometers and shot-holes in the experiment. 
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holes, to avoid the casing effect. 24 seismo- 
meters were disposed, 20 at the bottom of 
these holes and 4 at the ground surface. The 
depths of shot points are given in Table 1. 
Dynamite was used for the explosion. The 
amount of charge is given in Table 1. 


§2. Experimental Results 


The examples of records obtained in these 
experiments are shown in Figs. 4 and 5. 
What is remarkable about these records is 
that two kinds of surface-waves are found 
distinctly. Of these two waves, the first has 
the phase velocity of 500 m/sec. and the 
period of motion of 0.1~0.08 sec., while the 
second has the phase velocity of 60 m/sec. 
and period of 0.22 sec.. 

Both of them continued only for one or two 
periods. As seen in Figs. 4 and 5, for the 
shallower shot depth and obseved depth, we 


Fig. 4 Seismograms obtained in the experiment. 
(Experiments No. 3~No. 8) 
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M, 1 
Fig. 5 Seismograms obtained in the experiment. (Experiment No. 1) 
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got the larger amplitude recorded. As the 
observation of three components of motion 
was not made in these experiments, we could 
not detemine the trajectories of particle mo- 
tion. The observation of three components 
was carried out only at the ground surface. 
In the records at the surface, four wave groups 
were observed. They were called wave group 
I, Il, Ul and IV, respectively (Tazime 1957). 
Group I is a body wave and group II, III and 
IV are surface waves. However, among these 
three groups of surface waves, groups II and 
IV are predominant. Their phase velocities 
are shown in Table 2. The particle orbit has 
the same rotation in group II as that of the 


Table 2. 

Group | Phase Velocity 
II 500 m/sec. 
Ill 100~130 m/sec. 
IV 60~ 80 m/sec. 


gravity wave, and in group IV the same as 
that of the Rayleigh wave. 

The first and second surface wave in this 
section correspond to the group II and group 
IV, respectively. 

As is stated above about the particle orbit, 


Disposed Depth 


My 


Shot-Receiver Distance 
125m 


the first surface wave seems to correspond to 
M,-wave and second one to M,-wave. 


The first surface wave 


The first surface wave has the phase velo- 
city of 500 m/sec. and the period of 0.1~ 
0.08 sec.. 

At first, the maximum amplitude of this 
wave is examined in each trace in Fig. 4. 
Taking the value at the ground surface, when 
the shot depth is 1m, as an amplitude unit, 
we got the results in Table 3. The mark— 
in this table means that no first surface wave 
is found in that trace. 


Table 3. Maximum Amplidude of M.—-phase. 
Shot | | | | | 
Seismo: lm 5m 10m /}20m 30m 60m 
meter SS ON =f) ye sie | 
Om 1 |1.3 | 0.87 | 0.40 | small) — 
5m 0.70 | 0.90 | 0.55 | 0.25 | small) — 
10m 0.28 | 0.36) 0.27 | 0.13 | === 
20m 0227 | 0.38" 0°25 1 O13) eee 
30m 0. 21° / 0.270. 19 0.01 ate | fe 


60m | | | em 


Next, the amplitude distribution of this 
wave with the depth will be studied. In Figs. 
6, 7 and 8, the amplitude of the peak and the 
trough is shown against the depth of the 
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Fig. 6 The amplitude distribution of Mj-waves 
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Fig, 7 The amplitude distribution of M2-waves 
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Fig. 8 The amplitude distributiton of M:-waves 
with depth 
Shot-receiver distance; 75m 
Shot-hole depth Sapa 


observed point. As is seen in these figures, 
the amplitude is zero at the depth of about 


10m, and at the deeper points the phase is 
reversed. In other words, the vibration mode 
of this wave is such as that which has a loop 
at the ground surface and a node near the 
depth of 10m. But the node near 10m deep 
is rather unstable, as is seen in | tea o a Gal 
order to examine more minutely the amplitude 
distribution of this wave with the depth, 
choosing the case of Fig. 5, we study how it 
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Fig. 9 The amplitude distribution of Mj-waves 
with depth for every 1/100 sec. 
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changes at each 1/100 sec. and get the results 
in Fig. 9. According to Fig. 9, the place of 
the node shifts as the time goes on. It is 
near 8m deep from the begining to 4 x10~? sec., 
near 10m deep from 4x10-? to 8x10-*sec. 
and then it shifts to the depth of 12~13m. 

It is expected that when the explosion takes 
place near the depth of node (shot depth 10 m) 
the amplitude of the wave will be decreased. 
Fig. 4 shows that such is the case in our ex- 
periments. 

In short, this wave has the following character 
istics. The amplitude distribution is such as 
that which has a loop at the ground surface 
and a node near the depth of 10m. The 
absolute value of the amplitude decreases as 
the depth of observed point increases. The 
wave which shows such a type of amplitude 
distribution in the surface layer is generally 
called M,-wave. 


The second surface wave 


The 2nd surface wave has the phase velo- 
city of 60 m/sec. and the period of 0.22 sec., 
and seems to correspond to M,-wave. The 
amplitude of this wave decreases with the 
depth of observed point as seen from Fig. 5. 
In order to examine the general tendency of 
its amplitude distribution, the same method 
of analysis will be employed as in the first 
surface wave. 

Here again, the maximum amplitude at the 
ground surface for the shot depth of 1m is 
taken as a unit. The maximum amplitude of 
each trace is shown in Table 4. The mark— 
in it means that at that point the wave could 
not be observed at all. The waves clearly 
observed are limited for the shot point and 


Table 4, maximum Amplitude of M-phase. 
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the observed point, both from 0m to 10m in 
depth. 

Next, in order to examine the amplitude 
distribution minutely, the trace of Fig. 5 is 
taken as an example and its amplitude change 
at each 4/100sec. is shown in Fig. 10. This 
figure clearly shows that this wave has a loop 
at the ground surface and that its amplitude 
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Fig. 10 The amplitude distribution of Mj-waves 
with depth for every 4/100 sec. 
Shot-receiver distance ; 100m 
Shot-hole depth 2 Abin 


decreases exponentially with the depth, so that 
the amplitude is very small near 20m deep. 
But, strictly speaking, it does not decreases 
with the depth in an exponential form. 
Sometimes, the amplitude observed at the 
depth of 5m is larger than that at the ground 
surface. But in general, the amplitude may 
be considered to take an exponential form. 
This kind of amplitude distribution in the 
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surface layer corresponds to that of M,-wave. 


§ 3. Discussion 
1) Discussion upon the existence condi- 
tion 


From the experimental results, the two 
surface waves seem to correspnd to M,- and 
M;-waves, respectively. So, we shall discuss 
at first whether M,- and M,-waves can really 
exist in this field. 

R. Yamacucnt and Y. Sato (1955) have 

calculated the limit of existence of M,-and M,- 
waves. According to them, if two layers are 
both solid, M,-waves always exist. Therefore 
it is evident that in this field, the existence 
conditions for M,-waves are satisfied. 
_ On the other hand, M,-waves exist only 
when the velocity (V;, Vs;) in the surface layer 
is less or a little more than the velocity (Vez, 
Vsz) in the underlying layer. In this field, the 
longitudinal wave velocity in the underlying 
layer is larger than that in the surface layer. 
Though the transversal wave could not be 
observed, the relationship above may be con- 
sidered to exist. Therefore, this field satisfies 
the existence conditions of M,-waves, too. 


2) Comparision of the amplitude distribu- 


tions of M,-and M,-waves with those 
calculated using simple assumptions. 
K. SezAwa and K. Kanar (1935 a) made a 
study on the dispersive Rayleigh-wave propag- 
ated in a layer on a semi-infinite medium 
and got two dispersion curves corresponding 
to M,-and M,-waves. Later K. Kanai (1951) 
has calculated the amplitude distributions of 
these waves in the surface layer, using various 
assumptions. His calculated results will be 
compared with the experimental results 
obtained here. 
As this experimental field, which was con- 


M,-wave M,-wave 


Fig. 11 The amplitude distributions of M,- and 
Mp-waves with depth theoretically obtained. 
(after K. KANAI) 


sidered approximately to have a two-laryered 
structure, has in reality a more complicated 
structure, our comparsion will accordingly be 
rongh. The calculated results based upon the 
assumption most similar to the observation 
field are shown in Fig.11. As the experimental 
results are only on vertical components, we 
made calculations here only on them. From 
Fig. 11, we see that agreements between 
observed and calculated results are fairly well. 
So we may call the above two waves M,-and 
M,-waves, respectively. 


3) Wave type 


According to Sato’s (1953) classification of 
surface waves, we shall consider what kind of 
wave type M,-and M,-waves above correspond 
to. 

Generally, M, -and M,-waves have three 
types of amplitude distributions in a two 
layered structure. 

SS ES and EE EE 
EE’ EE IR EE 
meant the amplitude distribution which has 
exponential forms in the surface layer in both 
its compressional and distortional part. When 
the compressional part is exponential and 
the distortional part is sinusoidal, the distribu- 


They are By is 


tion is denoted by oo 
Table 5. 
| Cal. | Obs. 
Vp./Vpr 2.8 3.6 
Ale | 1 | Sail 
L/H | 2.23(My) 2.44(M) | 0.7(My) 2.6(Mz2) 


The relation among velocities corresponding 
to each type is as follows; 


SS 


EE V>Vei1> Vs1 (1) 
oe Vinge allied (2) 
EE V civ (3) 
EE PL 91 > 


In the case of M,-waves, 
Var2=500 m/sec., Vei=430 m/sec. 
So they satisfy the relation (1). 
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Therefore, M,-waves we observed are of a 
type. 
In the case of M,-waves, 
Vizi1=60 m/sec., Vei=430 m/sec. 
(Vs, is unkown) 
So it cannot be determined at present which 
relation between (2) and (3) they satisfy. 


Therefore both oe and EE type are possible 


EE 
in our case. 
If we have 
Vr1 ut Visi 
Ver Voz 


as usual we have Vs1>500 x 430/1550 m/sec., 
Vse > Vu. for Mz waves. 
Therefore, 
Vsi > Van 
This relation satisfies (3), so that M,-waves 


¢ LE 


we observed will be o type: 


4) Comparision with other experiments 


K. Tazime (1956) and H. Suma (1956) have 
observed M,-waves in the explosion by Seismic 
Exploration Group of Japan. G.V. KELLER 
(1955) also reported the observation of M™.- 
waves. All these studies are on the observed 
results at the ground surface. They drew 
time-distance curves from the arrival time of 
peaks and troughs in each trace and then 
classified the wave group according to the 
difference of the phase velocity. 

Seismic Exploration Group of Japan has 
classfied wave groups into group I~IV and 
G. V. KELLER into group A~C. Among these 
wave groups, I, HI, IV and C belong to the 
surface waves. They drew the trajectories of 
particle motions due to these waves and deter- 
mined M,-waves or M,-waves according to 
whether the motion was retrograde or not. 

After this, they made Fourier analyais of 
the wave form, drew the dispersion curves 
and compared these dispersion curves with 
theoretically calculated results. 

But they did not concern with the problems 
within the ground. 

The writer previously published a study 
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about the SToNELEY-wave generated by explo- 
sions. 
by the observation within the ground. This 
wave is a kind of interface wave, which has 
an especially large amplitude in the interface 
between two layers of clay and sand. K. 
Sezawa and K. Kanar (1939) formerly called 
this kind of wave also the StonelEy-wave. It 
has some difference from the definition of 
the SToneLEy-waves (R. STONELEY 1924). K. 
Kanat (1955) has recently demonstrated that 


this wave is M.-wave belonging to oe -type. 


So, the interface wave observed by the writer 


EE —type. 


: ; ’ f 
was, strictly speaking, M2-wave o EE 


The StoneLEy-wave also belongs to eas 
EE 
type. 
The wave types of the waves observed by 
K. Tazime and H. Sura are not clear. Those 
of waves observed by G. V. KELLER are 


EE “Sess 


a ee i t 
M,-wave EE EE —> EE ype 
, BSUIs SS. 
M,-wave ER EE type 


as some functions of the wave length. 


5) Dispersion 


The writer could not draw the dispersion 
curve. It is because there are a few records 
of different epicentral distance and because 
the waves in records continue only for one or 
two periods. K. Tazime (1956) reported the 
same difficulty. According to him, this difficul- 
ty is due to the fact that in the explosion 
waves the Arry phase is genrally emphasized. 

However, it is clear that the observed waves 
are the dispersive RAYLEIGH waves. It is con- 
sidered so by the fact that the vibration mode 
of the waves changes with the time, which is 
evident from the amplitude distribution with 
the depth shown above. It is known that the 
vibration mode in the ground changes as the 
wave length changes in the dispersive RAYLEIGH 
waves in general. 

If the spread length were long enough, as 
in G. V. KeLLER’s experiments, the dispersion 
curve might be drawn and we might observe 


It was based upon the results obtained — 
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other M,-waves than of = -type or other 


EE 
M,-waves than of EE ‘Pe: 


§4. Remarks 


By the experiments stated above, some 
characteristics of the dispersive Rayleigh-type 
waves have been made clear. 

This study will be important concerning 
the problem of ground noise in the reflection 
shooting method in seismic explorations. 
For the ground noise is mainly caused by the 
surface wave. 

Recently F. Press and M. Ewine (1955) have 
studied the surface wave in the crust or the 
mantle of the earth. Their purpose is to ex- 
amine the structure of the crust. Their studies 
are mainly on M,-waves, but by making the 
study of M,-waves at the same time, the 
structure of the crust will be better determined. 
But M,-waves have been found in the records 
of natural earthquakes only by K. Srzawa 
and K. Kanat (1935 b). To draw the disper- 
sion curves of M,-wave in natural earthquakes 
will be an important problem. 

For the sake of convenience, the assumption 
A= is made in almost all theoretical calcula- 
tions. But in most of our experimental fields, 
A>w. By this reason, it is inevitable that 
the comparison between the theoretical results 
and the experimental ones, is somewhat 
qualitative. 
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Minimum Group Velocity, Maximum Amplitude and Quarter 
Wave-length Law. 


LOVE-waves in Doubly Stratified Layers. 


By 


Kyozrt TazImMe 
Department of Geophysics, Faculty of Science, Hokkaido University, Sapporo. 


Summary 


The period characteristics of the amplitude function of Love-waves were obtained by 


the use of the dispersion curves. 


The quarter wave-length law may be approximated by 


the period at the AIRY phase, but it can be better approximated by the period at the phase 


corresponding to the maximum of the amplitude function. 


These approximations, more- 


over, become respectively all the better with the increase of the rigidity ratio between 


superficial and lower strata. 


Differences between the maximum of the amplitude function 


and the maximum amplitude are also considered. 


§1. When many wave groups are generated 
by an explosion, the seismic energy of each 
wave group will be progagated with its own 
group velocity. It seems to be natural that 
the maximum energy may be propagated with 
the minimum group velocity, though this 
conception has been arrived at somewhat 
implicitly. 

It is known, for example from JEFFrREY’s 
(1931) book, that the phase corresponding to 
the minimum group velocity, called the Arry 
_phase, has a special expression, including the 
Atry’s integral, for the amplitude. Amplitudes 
of the other phases, on the contrary, are ex- 
pressed by another expression having the 
derivative of the group velocity in its demoni- 
nator. 

The avove mechanism of wave propagation 
was successfully investigated by PEkErIs (1948) 
in detail and his theoretical consideration was 
confirmed by the observation by WozeE.L and 
Ewine (1948). Though the discussion by 
PEKERIS was confined to water waves, the 
present author (1956a) has observed that seis- 
mic waves in the earth are also apt to have 
large amplitudes near the Airy phase in the 
wave group. 

Therefore it can be considered, qualita- 
tively at least, with more reliability than ever 


amplitude at its Airy phase. But, quantita- 
tively, it must be ascertained whether the 
amplitude will reach the maximum mear or at 
the Airy phase. In order to make clear this 
question, LovE-waves will be taken up, as the 
subject of discussion, in this paper. 

§2. The line source of SH-waves at z=£ in 
the superficial layer resting on a semi-infinite 
layer, such as is shown in Fig. 1, can be 


Fig. 1. The depth of a source is smaller than the 
thickness of the superficial layer. 


expressed by 
do=r Hy (ki nexp(iot), 
= \" exp (jut iE Fin (B—2)y 
ke! 1 


2z2E, Goa be 
where 
Me per a OE ee ae) 
ki =o/0,, k,=0/V2, P=g+(h—z)*, 


that the dispersive wave has the maximum and the following notations are used: 
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thickness of the superficial layer, 

po; density, 

zw; «rigidity, 

v; velocity of the bodily SH-wave expres- 

sed by Vy/p, 

w; angular frequency, 

&; wave number in x-direction, 

H,®; Hanxew’s function of the second kind 

and zeroth order, 

t; V-—l1 

%; displacement of the SH-wave generated 

from the source, 
and the suffix 1 means the quantities within 
the superficial layer but the suffix 2 means 
the quantities within the lower layer. 

The initially generated SH-wave will be 
influenced by the boundary at z=0 and z=H 
in Fig. 1. Then the secondary SH-waves come 
respectively into existence, within the super- 
ficial and the lower layers. Thus the total 
displacements in the two layers become respec- 
tively as follows: 


br=2l ef(t—-Ex) COS 912 


e-' 8 + ei 8 e-2i A dé 
So 


1—Ke-2'14 "1 
2B, 443) 
=2\" e'(t-£2) cos7,E 
e-'2 + et Ke- 24 dé 
1— Ke-2!n4 > n ) 
2>E, (4) 
and 
aol” ef (wt —Ee)e-tn2(2-F) cosy, 
Sa + int Keim dé 
1—Ke-171 nn , (9) 
where 
He(1—x?) /(1+7) (6) 
1 11 


is considered to be the reflection coefficient of 
the SH-waves, and the following new nota- 
tions are employed: 


$1; total displacement in the superficial 
layer, 
ys; total displacement in the lower layer, 


%;  Yvigidity ratio expressed by s19/s1. 
Now the above integrals have poles «’s at 


£’s given by the next equation 
F(E)=1—Kexp(—277,H)=0. (7) 
Love waves, as is known, correspond to the 
residues semi-around the above poles. Con- 
sequently any order of Love-waves will be 
expressed as follows: 


[¢:]Love= 72 X (residue), 


= Ans | etcet-e ens cosy 12 


I 
nF’ ce 
= —27A(x)[e*-§=)cosyicosyiZ]e-«, (8 ) 
and 
[Paltova= —27 A(x) [etCot-E)@- tna“ #) 
x cosy cos A]e- «x, (9) 
in which 
A(«)=[—22/{9F’ (E) He-« (10) 
is called, by Saro (1952), the amplitude func- 
tion of LovE-waves. 
Introducing new notations again, the ampli- 
tude function will be expressed by 


Li lees =sing, {EOF OSeoe 


a Woke ZA+Z2) 
oH 
a (11) 
where 
Z=7secO;(sin?d,—sin29,)'/2, és 
sind,=2;,/c, sinO =2v,/v, ( ) 


and ¢ is the phase velocity of Love-waves. 

Moreover, the characteristic equation of 
Love-waves, that is (7), can be expressed with 
these notations 


oH]/v,=secd, tan“!Z. (13) 
But one has from (12) 
dZ/d0,=tan,7?cos*Osec?6,/Z; (14) 


the next relation will be derived from (11) 
and (13) 
HH do _ sec, 
v, dO, A(k)’ 
Considering the group velocity U, on the 
hand, one has 


ie aie Sa Se do } 
= fon, (sin 9,+-wcos8; ral (16) 


(15) 


using the relation 


J 


} 


Minimum Group Velocity, Maximum Amplitude and Quarter Wave-length Law. 45 


kK=(o/v,)sing,. 


At last, eliminating 0, and d0,/dw from (12), 
(15) and (16), one has 


2a A(c)={T/(H/v1) }{1—(oi/c)?}-{(v1/U) 


—(v,/c)}, (17) 
= (A4/H){(c/v1) —(v/c)}-{(v1/U) 
—(w/0)}, (18) 


where 
T=2r/w and A=cT 


mean respectively period and wave-lengths in 
x-direction. 

Here c/v; and U/v, vs. T/(H/v;) or 4/H are 
considered to be given in the dispersion curves, 
or can be rather easily obtained by means of 
the characteristic equation (13). When the 
dispersion curves are given, the amplitude 
function of Love-waves will be obtained at 
once by (17) or (18). These relations may be 
satisfied by any order of Love-waves, so the 
amplitude function of any order of them can 
be calculated in the same manner if the dis- 
persion curves of that order are known. 

Two examples of the amplitude function 
calculated by (17) are illustrated in Figs. 2 (a) 
and.2 (b) when p.=p; but ~=4 and 30 respec- 
tively. The dispersion curves corresponding to 
these were calculated by means of the nomo- 
graphs drawn by Saro (1956). 


— T/H(H/M) 
(a). y=4. 


70 
3,0 


3.0 


20 


S 7 I yrs A 


2 3 
— T/(H/vi) 


(D) 30; 
Fig. 2 The amplitude functions of the zeroth 
order of LOVE-waves and the dispersion curves. 


§3 As the peresent author (1956b) has al- 
ready pointed out, the period T./(H/v:) corres- 
ponding to the Arry phase does not coincide 
with the period T,/(H/v,)=4 which represents 
the quarter wave-length law, though they will 
approach each other when the rigidity ratio 
becomes large. In Fig. 2, it is seen moreover, 
that the period T;/(H/v1) corresponding to the 
maximum of the amplitude founction does not 
coincide with T./(H/v,) as well as Ty/(H/v:). 
But it must be noticed that TJ, is nearer to 
T; than T., especially at small rigidity ratio 
when the quarter wave-length law is hardly 
satisfied by Tu. 

T./(H/v;) and T;/(H/v1) for more wide range 
of rigidity ratio are illustrated in Fig. 3 by 
reference to which one can more easily under- 
stand the explanation above given. 


0 (es ay See aT ee ern 
\ Sts) 19S si 7161 10, 20 30 50 70 100 
————— f/ 


Fig. 3. -The periods corresponding to the minimum 
group velocity and to the maximum of the 
amplitude function. 


46 Kyozi TAZIME 


But, returning to eqs. (11) and (16), the 
respective stationary phases of A(«) and U(«) 
will be considered, apart from the graphical 
representation. A(«) and U(«) are respectively 
the functions of @; alone through (12) and (13). 
From (11) one has 


aA GQ) Nae dw) au 
inj a LN Te) EGE 
fae do) sec {tanti( in) ae ot 


Therefore #, and w, corresponding to 


dA/d6,=0, (20) 
will be given by the relation 
@w/d0,2—tan0,(dw/d0,) =0, (21) 
except when 
sec0,(d0,/dw)=0. (22) 


On the other hand, one has from (16) 
—2 20 2 
ae ee cost | oe a el 


dé, © \db, db w\d0, 
stan |. 
Therefore 0; and w, corresponding to 
dU/d0,=0, (24) 
will be given by the relation 
@a 2/(dwo dw 
ee 6, = 26 
dO ee) i ee 
except when 
Uoacos6\(d0,/dw)=0 ° (26) 


As we have no interest, in the present dis- 
cussion, in the cases when c=v, and v2, cases 
(22) and (23) will be excluded. Then @, and 
© common to (20) and (24), if they exist, must 
satisfy the following relation 


(1/w)(do/d6,)=tan 0, , (27) 
from which it follows that 
w=constant x sec 0, , (28) 
which is rewritten by the use of (13) 
tan="Z= constant. (29) 


Differentiating this with respect to 6,, one has 
at last 
tan 9,7? cos?@ sec?/{Z(1+Z?)}=0. (30) 
But there is no 6, satisfying this equation 
except #,=0. Yet it is known as to Love- 
waves that 


r/2>0,>0 . 


(31) 


Thus eq. (30) can be satisfied only when 
(GEAW (ir 9= 82 5 
The analytical process of discussion just 
given has arrived at the same result as that 
which was already understood by the graphs 
MUN IS Toe 
$4. If x is large but #, remains different, 
the larger the better, from 7/2 as well as 
from O, Z in (12) will also be large. In this 
case, (13) will be reduced to 
(wH]/v;) cos 0,= (2n+1)z/2 , 
or 
H={(2n+1)/4}(w,/cos 6,) T. 
The quarter wave-length law in general can 
be expressed by 
H={(2n4+1)/4}0iTy . 
Therefore it follows 
Pho— 7 /cos re 
in which cos 6, for T, and TJ; are estimated 
as in Table 1. As cos@, is larger than 0.90 
for x>10, (32) may be approximated fairly 
well by 7: in place of Ty. 


(32) 


Table 1. Cosé@, at the phases corresponding to 
the stationary values of group velocity 
and the amplitude function. 


ll Gee a es 
PG" L352 hy a7% oleae ser Oaes 
Pee ee eee ee 

Pe | cosa 0.47 0.79 0.88 0.96 


0.68 


T: | COS 41 OD 10.88 0292 0.94 0.96 


As explained already, the present discussions 
need not be confined within the zeroth order 
of Love-wave alone, but is common, at least 
qualitatively, to all orders of them. Several 
values of T;/T, are shown in Table 2 where 


Z in (12) is not restricted to a large value 
alone. 


Table 2. The values of T/T, for various orders 
of LOVE-waves. 
xX | 1) In Ly Ly Tn 
4 | 0.90 0.86 0.86 0.86 0.86 
8 | OFOZ ORL 0.91 0.91 0.91 


ieee 
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If the phase corresponding to the maximum 
of the amplitude function is observed, the 
calculated error of A, i.e. the thickness of 
the superficial layer, will be smaller than 5 
~10 %, owing to the neglection of the diffe- 
rence between JT; and Ty, even for not very 
large values of x. 

Investigating the dispersion curves, it may 
be recognized that v,/c is much smaller than 
unity and v,/U—v,/c is almost constant for 
any order of Love-waves if the stationary 
phase of the amplitude function or the group 


——> c037z 


—— T/ (H/vi) 


—» exp{-it(z-4)} 


Omblige 2a Seo 40.5 
—> T/(H/%) 


velocity is concerned. Thus roughly the 
maximum of the amplitude function can be 
taken as proportional to 7; in (17). Very 
rough amplitude ratios at stationary phases for 
various order of Lovr-waves are estimated in 
the lowest line of Table 3 in which two examples 
calculated by Sato (1952) are also tabulated. 

Until now, the depths of the source and the 
observing point have been neglected. How- 
ever, to tell the truth, the effect of the depth 
is again a function of the period, because, 
from (2) and (11), 


N 
iS) 
a 
° 
° 


(a). 


fi/I=A 
j1//A=30 


10 


(Z-HY/H 


(b). 
Fig. 4 The effect due to the dépths of the observing points. 
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= (w/v) Cos A, eo —i(w]v,)V sin? 6,—sin? 9 e 


es 28 COS Ui ee to TT ee {_2AV sina, sine : Z-dtt 
cos m2=0054 T|(H]01) at a ee to T|(H]01) H 


(33) 


ble 3. Amplitude ratios at the maximum of : mn ior be [ ; 
Ae fate eye =f Pieeerains corresponding to the maximum amplitude may 


a teens AOL ber cistituedne ha aD van ne effect of the 


x ie In Ly T3 ly depths mentioned above, though it will become 
r y 7 1) Tene somewhat large when the depths approach or 
12 3.6 6.0 8. : j 
. : : hickne f the superficial 
8 7 ee ey ees | i | ee is fhe) CUE heeanes B 
a 1 1s 6 Wy Mg NF 


§5. If a line source lies within the lower 
layer, as shown in Fig. 5, 


Here @, and T/(H/v,) are related by the 
characteristic equation or the dispersion curve. 
Four calculated examples as to cos7,z and 
exp{—7.(z—A)} are shown in Fig. 4. Repla- 
cing z by £, cos7,# and exp{— 17,(E—H)} 
can be also estimated by the formula (33). 
Now comparing Fig. 2 with Fig. 4, one sees, 
the gradient of the corresponding curve is 
more sharp in the former than in the latter, 
especially near the period corresponding to 
the Amy phase and the maximum of the Fig. 5. The depth of a source is larger than the 
amplitude function. Therefore the period thickness of the superficial layer. 


expression (1) must be replaced by 
by=n2Hy©(kzr) exp (twt) 
-| exp (inte Fin(E—a) , 2z2E. (34) 


Consequently total displacements in the two layers will be respectively as follows: 


=o" gitwt——ne 2G" cos Ue sae ‘ dé 


ae 1— Ke-2#112 ent x,’ 


=|. pttategey, 2. 8 COS Wiz d. 
22 "ye “Oe (ctr nee Elan, (35) 
and 


e-'ng(k-A) 


in| ei(wb-Ex) p 4 ae {eln2@-4)(] — Ke-214) + e-ina@-2)(g-2in, # — K)} dé ' 2<E : (36) 
ne ai Ne 


oe i(wl-fm) ena ing(K-A | i dé 
= e ° [Ligoeiae {ein(e- (1 — Ke-212) + e-t2(8-4) (e211 HK )} VZle G7 
<0 th 
where 
482% , rig 
% m 1-K~ 


Love-waves, as before, correspond to the residues semi-around the poles, given by (7) 
of the above integrals, : 
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[1 Love = —27 A(R) z[e!@'-e-2F-) cog mA cos mZ\¢-x , 
[Peliove = —27 A(R) z[e! tf) e- t(D e-tn2@-#) cos?y, A] en « 


It is seen that (8) and (39) will remain the 
same even if z and EF are exchanged. Thus 
the complete reciprocal theorem is concluded 
when the source and the observing point lie 
within the same layer. In order to conclude 
the reciprocal theorem between different layers, 
on the other hand, one must take, by com- 
parison of (9) with (38), “4.4)@(kr) for 
Ay(Rir) in (1) and 4,A)©(kr) for Hy (kr) 
in (34). 

Now it has been found that amplitudes of 
Love-waves in this case are expressed expli- 
citly also by the same amplitude function as 
that obtained before. Therefore the discussions 
as to the present title will yield nearly the 
same conclusions with each other within 
whatever layer a source may lie. 

But it must be remarked that the quarter 

wave-length law is a nature of the stationary 
phase of A(«) which appears directly in the 
case of no depth of the source and the observ- 
ing point. 
§6. StToNnELEY (1937) estimated the depth of 
20° discontinuity layer in the globe by means 
of a similar formula to (32). Srzawa and 
Kanar (1937) discussed his estimation. How- 
ever STONELEY used the formula merely impli- 
citly having no strict physical meaning. 
SezAwa and Kanal have reached already the 
theoretical conclusion similar to that described 
in this paper as to the estimation of the 
thickness of the superficial layer. But they 
were too hurried in numerical calculations and 
too nervous about practical quantities in natu- 
ral earthquakes to get the physical meaning 
of the formula or the quarter wave-length 
law. In this paper the behaviour of the 
normal mode propagation of seismic waves in 
general has been kept in mind, although Love- 
waves alone are treated. 

The sharpness of the maximum amplitude 
will grow, as shown in Fig. 6, with increase 
of the rigidity ratio which may reach con- 
siderable values near the surface of the earth, 
as in the case of seismic prospecting. Thus 
a layer may play the part of a considerably 
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(38) 
(39) 


narrow band pass filter for Love-waves. The 
quarter wave-length law, as the author (1956b) 


H/4=30 


—> T/ (W/) 
eel “a eee noes eee 
9 2 4 6 
Fig. 6. Amplitude functions for various rigidity 
ratio. Lo, 02=(1. 


has described, can be better approximated at 
the Arry phase of RayieicH and SEzAwa- 
waves than at that of Lovre-waves. The 
present author believes that a layer will serve 
as a narrower band pass filter for RAYLEIGH 
and Sezawa-waves than for Love-waves, though 
full investigation is left for the future. 
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Relations between Charge Amounts and Periods 


in Resulting Seismic Wave Groups. 


By 


Kyozi TAzImMe 


Department of Geophysics, Faculty of Science, 
Hokkaido University, Sapporo. 


Summary 


Variations of periods in superficial seismic waves were observed with ‘various charge 


amounts. 


Those seismic waves were considered to be the zeroth and the first orders of 


RAYLEIGH and SEZAWA-waves respectively, due to the superficial layer in the experimental 
ground. Each period was found to have its own asymptotic value which was seen to 


satisfy the quarter wave-length law discussed previously by the present author. 


Frequency 


characteristics at the origin of seismic waves were also considered. 


§1. Relations between amplitudes in seismic 
waves and charge amounts have already been 
investigated by many persons, such as HapBeEr- 
JAM and WHETTON (1952), RicHarD and PieEv- 
cHoT (1956), GaskELL (1956), Suzuki and Sima 
(1956) and so on. They have respectively 
offered experimental formulas for the relation- 
ship between amplitudes and charge amounts, 
not always coinciding with each other. How- 
ever, few investigations have been made for 
relations between periods and charge amounts. 

Matuzawa (1939) found that larger periods 
were contained in seismic waves generated by 
large earthquakes than in seismic waves by 
small earthquakes. Kasanara (1956) has come 
to the conclusion that the maximum of the 
period spectrum near initial motions in seismic 
records tends to be a large period with the 
increased scales of the origins, from a very 
small explosion to a very large explosion and 
to a natural earthquake. 

But their scales of comparing the origins 
are too rough to enable one to find a quantita- 
tive rule between period spectra and charge 
amounts, especially. in the case of comparing 
the seismic records obtained by small explo- 
sions. It has generally been considered, that 
few changes in periods may be observed by 
such a small explosion as that used in field 
model seismology. Some students consider 
that no change in periods will be observed, 


unless charge amounts should be varied by 
10 times at least. 

The present author (1955) has pointed out 
that an observed period of the initial motion 
is apt to be greatly influenced by the nature of 
the recorder due to its inertia. Consequently, 
few changes in periods may often be observed 
at the initial motion on seismic records, though 
periods of the earth movement themselves were 
varied. 

In order to observe systematic relations 
between periods in seismic waves and charge 
amounts, the present field experiment has 
been executed employing small explosions. 
Seismic records obtained by such small scale 
observations as this will be greatly influenced 
by very small subsurface structures. There- 
fore superficial waves, especially the higher 
orders among them, are apt to predominate in 
the direct rear of the initial motion, reminding 
the result of the previous reports by the 
author (1956a, 1956b). It will be difficult, for 
this reason, to analyze the wave form of the 
P-wave itself, even if several peaks and 
troughs near the first break only are concerned. 
Yet, as mentioned above, the first break must 
be severely influenced by inertia proper to 
the recorder. On account of these items of 
information, changes in periods of all waves 
predominant on the seismic record will be 
studied in this paper, without discussion on 
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the forms of P-wave. 

§2. Many seismic records were obtained by 
the members of the Seismic Exploration Group, 
in the summer of 1956, at a suburb of Sakata 
city in Yamagata prefecture. 

The E. T. L.-geophones for a: vertical com- 
ponent were used; their natural frequency 
was 7c.p.s. They were set in a straight line 
on the ground every 0.50 m from an epicenter. 
Total length of the experimental span was 
26.5m. Two types of recorders were used: 
one was an E. T. L-oscillograph with its own 
amplifier and the other was a Sanei-oscillo- 
graph with no amplifier. The former will 
be called A which is the initial of Akita 
University, while the latter will be called K 
which is the initial of Kyoto University. K has 
an attenuator to allow various attenuations, 
keeping the damping factor constant. The 
natural frequency of the galvanometer A is 
125c. p.s. and that of the galvanometer K is 
30c.p.s. The total frequency characteristics 
of those recorders connected with the geophone 
are jllustrated in Fig. 1. The broken line 
corresponds to recorder A and the full line to 
recorder K. 

The geophones connected to A and K are 
set one after the other, as shown in Fig. 2, 
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Fig. 2. The arrangement of the geophones. 


every 0.50m. The depths of shot points were 
1.0m in one case and 2.0m in another. A 
series of charges of various sizes were blasted 
at each depth. Starting from a cap alone, 
3x2" grs dynamite were respectively added to 
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Fig. 3. Recorder A, shot depth 1.0m, epicentral 
distances 0.5 to 11.5m, charge amounts cap 1, cap 
2, 1/720, 1/360, 1/180, and 1/90 x 2.25 kgs. 


a cup when 7=0, 1, 2,3,4,5 and 6. In order 
to get constant amplitudes on the seismic 
records, the gains of the recorders were ad- 
justed properly in accordance with the diffe- 
rent charge amounts and different epicentral 
distances. Considerable technical difficulties 
were experienced in making the gain adjust- 
ments. 

The seismic records thus obtained are 
shown in Figs. 3 to 10. Time marks in the 
records indicate 1/100sec and a shot mark is 
seen in the lowest trace in each record, 
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‘Fig. 4. Recorder A, shot depth 1.0m, epicentral 


distances 15.5 to 26.5m, charge amounts cap 1, 
1/360, 1/180, 1/90, 4/90 and 8/90 x 2.25 kgs. 


§3. The time-distance plots of every trough in 
the seismic records are shown in Figs. 11 and 
Table I 
wave phase classifi- 
group velocity cation 
I 4.3 x 10? m/s direct P 

15.5 refracted P 
I Bppay ie) PAW) Mo. or My. 
Ill 120) toi 422 Mo; 
IV 0.6 to 0.8 Mu 


12. Keeping the classification of wave groups 
in mind, predominant wave groups in the 
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Fig. 5. Recorder A, shot depth 2.0m, epicentral 
distances 0.5 to 11.5m, charge amounts 1/720, 
1/360, 1/180, 1/90 and 2/90 x 2.25 kgs. 


present seismic records are classified by dis- 
tinct phase velocities, as tabulated in Table I. 
The I wave group in Table I is a direct or 
refracted P-wave. The other three groups are 
all considered, on the basis of previous ex- 
periences of the author (1956a, 1956b), to be 
RAYLEIGH and SEzaAwa-waves, i.e., M-waves. 
The great branches and the orders of M-waves 
appears respectively to be as shown by the 
classification on the right hand column in 
that table. 

Indeed the phase velocities of the above 
superficial waves have tendencies to increase 
with increase of charge amounts, but these 
tendencies are not so significant that phase 
velocities can be considered useful for classi- 
fication of wave groups in time-distance plots. 
As to the relation to charge amounts, the 
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distances 12.5 to 23.5m, charge amounts 1/360, 
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periods of the waves will be taken up at 
first. 

§4. The observed relations between the 
periods and charge amounts are obtained by 
the aid of the above classifications of the 
wave groups, as shown in Fig. 13 where the 
observed periods by recorder A as well as by 
recorder K are plotted directly at the same 
time. In spite of the different frequency 
characteristics between A and K, few diffe- 
rences were found between the periods observed 
by A and by K. It seems to be certain that the 
period may increase a little in accordance with 
increasing distance from the epicenter, but this 
effects is too small to be taken up in the 
present problem. 

Fluctuations of the plots in Fig. 13, on the 
other hand, are mainly due to fluctuation of 
the effect of an explosion. In the present 
experiment, amplitudes were often observed 
quite different from each other, in spite of the 
use of the same charge amounts. Technical 
troubles, as noted above, were experienced very 
much in making the gain adjustment,  Fluc- 
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Fig. 7. Recorder K, shot depth 1.0m, epicentral 
distances 1 to 12m, charge amounts cap 2, 1/720, 
1/360, 1/180 and 1/90 x 2.25 kgs. 


tuations of the effect of an explosion must 
cause fluctuations of amplitudes as well as of 
periods. 

In the present seismic records, separation 
of the III wave group from the IV wave group 
is often ambiguous, especially near the epi- 
center, and the former is apt to be covered 
by the latter which is very strong everywhere. 
The periods of the III wave group were not 
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Fig. 8. Recorder K, shot depth 1.0m, epicentral 
distances 16 to 27m, charge amounts cap 1, 1/360, 
1/90, 4/90 and 8/90 x 2.25 kgs. 


plotted in Fig. 13, though they seemed to be 
almost equal to those of the IV wave group. 
Now it must be evident in Fig. 13 that the 
period of each wave group excepting a P-wave 
is increased by the increase of charge amounts. 
Moreover it must be noted that these increase 
in period become gradually slow as the charges 
‘grow larger, and each period will reach its 
own asymptotic value. Further it must be 
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Fig. 9. Recorder K, shot depth 2.0m, epicentral 
distances 1 to 12m, charge amounts 1/720, 1/360, 
1/180, 1/90 and 2/90 x 2.25 kgs. 


remarked that the period of the II wave group 
will reach its asymptotic value sooner than 
will the IV wave group. Further it will be 


found in every wave group that 
((ipacymp =. ( S10" SECs) 
and ( Tyv)asymp == VA 10-2 SEC, 


therefore 
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Fig. 10. Recorder K, shot depth 2.0m, epitentral 
distances 13 to 24m, charge amounts 1/360, 1/90, 
2/90 and 4/90 x 2.25 kgs. 


(Tyv)asymp/(T1pasymp =3 . ( 2 ) 


On.the other hand, the author (1956 a, 1957) 
suggested previously that the period corres- 
ponding to the maximum amplitude was almost 
equal to the period T, at the Airy phase which 
might be approximated by the following for- 


mula: 
4 (3) 


H=(3/4)(UmTa) , 
or AH=(j4)\@n fa) ’ 
where A and vp, mean the thickness of the 
superficial layer and the velocity of the P- 
wave in that layer. 

By the comparison of (1) with (3), it will 
be found that (7ppssymp and (Tyy)asymp in (1) 


for Mi. and Moy 
for Mu and Ms, 


correspond respectively to J, in the upper 
and the lower lines of (3). So one has ar- 
rived at the conclusion that the stationary 
phases of the wave groups will be observed 
as a result of a large explosion. On the 
contrary, any phases having smaller periods 
than those which correspond to the stationary 
phase will be observed as a result of a small 
explosion. 


$5. In order to get more concrete comparison 
of (1) with (3), the subsurface structure of 
the experimental ground must be known. The 
time-distance plots of the first break are ob- 
tained by the refraction shootings, as shown 
in Fig. 14. The part indicated by a thick line 
in this figure corresponds to the spread in the 
present experiment, The subsurface structure 
was analyzed by usual refraction method as 
follows: 

H=19im? 

Vy = 4.3 X10? m/s , (4) 


Vp2 = ele m/s ’ 


the boundary planes being almost horizontal. 
From (3) and (4), one has 


Lea=18 *10-2'séc’, 


coinciding approximately with (Tyy)asymp in 
(1). In discussions on numerical values in a 
field experiment, errors of 10% must be allow- 
ed in general. Thus the above coincidence 
must be considered rather good, though 
Tasymp Might be increased to a larger value 
than J, as a result of the increase of charge 
amounts. 


§6. The wave form R(¢) on seismic records 
and its frequency characteristics Y(w) are 
related by the following relations; 


R(t)= ve Y(w)do , 


Y= 75. | e- RL) dE , 


and 
by means of Fourier transformations. 


The same relations would be satisfied at 
the origin of seismic waves, 
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taken respectively as Y2(w) and Y3(m), the 
following relations will be obtained, 


St)= age) ee Vileddo , 


oli co 
and Wi@)= =| ett f (t)dt ‘ RA = Y.(o)- Yalo)s7 (@) ’ 
na Y(w)=YV,(w): Yx(o)- Ys) . 
If the frequency characteristics, i.e., the 
transfer functions of the propagation mecha- . In this experiment, charge amounts alone 
nism through strata and of the recorder, are were variable, the conditions of the strata 
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relations between the periods and charge amounts. 
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Fig. 15. The period characteristics considered at 
the origin, 
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The time-distance plots of the first break. 


and of the recorder being kept unchanged. 


Therefore the observed variation on wave. 
forms must be attributed solely to the variation 
of wave forms or frequency characteristics at. 
the origin, because not only Y3(w) is proper 
to the recorder, but also Y.(w) is considered. 
to be proper to the propagation mechanism 


through strata. If Yi(w) was also unchanged. 


by the variation of charge amounts, the wave: 
forms observed at some fixed distance from 
the epicenter should be always equal to each 
other in the different cases; the larger gain | 
of the recorder was used with smaller charge 
amounts and the smaller gain was used with 
larger charge amounts. However the observed 
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fact was contrary to this, that is the recorded 
periods of the wave groups were found to be 
increased with the increase of charge amounts. 

Because Y;(w) illustrated in Fig. 1 has no 
peak near (Tyy)asymp aS well as near (Tivasymp; 
Y, (@) must have considerably sharp peak 
Beate t Pii)asgmp. Or. (Lry)acsymp.. On the other 
hand, the breadth of the band of Y\(@) would 
be broadened towards large period with the 
increase of charge amounts, as considered in 
Rice 15, 

If Y.() is known, as in the case of Loveg- 
waves, Y,(w) or f(¢) may be acquired, because 
R(t) and Y;3(w) are considered to be given. 
Yet in this experiment it was not possible to 
observe Lovse-waves. Unfortunately Y.(w) of 
M-waves in general has not yet been theore- 
tically obtained. 

Finally, a note must be added that observed 
amplitudes are due to two factors: one is 
absolute height of Y,(w) in Fig. 15 and the 
other is |Y,(w)-Y.(w)|. Thus a relation bet- 
ween amplitudes and charge amounts should 
be also dependent upon subsurface structures. 
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-Summary 


temperatures up to its melting point. 


\ 


For the purpose of investigating elastic properties of rocks at the melting point, the 
velocity of ultrasonic dilatational wave in polycrystalline sodium»was measured from room 
Very marked drops in elastic moduli with the 
corresponding decrease in the velocity were observed as the melting point is approached. 
The marked drop in the velocity near the melting point was compared with the order- 
disorder transition phenomena of a metal alloy. 


The nature of the core-mantle boundary 


(Bullen’s region D’’) was discussed from such experimental evidences. 


§1. Many of the problems concerning the 
earth’s interior which seismologists are now 
confronting with require for their solutions 
experimental informations about the mechanical 
and thermodynamical behaviours of rock- 
forming materials under high pressures and 
temperatures. Elasticity of rocks at high 
temperatures up to the melting point is one 
of the essential subjects of study in geophysics, 
especially with reference to the nature of the 
earth’s core. 

As the specimen material for this. sort of 
study, a chemical element, in stead of rocks, 
of which physical properties are accurately 
known is desirable for the simpleness of the 
related theoretical treatment. Thus, poly- 
crystalline sodium was adopted as the specimen 
material and a test piece of this material was 
prepared, being cut from the interior of a 
large coherent block of sodium, The ex- 
perimental procedures were described and the 
results discussed from thermodynamical point 
of view in my previous paper. (D. SHImozuRU 
1956). / ; 

, The measured velocities of the compressional 
pulse. wave in sodium are plotted against 
temperatures in Fig. 1. The velocity decreases 
almost linearly with temperature up to the 
point several degrees below the melting point 
(=97.6°C). The-velocity begins to drop very 
sharply as the melting point is further ap- 


proached, and at the melting point the velocity 
decreases discontinuously from 2735 m/sec to 
2542 m/sec. This sort of discontinuous jump 
in the velocity at the melting point was:also 
observed in the author’s previous experiment 
(D. Suimozuru 1954) in the case of Woop’s 
alloy. In both experiments, the discontinuous 
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Fig. 1. Velocity of dilatational wave in sodium 
as a function of temperature. 


jump in the measured velocity occurs in 
almost the same manner. It is very interest- 
ing that the large drop in the velocity begins 
from several degrees below the melting point 
and this phenomenon will be the main subject 
with which the present paper deals. 

Thermal dependence of the adiabatic bulk 


modulus was deduced using the following 
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thermodynamical relations and empirical values 
of involved constants: 


f*oR\* oe 

(Se) = ae OT Js 
1 /OBr 

mes raat OT ) 


De ake (2) 
BOY E 


1 (O87 
——) =4.4, 2) 
ay oT ) 
ae = —159.3x 10-7H+3 x 10-!%, 
0 
+1.6 x 10-93 , (3) 
Ta Tae. (MSHA Sey (esl Ss 
Kr 
where 


7-+++GRUNEISEN’S parameter, 
8----compressibility (=reciprocal of bulk 


modulus), 
a----coefficient of linear thermal ex- 
pansion, 
Ky-++-adiabatic bulk modulus, 
Kr----isothermal bulk modulus, 


4V/Vo----volume change due to hydrostatic 
pressure P. 


Thermal variation of the adiabatic bulk modu- 
lus is found to be expressed by: 

OK; 
(or 


Calculated bulk modulus is plotted against 
temperature in a solid line in Fig. 2. In my 


) = 1.44 x 10!° dynes/cm?, deg. (4) 
p 
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523% 10° dynes/cm? 


(00 120°C 


Fig. 2. Change of bulk modulus of sodium with 
temperature, 


previous paper, the above value of (0K;/0T)p 
was taken to be good up to the melting point) 
From the measured variation of (K+4#) and. 
the variation of K which was calculated by 
the above procedure, we get the thermal 
variation of rigidity modulus # up to the 
melting point as shown in Fig. 3. . 


2-0x 10"° dynes/cm* 


Fig. 3. Change of rigidity modulus of sodium 
with temperature. 


Rigidity modulus which decreases almost 
linearly with temperature begins to drop very 
sharply from the point several degrees below 
the melting point. Table I shows the values 


Table I 
Materials —0u/0T x 10-8 
Na 1.0 + 
Pb Hel: 
Zn tel 
Al 4.7 
Ag 2.4 
Fe De 
Cu hey 


of (0@“/0T) for sodium in the range of linear 
decrease, as well as those for other materials, 
obtained directly by T. Krxuta (1921). 

It must be noticed that the apparent rigidity 
modulus at the melting point does not vanish 
to zero but seems to remain to be finite. 

§2. Abnomal rise in the specific heat and in 
the coefficient of thermal expansion within a 
few degrees below the melting point have 
been obseryed by several authors for a crystal- 
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line solids. This abnormal rise in the specific 
heat is not the influence of impurities involved 
in the examined substance, but has been re- 
cognized as the phenomena of “ pre-melting ” 
by Usseionpe (1938). For all alkali metals so 
far investigated, there has been found a marked 
increase in the specific heat as the melting 
point is approached. Such abnormal rise in 
the specific heat is shown in Fig 4 for the 


cal K" mat 
Specific neat 


Fig. 4. Specific heat of sodium versus temper- 
ature. 


case of sodium (T.M. DaupHinesg, et. al. 1954), 
As for the elastic modulus, very marked drop 
in apparent modulus near the melting point 
would be the consequence of such pre-melting. 
The pre-melting phenomena are believed to 
be the consequence of the liquid phase which 
is formed amidst the metastable solid phase 
in the form of “nuclei” or “embryo” as the 
melting point is approached. 

For describing quantitatively the nature of 
the liquid state and of the mechanism of 
fusion, a certain parameter was introduced 
by some authors (FRENKEL 1946, LENNARD- 
Jones and DervonsHireE 1939) which specifies 
the degree of order in the arrangement of 
atoms in the substance. The perturbation of 
order in the liquid state is considered to be 
attributable to the. formation of holes, i.e. 
vacant sites of the crystal lattices. Since 
melting of solid is a co-operative phenomenon, 
the theory of order-disorder transition in 
metal alloy may possibly be applied also to it. 
The transition from the ordered to disordered 
state can usually be seen in many kinds of 


alloy for instance directly by X-ray diffraction 
methods. S. Stecen (1940) investigated the 
effect of such order-disorder transition on the 
elastic constants-of single crystal of Cu;Au. 
Variations of the three principal elastic con- 
stants Cy, C2, Cay with temperature are shown 
in Fig. 5. The three elastic constants decrease 
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Fig. 5. Variation of the principal elastic con- 
stants of CusAu with temperature. 


in magnitude linearly with temperature from 
0° to 250°C, and beyond 250°C, begin to do 
so very rapidly with a culmination at the 
transition temperature. 

Fig. 6 shows the temperature variation of 
bulk modulus of Cu;Au as calculated from 
the above Siecen’s results. Bulk modulus 
also decreases notably as the transition temper- 
ature is approached. 

The theory of order-disorder transition has 
been developed in various degrees of approxi- 
mation. But, for the qualitative evaluation 
of the problem, it will be sufficient to apply 
the Bracc-Winuiams theory which deals with 
the long-range order; i.e. to assume that all 
atoms interact with one another. The free 
energy expression according to Bracc- WILLIAMS 
theory is, after some simplifications, 


Fo i NW—N( oe eT) 
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where N----number of atoms, 


W.--+-energy per atoms. 


€ is a parameter specifing the degree of inter- 
change of positions between sodium atoms 
and holes, and is called “degree of order”. 


Bulk Modulus 
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Fig. 6. Variation of bulk modulus of CusAu with 
temperature. 


€=1 corresponds to the state of order and 
€=0 to the state of complete disorder, i.e. the 
liquid state. & decreases rapidly with temper- 
ature in the neighbourhood of the transition 
point as can be seen by the following equation, 
WE 
Py pes = 
From thermodynamical relations, the bulk 
modulus AK can be expressed as follows, 


un Op \ __ 02h O2F \ ( Or 
K=—VWau)= ni Sahel 
(7) 
The above equation shows that the bulk 
modulus AK decreases with temperature con- 
spicuously as the transition point is approached. 
This figure corresponds to the experimental 
evidences shown in Fig. 6. 
Accordingly, in our case, thermal dependence 
of the adiabatic bulk modulus, deduced in the 


€=tanh 


preceding chapter, is no longer valid as the: 
melting point is approached owing to the 
formation of vacant sites of the crystal lat- 
tices. The ‘bulk modulus should decrease 
sharply with temperature as shown in the 
dotted line in Fig. 2. Hence, the observed 
drop in the velocity of elastic wave near the 
melting point is partly due to the decrease off 
bulk modulus and partly to the decrease of 
rigidity modulus. 

According to Bripeman (1949), the change oft 
compressibility of sodium on melting is 


AB= B.=8s=34 x 10-18 cm?/dyne , 
where #8; and 8s are compressibilities of liquid 


and solid sodium, respectively. As Kia 

s! { 
=6.35 x 10!° dynes/cm? at the melting point om 
neglecting the pre-transition effect, we gett 


Kis = =5,.23 x10! dynes/cm? 


Bi 
sodium at the melting point. “ 
§3. On seeing the conclusion in the preced- 
ing section we are tempted to calculate rigidity, 
modulus of liquid which is of great interest 
not only from physical but also from geo-: 
physical point of view. The following calcu-: 
lation is based on the experimental evidences: 
and consideration on free energy at the melting: 
point. 

Free energy of solid (Fs) 
phase is defined as, 


F;,1=—kT log Z:, (T), 


for the liquid: 


and liquid (F,) 


(8) 
where Z;,7 is the partition function expressed 
by 


Ze (N=) em sr ikr . 


hy 
kT ) 
/ BR ceil : 


From the stand-point of formal thermodyna- 
mics, a transition point in a substance is 
given by the intersection point of the free 
energy curves of the two states concerned. 
Therefore, equating the two free energies at 
the melting point, we have 


Ur—U,=3NRT log (vs/vr) . 


Hence, 


Fs,1=Ur,1—3NKT log ( (9) 


(10) 
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Since the left-hand side of the above equation 
is nothing but the latent heat of fusion Z, we 
get. 

(vs/v,)=exp (L/3NRT,,)=exp (L/3RT») 


where L----latent heat of fusion, 
R----gas constant, 
Tn: -+-melting temperature, 
vs,vz+++-frequencies of the lattice vi- 
bration for solid and liquid, 
respectively. 
Putting R=1.9869cal/deg., Z=27cal., T 
=371°K in the above relation, we obtain 


Ye=1.0123y,. 


This implies that liquid sodium has its fre- 
quency of lattice vibration somewhat smaller 
than that of solid. 

In general, the frequency of lattice vibration 


(11) 


is connected with the velocities of elastic 
waves by the following expression; 


9No 1 2 
Mare) meri. hi 
where N--+-++Loscumipt number, 
p-+--density, 
M:---atomic weight, 
Up, Vs+++Velocities of longitudinal and 


transverse elastic waves. 


If vz of the liquid sodium is known, the 
velocity of transversal wave and _ rigidity 
modulus could be evaluated from the above 
equation, since the velocity of longitudinal 
wave in liquid sodium at the melting point 
was measured in the present experiment. 
Values of v, calculated by two independent 
methods are tabulated in Table II for a number 
of metals. 


Table II 
Li Na K Al Gael sd 
E 10-18 partition function 366 Dees pe 3.9 ee 1.6 
1 i 4 | | 
vs exp (L/3RTm) 5.2 2 0.aae 216 3.9 Oo hems 


The upper row is the value determined by X- 
ray diffraction patterns, and the lower row 
is that calculated from Drsyr temperature 
and equation (11). Good agreement was ob- 
tained between the two values. For sodium, 
we take 2.110! for vz, arithmetical mean 
yf the above two values. 

Putting p=0.9287gr., M=22.997 gr., N 
= 6.023 x 1025, vp=2542 m/sec, vz=2.1x10" in 
12), we get vs=1032 m/sec., as the velocity 
yf transversal wave in liquid sodium at the 
nelting point. From vp» ond vs above, the 
‘ollowing quantities are obtainable. 


_ Poisson’s ratio; o=0,402, 
bulk modulus; K=4.70 x10" dynes/cm?, 


rigidity modulus; =0.99 x 10! dynes/cm?. 


t should be noted that the above moduli are 
10t true elastic moduli, but of complex 
sharacter owing to the visco-elastic nature of 
he substance. But, both the experimental 
esults described in §2 and the above treat- 
nent imply that the apparent rigidity modulus 


does not vanish to zero at the melting point, 
but remains finite. Therefore, it is very likely 
that the impossibility of propagation of shear 
wave in liquid is not due to. the vanishing 
rigidity, but to extremely large absorption of 
wave energy. This result may be useful in 
interpreting the opacity for S wave of the 
earth’s outer core. 

§4. In the present section, some consider- 
ations will be made in connection with the nature 
of the mantle-core boundary. Buen (1949) 
has subdivided the lower mantle D into two 
parts, i.e. D’ and D” on the basis of the 
seismic wave velocity. The region D” is 
between the depth of 2700 and 2900km, just 
outside the core. In the region D’ (residual 
lower part of the mantle between the depth 
of 1000 and 2700km) the P and S velocity 
gradients are positive and almost constant. 
Thus, we can conclude that the region is 
homogeneous and chemically uniform. 

On the other hand, Jerrreys and GUTENBERG 
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agree on the smallness of the P and S velocity 
gradients in the region D’ as illustrated in 
Fig. 7. Some other authors suggest even 
negative velocity gradients in this region. 
BuuLen has interpreted this phenomenon as 


Fig. 7. 
of the mantle. 
by Jeffreys, and dotted line is that of Gutenberg. 


Velocity distribution in the lower part 
Unbroken line is those given 


due to either (a) a sharp increase in the 
density gradient or (b) sharp decrease in both 
the AK and yw gradients in the region D”. 

In a homogeneous region, we have 


do Nee nt dp _ Gmo 


4: 
o= i fos == ’ 5 
Fe de Semi g e ire ar ot) 
Hence, 
dK _ OK 90 =14?: 00 
dp Op “Op 0.00 
dO dO 
=] fax ts -1 
a Ai l+g9 pA (14) 
Table III 


Depth 


oo) | (c.g.s.) 
2500 2.2x 1038 Se 
D' 2600 Po) DRE 
2700 232 Oe 
2700 0.0 x 103 1.0 
Lge 2800 0.0 1.0 
2900 1050 1.0 
2900 | 2.0108 = gure 
1H, 3000 20 Sal 
| 3100 Ziel a2, 
LO LE 


Using values of gy and @ in his model A, 
BuLuen obtained the results tabulated in Table 


Ill, on applying (14) to: the depth between 
2500 and 3100 km. 

It must be noticed that the values of aK ap: 
in the region D’ are very small compared! 
with those in D’ and E. This implies that 
D” is different from D’ in some important 
physical aspect. It may follow that the region) 
D’ is not homogeneous, and there is some: 
change in chemical composition or polymorphic: 
transitions. 

Using his second earth model, model B, | 
assuming the compressibility-pressure hypothe-- 
sis, BuLLEN (1950) concluded that the density? 
gradient in D” is about three times that i 
D‘’. This implies the accumulation of a denser: 
matter in the region, the change of density? 
inside D’’ being at least 0.3g/cm*. Since: 
BuLLEN’s compressibility-pressure hypothesis, 
from which his second earth model is derived,, 
is based on the assumption that A and 4: 
chane continuously in the region between 0) 
and 5370km depth, the smallness of the P 
and S velocity gradients in D’ must be due: 
only to the change of density gradient. 

As for. the sharpness of the mantle-core: 
boundary, EuckEen (1944) said that the data: 
on the reflection coefficients of seismic waves: 
does not seem to permit any other inter-. 
pretation than a genuine boundary between a: 
solid and a liquid. On the assumption that! 
the core is in a melted state of denser material! 
such as iron as usually conceived, it is very: 
probable that in the region D’’, just outside’ 
the core, the liquid core phase would nucleate 
as a pre-transition phenomenon, The density. 
of the core material is greater than that of 
the main material composing the said region. 
The nucleation of such liquid core material 
would also make K and yz small as described 
in §2. If this is the case, the observed small- 
ness of the P and S velocity gradients in D’” 
may be due to both (a) sharp decreases of K 
and fe as in the present experiment and (b) 
the accumulation of denser matter. There- 
fore, the nature of the mantle-core boundary, 
in the author’s opinion, is associated not only 
with the accumulation of denser matter as 
discussed by BuLLen but also en the marked 
decrease of AK and yz. 
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‘ Results of high pressure experiments so far 
made show that the difference between the 
compressibilities or the bulk moduli of different 
substances tend to diminish as the pressure 
increases. 

Such experimental edivences seem to be in 
favour of BuLLen’s conclusion. But, there 
is no reason why the bulk and rigidity moduli 
ought to have the monotonic increase in the 
region D’’. More complete experiments are 
needed including the measurements of shear 
wave velocity at the transition point. 
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Abstract 


New general differential formulae for the change of the geographical coordinates system 
are derived under the simple assumptions that the datum point and any other triangulation 
point are fixed in space and that the small quantities such as products, squares and so on 


of da, de®, Ogo, OA, Sag etc. are all negligible. 


The formulae for dy, 6A, dh are expressed 


in terms of geographical longitudes and latitudes, including those given in the text-book 


of HELMERT and of BONFORD as special cases. 


With respects to 6s, dai, dan, the usual 


formulae are available as they are, provided that d¢s’, dds’ are previously evaluated by 


the present formulae. 


The astro-geodetic conditions for the azimuth and the zenith distance are derived, and 
in particular the Laplace condition corresponding to a special case of one of the astro- 
geodetic conditions is examined from a general view point. 


1. For the changes of the adopted coordi- 
nates of the datum point and the reference 
spheroid, several differential formulae have 
hitherto been given. We meet with ambiguity, 
however, in applying these formulae to such 
problems as of connection of different triangu- 
lation systems and of determination of the 
earth’s ellipsoid. For, aside from a question 
of approximation, the respective coefficients in 
the formulae derived from the different princi- 
ples are not necessarily equivalent to each other 
as readily seen by comparing the formulae in 
Hetmert’s text-book (1880) and those in Bon- 
rorD’s (1952). And further, the minor axis of the 
reference ellipsoid is defined, as usual, to be 
parallel to the rotational axis of the earth, 
but the reference ellipsoid postulated at first 
does not necessarily satisfy this definition be- 
cause of the observational errors and the 
deviation of the vertical at the datum point. 

Taking the facts into account, I have derived, 
from a quite general point of view, the new 
differential formulae, which are accurate so 
far as the first differentials are concerned, and 
which are expressed by the geographical co- 
ordinates not in the forms of infinite series but 
of finite terms. Assumptions in their mathe- 
matical treatments are simply that the datum 


point and any other triangulation point are fixed 
in space, and small quantities such as products, 
squares and so on of the first differentials 
Oa, 0e?, 02, OG, Gay etc. are all negligible. 

2. Let two systems of the rectangular 
coodinates, the old and the new, be O-XYZ 
and O-X’ Y’Z’ and the corresponding axes 
deviate from each other by small angles 
(Fig. 1). Regarding these angles as resultants 
of small rotations 40z, 40,, 40. about the axes 
OX, OY, OZ respectively (clockwise rotation 
seen from the origin is assumed to be posi- 
tive) and neglecting second orders of small 
quantities, the direction cosines between the 
two systems are given as follows. 


| x, Y Z 
x! 1 40, — Ay 
Ye — 46, i 40x 
Z! | Ab, — Az 1 
Now, put 


A0z=Aacos¢%, 40y=4¢, 40:=42+4asin go 


and denote the coordinates of the old origin O 
with respect to the new system by Xo’, Yo’, 
Z). Then, transformation of the coordinates 
beteween two systems are expressed by 
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Fig. 1. 
X’-X)’ 1 (4A+4asing) —dAg xX 
Y’-Y,)’ | =| —(4A+4a sin ¢o) il Aa COS¢y Y: (1) 
Z’- Ly AG —Aa COS Py 1 eZ 


It is noted that a small rotation 4a around a line X/cos = Y/0=Z/sin ¢) is equivalent to 
40,=Aa cos $y, 40y=0, 40.=Aa sin Yo*. 


3. The datum point P) and a point P representing any triangulation point being fixed, 


each coordinate system is assumed to have a respective reference ellipsoid (spheroid) such as” 


for O-System aN ls a =1., 
@ a a*(1—e?) (2) 
Py aes ai ? cos (A—Ay) , y= g w id sin(A—Ap) , — ae sing, 
p Ze 
for O’-system a eas A . 
on 2 a?(1— ae he (2) 
,_@ COS r . e a 
x se cos(a’—Ap) , y= z oes sin(a’—Ay), 2’ oe) sin gy’ 


Let the normals from Py and P to both the ellipsoidal surfaces be PoQo, PoQo’; PQ, PQ’, 
respectively, their lengths be ho, h’y; h, h’, and for the sake of simplicity let the line P)Qo 
lie in the plane XOZ (Fig. 1). The coordinates of P) and P become 


P 4 X=2)+ho COS Po, Y=y%=0, Z=z+hy sin ¢o, ' 
0- 
X"=29' thy’ Cos Yo’ COS (Ao Ao), YY’ = Yo’ +ho’ COS Go’ Sin (Ap’ —Ay), Z’ = 20’ tho’ sin gy’. 


* The same holds even if the line may be X/cos (v)+0¢0)= Y/O=Z/sin (ep 8¢0) whithinan order of 
approximations admitted here. 
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P- { X=a+hcos ¢ cos (A—Ap) , 


Y=y+hcos gsin (A—A,) , 
X’=2' +h’ cos g’ cos (A’—A,), 


Z=z+hsing, } (4) 
Y’=y' +h’ cos g’ sin (A’—2)), 


Z=2 +h’ sing’. 
Putting 


=a+0a, e?=¢e?+ de? ; ho’ =hy+0hy , h’=h-+0oh, 


: , (5) 
g =9t+0¢, M=A+O2 | 

Ne assume that the quantities with 6 and 4 to be so small that their products, squares and 
© on can safely be neglected. 


4. At first, let us determine the position of the origin O with respect to the new coordi 
lates system. Substituting (3) into (1), we have 


Yo =G0t 9G, A p=Ay td ; 


xo +ho’ COS Go — Xo’ =(ao+hy COS Yo) —(Zo +My Sin Y)AG , 


Yo +ho’ COS Go’ 0Ay— Yo’ = 


—(ap+ho cos Gy) 4A+{(2o+ho SIN Yo) COS Gyo—(ap+ho COS Yo) SiN Yo} 4a 


Zo tho’ sin Go —Zo’ = (Zo +ho Sin Go) + (ao thy COS Yo) AG . 


Tence, 


Xo = (ay —20) + (Ao COS Go’ —h COS Go) +(Zo+ho SiN Go) AG , 
Yo’ =Yotho’ cos Go’ O29 + (wo +ho COS Go)4A —(Zp COS Gy —2%q SIN Yo) da , (6) 
Zi = (Zo —Zo) + (hy sin Qo —hy sin Go) — (ap tho COS gv )Ae : 


t is evident that 


a( as sin? 2 502 24% Sin Y COS? 5, 


W. Were ws 
ee a smti-sones tine cos (A—Ay)ae# 
“s+ ) sin g cos (A—Ay)ade gs sin (A—2))add , 
ts 
y’ —y=dy= a} gee Se sin (A— Aydap Sn Teese sin (A—A))ade? a 
a ) sin g sin (A—Ay ade asa cos (A—Ay)ad2 , 
2x —z=0z= eS on o} : =a sin esa (1 + sae aes a COS GANG . 
n a special case when y=, A=Ay in the above expressions, we get 
wy! y= 78 Ba a Pe tate ade—! ee sin goa0%o , 
Yo y= ada , (1) 


1 
Zo SS 


—e . SiN Yo COS PO Ve sone IG 5 
wetted pe ie aoe? + COS Pod Po « 
W SiN oa 2Ws ( WwW? ) Wo? poegy 


ntroducing (7’) into (6) and referring to (5), we have 
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2 h 
Xy/=( ae + Hu) COS Yo— pe A SS We sam) 2 ) SIN Ppa Go 


+") n gad 
(ape We sin goad , 


Yi =(qy+ +) COS Gpa(OAy+42a)+ oa SiN Gj COS Yoada , (8) 
sin cos? " 
Zi =( 2 al sin Yg— 5 = (1 + aa ) ave 
+(e fo “cos YoadGo— (7+) cos Pade . 
Wo we Wi 


These formulae are useful for estimating how far the center of the old reference ellipsoic 
is located from that of the new one as discussed by K. LeprersTecerR (1955). It is to be 
noted here that even if the reference ellipsoid is kept unchanged, its new center canno 
occupy the same position as before when the coordinates of the datum point are changed 
For example, even in such a special case where da=de?=0, 4g=0¢, 44=—OAy, da=0, we 
have X,/=0, Yo’ =0, Zo’ = —(e?/ Wo?) cos? gp0Go, thus Zo’=<0 in general. 


5. Now, the required differentral formulae for the change of the geographical coordinate 
systems will readily be obtained. Substituting (4) into (1), we get 
x’ +h’ cos ¢’ cos (A’ —Ay) — Xo’ =a +h cos ¢ cos (A—Ay) + {y+h cos ¢ sin (A—Ay) }(4A+4a sin go) 
—(z+hsin ¢)4¢, 
y’ +h’ cos ¢’ cos (A’ —2))— Yo =y th cos ¢ sin (A—Ay) +(z2+h sin v)4a@ cos Yo 
—{x+hcos ¢ cos (A—Ay)}(4A+ 4a sin go), 
2’ +h’ sin gy’ —Z) =z+hsin g+{x+h cos ¢ cos (A—Ay) }49—(y +h cos ¢ sin (A—Ap) } 4a Cos Go . 


Referring to (2), (7) and (8), the above equations are reduced to as below. 
1-2" kh\ . eh : 
-(ao+ ‘) sin 9 cos (A—A,) ade = wr g eos g sin (A—Ap)a(dA +42) 


da 


+cos ¢ cos(A—A))Oh= \- W608 gy cos (A—Ay) ‘wes ie cos vo} 


sin? g cos 2. Sin? Qo C ee Po ah 
ed pase ee SUS P net) ekg e de2— uy 


l-e h\. l—e h : 1 , : 
+\-(Fe+G _ o+( w,t ) sin sire 7 COSY SIN Gp SIN(A—Ay)ada@ , 


l—-e h\. ; i. We 
-( ——— - )sin g sin (A—A)ado-+( wr =) cos g cos (A—Ay)a(dA+ 42) 


; Oa : 
+cos¢ sin(A—A,)dh= — woes gy sin (A—Ay) — sin ee ? sin (A—2) ade? 


WI bey oe: ‘ 
+R W, “|. ~ ) COS Pod (OAy+ 4A) + (tg isin Y COS Py—COS ¥Y SIN Yo COS (A—Ay) | 


, sing Sin @o 
—e? cos onl = hase : 
W Wo 


Lok Le". 
( We ait 1 es v0g-+sin voh= \- = sin goa + (7p, a+ ha) sin vo} 
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sin g Cos? y SIN Go COS? @o > h 
rapes 1 assy is tes a 2 9 
: ( an : )- OW, (4 3 fader + (7 We a a) 608 $Y) A0gy 
DEMS <5) GET Ls hee 
+4( aS 7 £08 #608 (2 Ao) Gat 1 60s vo bade 


Lepahk . 
aa: =) COs ¥ COS gy Sin (A—Ay)ada . 


Now the three unknowns dg, 62, 6h are easily determined from the above equations. In 
act, 


Byes) C= ‘ca Be [cos ¢ Sin Yy—SiN Y COS Yo COS(A—Ay)] 
0 


“al ba i; eae eos g =i, +| 9 a sin g cos g— OW a 7 SIN Yo COS Yo[COS Y COS Yo 
0 


+sin ¢ sin gp cos (A—Ay)] lor l(et!) [cos g COS Yo +SiN ¢Y SiN gp Cos (A—Aj)] ap 


W.° 


+(qty h cos a y= (w+ W, +2 )ieos Y COS Yo+sin Y SiN Yo Cos (A—Ay) 


—¢ en We ing cos(i— 2) bg iat sin ¢ cos gp sin (12) (ay + 42) 


Ww Wo Wo 
of” Rae oa) - heos ¢ sin (A—A,)da (10) 
H-Grta)t4( Wei 
ae x 1 da, oho in (A—Ap) Si? Po. COS 0 in (QJ. 002 
aaa — (Gi a tg) 608 #e sin A=) SEE sin (A— 2a) 


h dea BN ae fee 4 A—Ay)d 
+(H Ws i 7.) sin go sin (A— 1) 890+ \( W +7) sing tes 5 SIN Yo¢ Sin ( Ay 


+( a ho wt.) cos Yo COS (A—Aj) (049+-42)—{(F-+7 Jlcos y Sin gy—SiN ¥Y COS Yo COS (A—Ay)] 

0 

+(e ooo) COS Yp COS (A—Ap) | dex, (11) 

0 
sing SiN Go 

oh= \- este + 0h sin Y SiN Yo+COS —Y COS Ya COS (A—Ag)] +e? Lene aor ) sin goal 

as sin oO Se aera Por sin g COS Yo—COS ¥Y SIN Yo COS (A—Ay)] haves 

2, 0 0 


1—e? , ho ) [sin g COS ~y—COS ¥Y SIN Gp COS a—4)]} A0Qo 
a 


: 2 sin ¢ sin ¥o — 0 
—s +) [sin g COS Yy—COS Y SIN Yo COS (A—A))]—e ( we Wy, ) cos g cos (a a) ha g 


sing SiN Qo 
+1(54 - +f ) COS Y COS Gp Sin (A— 4) a(dAy+ 4a) — ae SC) 
12 
X COS Y COS Y Sin (A—Ay)ada . (12) 


nese expressions, however, may be simplified by.using the following notations, 
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A = sin gp Cos Y—COS Gp Sin Y COS (A—Ap) , B=sin 9 sin ¢+cos gp Cos g cos (A— Ap) , 
B’ = cos gp cos y+Sin gp sin gy cos (A—Ap) , C = COS Go SiN Y—SIN Go COS Yo COS (A—Ap) , (13) 


: __ sing SiN go 
D= cos @p sin (A—Ap) , E= WW, k 


Geometrical meanings of the notations may be inferred from the fact that these are re- 
duced to the following, when the old ellipsoid (spheroid) turns out to be a sphere. 


A=sino cosa, B=cosa, C=sinocosa,, D=—sinosinay, E=sing—sin go. 


By the aid of (13), (9) (10) (11) are rewritten as follows 


l—-e@ h - 1 da, bho 4 0a os sin g cos Bi go OS ool aes 
ae + i) ee = A(= a ;: \ +28 cos oT 5 COS et Ws Ws 


—e n\ » Beek: ~ -(4. te.) ‘Mae , 23 la 
Gr + He B dgot |( apt Geos ( Ao) wt a B’—e’E sin ¢ cos ( 0) r4e¢ 


-{( 7 ee »\D a (Bay +.42)-+ ~(yta)t eE sin o| Baas 10 


il h 1 da. bho Sin? Qo 3 
=) cos @ (6A-- ASD | ee ey) Sp oe 


: 1—e?, h\ . 1l—e? h : ‘ 
(ee ho “sin Yo sin (A— —4)} OGot \(“yr + =) sin o-( _ +") sin eatsin(A—A)de 


aie he a COS Gy cos (A — 4) + 42) {pt g Ate Boos @o COS (A) Le ; ai 
Wo a Wa 


a, E 
ae A ae Esin ial +} 3 


oh -sin g— SIN @p COS voh ade? 


2W.3 
oat LOOK} 
) ade +4 — oe = i.) + Ecos ¢ cos (A— i) bade 


Wo 


H ibe In "YD cos va(d29+-42)—eED cos pada. (12/ 
0 


The expressions (10)-(12) or (10’)-(12’) are the new general differential formulae. 


6. The formulae obtained in the preceeding paragraph are considered to represent the relatio1 
between the datum point and any other point. As is seen in §2, the datum point play 
a special role in the transformation of the coordinates, since the rotation axes of 4g and de 
have been defined with respect to its geographical coordinates. Hence, for arbitrary pair o 
triangulation points, these formulae cannot be available. We are going to introduce them. 

Solving the equations (9) by regarding dg, 649, Sho as unknowns, we get 


1l—e?_ hy 1 da. oh Oa (E sin g cos g 
dy ot eee 23 pee ees 
( ot =) Qo =C( ata) eB COS Yo {9 608 #0 +B Sic 


SIN Pp COS Po 3 fi—e. h a ( =) ( ee : 

a Po COE SaaL| bc2 2B (ee) pines ay ed 2 
a 2W,3 er ( ws + Z| = wt He if ee Esings} de 
+( a s) cos.p sin vo sin(A—ApOAT AA) (ee 7 808 ¢ Seay 

W'a Wa rere 


1 fo m/l Og. On sin? g cos¢ .. 
Ge iB ; ) cos pride dd)=(45 r — 7, ) 698 gy sin (A—A,)— ops sin (A—Ay) de? (14 
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“Ci 
ie +2 )e —e°E cos go} de, 


1 9a, oh 1 Oa singcosg  E . ) 
2 Saat 2B 8 : ‘ 
ae. +" )- (Gyte sin ¢) \A 2 Ws 4 9 pr ao eek 


42) sin y sin (A=A,) do+( pts) cos g cos (A—Ay)(BA-+.da) 


“paQie a 1 )ee+4 al; : +7) aE cos val de —( vty)? cos 9(84-+42) 


Ww Ww Ww 


In these formulae, the coefficients except those for 4g and 4a are essentially the same as in 
the formulae (10’) (11’) (12’), since the notations defind by the right-hand side of (13) are 
replaced in such ways as 


A-C, B-—B, B’=>B’, CoA, D— —cos ¢ sin (A—2)) , 3p 


as the suffixes are changed. Hence, it is sufficient for the present purpose to revise only the 
coefficients of 4g and 4a in the formulae (10) (11) (12). 

As the differentials for the datum point dg, 04), day) are expressed by those of the first 
triangulation point 0¢,, 64,, dh; using (14) and those for the second triangulation point dgz, 
OA,, Oh, are expressed by dg , dA9, Shy by means of (10’) (11’) (12’), the formulae giving the 
relation between P; and P, are now obtained as follows. 


1l-e h, 1 da. oh, 0a Px Suna 
(pS teen aa( ft) ce tac, Ms Pomentfs 


— Byy/SiN P1008 G1) 5604 { B (toa 
Br WwW f é a 21 W 


au {( = pa hs 4 cos (A,— (+2) [cos y, COS Y2 COS (A — Ag) +SIN G1 SIN G2 COS (A2—Ay)] 
2 1 


cee (2 W 5 ) 
6¢:—4( ——+-— ) Da sin g2¢ (0A +42 
= - Pi =e 218 oak 1 ) 


—e? Ex, sin g2 COS (.— 2) |e a {| - C= hs 2) +e? Ey; sin os feos go sin (A, —Ag) 
2 


+(5y,+ i )(cos 91 COS G2 Diotsin ¢1 SiN G2 Doo—Sin Yo SIN ¢Y2 Dz) 4a , (15) 
1 #2 


il h, 1 0a a) sin? ¢1 é 
ae =—]),,{ —- — —D = + j¢ 
( = *) cos g2(02, +42) ( ce + he 21 OW! € 
a 


—e? . : 1—e fp .. 
ioe +i) sin ¢g; sin (s—2)} do. +4( oe ab. >) SIN @2 


1 hy 
lew. i +7) sin oi} sin(s— 2) de+(+ +) cos ex cos (4s Aydt AY) 


il 


—| Aw (p+ hy 2) +eEn cos 0008 (da Aa)t (p+ [eos Yo SiN G1 COS(A,—Ap) 


+sin Gp COS ¢1 c08 (22a) (16) 


0a 0a 
W, Wi 


heal in 
ae Cu (uae +" Vade, sia {( W, +7) De cos esha a(di,4+ 4a) 


; (Sate Cy ‘ y ip 
Oey \- TO + Bal +h, )+6°Bs sim vba +45 SIN P2— 9 Ws pr Ce on ase 
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ea dCi cosa sin eee nee at ge 
aly Woo [sin g1 COS Y2C 2— Ao P1 2 1—4o 2 Waa 2 
+e? [E1o(Cap Sin Yo —Bao COS Go) | +e? E10 (C20 SiNGo — Buy COS Go) + Ex9 COS G2 COS (Az—Ay)] lade 


HGnt EY Dw COS 91 SiN Y2—Dzy Sin Y1 COS Y2+SiN Pp COS Y1 COS Y2 SiN (Az— Ao) 
5 
—e? Fy, Da) COS sada : (17) 


where 


Aj: =sin gi cos gj—cos gi sin gjcos(4j—4:), Byu=sin gi sin y;+COs Yi COS Yj COS (A;—A1), 
Byi’ =cos gi cos g;+sin gi sin gj cos (Aj—Ai), Cys = Cos Yi SIN Yj—SIN Yi COS Y; COS (Aj—A:1), (13’) 
sing; sing: 
Ex= — : 
ae: Wi 


Dy, =cos gi sin (Aj—Ai), 


7. As for the differentials concerning a geodetic line, namely, the differential change of 
its length 6s, that of its forward azimuth da. and that of its backward azimuth da, the 
usual formulae as in HeLtmert’s text-book (1880) can be available as they are, if dg’s and 04’s 
are evaluated by means of (10)-(11) or (15)-(17). It is needless to say here that a geodetic 
line is defined by two end points Q:, Q. which are projections of a pair of triangulation 
points P,, P., on the reference ellipsoid. Taking the eastward longitude as positive and re- 
ckoning the azimuth from the north, HeLmert’s formulae are written 


l—e?) . dm a(l1—e?) .. 
Barge oe CO eee aa OR Beh & len) 8 waiters 
ue We 2(0A,—021) + AWA Pas RET eee te ae 
__ Se? cos 91.C0s #2( sin Q1 COS @1,——>- Sin ayy cosa 
2(1—e?) WW, a om ses eae B 
1—e?) . a(1—e?) . dm 
0Ay= COS 1 COS@}2(0A,—0A alle") N20 ao : (eo) OQ, 
mW ak a mW ait mW;$ Tent S PEACE ee i 
_ Se? cos 1 Cos at sin ayy COS a1 —>- sin a2, cosa ra 
2(1—e?) WW, 12 21 ae 21 COS «) 4 
_ 9a, aV1—e a(1—e?) a(1—e2) 
OsS=s laypae COS $1 SIN Ay2(02,—02,) — — We COS @12,091—— We COS A210 Gy 
be ade? 
2(1—e) 


In most text-books, each of the above formulae is given as a result of combination of two 
sets of expressions, the one for changes of dy’s, 64’s and the other for changes of 6a, de? 
(or 0%), but it is possible to derive (18) at the same time by taking into account simultaneous 
changes of 6y’s 04’s, 6a, de, although this is self-evident and such treatments are con- 
siderably complicated. The reason why the formulae (18) are valid irrespective of orientation 
of the ellipsoid is evident. As these are derived by supposing an auxiliary sphere with 
radius a, on which the reference elliposid has been projected, the orientation of the ellipsoid 
does not matter explicitly. In this case, however, its indirect effects should be included 
in terms of dy’s and 04’s by using the formulae (10)-(11). 

i a way, a formula for 4a is obtained by elimination of dh) from two equations (107) 
and (11’) 
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MLS ANG | oi l—e ho l-e 
Ge ya B*)—e@ EC £08 v» [dae= {( Ws oes 7) 008 #804 (Cy + )eos eo de sin (a—2o) 


ho 
+Al(— We a nos o COS (A—Ag) (Blot Ai)—( +7) cos (i442) 


ecg 1 h co é 
+} B( we : )+( W, le 3) —eEsin eo}cos g sin (A—2£))dg+e2ED cos oo 


in3 in3 
+ DIE +e ( a We) cos gde* (19) 


8. Let us reduce the general formulae (10)-(12) and (15)-(17) in some special cases. 

(i) A case when the coordinates axes of the new system are parallel to those of the old 
system. Referring to (2), this corresponds to 4g=4i=4a=0. 

In this case, the formulae (10), (11), (12) coincide respectively with (15), (16), (17) and are 
-educed to 


—p2 ° 
a 


W,? y Wia Z, 2W3 
/sin 2 pile ; 
— By : are (8 + ae) de—\(G+ 7) Da sin os} OFn, 
il | 1 0a. oh sin? ¢1 
0A. = — 2 u i ae ere =r) 2% 
Gate oe eae en) mye ; 
(20) 
+{(55 +%) sin ¢; sin (a,—2)} OG1 (Gta) COS 1 COS (Ap—A1) 02, , 
Oh,= 1-7 W. es Oo 4 dh: \+eBn sin eatal + +| im SiN G2 0, sin ¢, COS or} aes 
2 ff i 
2G ye hy prey f 
21 — +o a Ort) eee Dy COS es}aca i1t4a) . 


[The formulae given in Bonrorp’s text-book (1952) correspond to this case. Nevertheless, the 
xpression for N therein equivalent to the above 6h (i,=0 is assumed) is transformed into 


N=| +f —T+ at e’E sin o| a+] ww alice ar sin ) hy 


ini Mele e “U3 E sin g- aay a { ade 
W 1—e? < [Wolf D 
+L weal ane ht Le dae, cos foe 


Somparing this with the third formula in (20), it is readily seen that the differences in the 
corresponding coefficients are of orders of e. As the original literature has not been 
ecessible, I cannot now state why the formula N deviates from my accurate d/-formula. 
Phe other two expressions for dg and 04 therein have been reduced from the N-formula by 
neans of differential relations, which are correct so for as a rotation 4a is not taken into 
count. Consequently the same order of inaccuracy is also expected with respect to the 


ther two formulae.* 
(ii) The minor axes of the both systems are parallel but with a small rotation of longitude Ad. 


* In BONFORD’s text-book, a mistake is found in oe the variable from Q to U, accordingly 
. term in dy must be (1— f cos 2¢) instead of (1+ fcos? ¢) 
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) 
As 4g=4a=0, the only difference from the case (i) is that the terms with 6/1 and 0A, are 
replaced by 62,+44 and 64,+42 respectvely. If 44=—6A, is assumed, the terms including 
(84g +42) in (10) (11) (12) disappear, whereas the terms with (64,+44)=dA—6A) in (14) (15) 


(16) remain as they are. | 
(iii) Corrections for the geographical coordinates at the datum point are compensated by 


small rotations of the coordinates axes, namely 4g=0g, 4A=—OAg, da=0a. On this condi- 
tion, (10) (11) (12) are reduced to 


Gas: 7) be A( leu oe + 6E cos ot 
a 


Wo a 
E sin g COSY _ , SIN Po COS Yo be e+\(w *) cos (A— 2p) 
ae Siri yarn eae Cy hans. 


=F [cos? yo cos y—(1— e?)sin? gp sin g Cos (aa) Oo 
0 


{(tG ree sin | Dias ; 


W a 
(ga) cos (6A — 3i)=—D( 1 ee ame P02 


Wo a 2 W,? 


ae le +7) sin —— sin? sin (A—Ay) bo 


(21) 


oh= es + dh \te Esin eval +} E sin g - SiN Gp COS vo} aoe” 


2  2W,3 


w 


© "sin gy cos g cos(A—Ay) 


Ww 


\(y+ a: i) A+@E cos gy 008 (A— 4) da, 
W 

4 
al 


= as Yo Sin g+(1—e?) sin’ go cos g cos(A —4,)]}adgo—e* ED COS Y Ada, . 
0 

The formulae (15) (16) (17) should be transformed into other forms, but their detailed descrip: 

tions are omitted here. 

In a limitted case when the spheroid is kept unchanged and hy=h=0, V. R. OLanpER (1951) has 
derived very elegant formulae with a slight reduction of accuracy from Heumert’s origina 
ones. He estimated the corrections to be applied to his formulae by refering to W. K 
Hristow’s series (1942), which are accurate up to the order of o3, e262, eto so far as_ the 
principle accepted tacitly is true. The very principle in deriving the series may be 
regarded as the same as that for Hetmerr’s differential formulae, viewing from a_ poi 
that the straightforward differentials are only admitted without attention to the changes 0: 
the spatial locations. Taking the fact into account, I proceed, for comparison, to evaluat 
the corrections of his formulae on the basis of the above formulae (21), by letting h=h= 
Kee a de*=0. The results are quite the same as those given by OLANDER himself as show! 

elow 


pt 2 
(d9)= =(4 P cos* eo) Boot ( : b? cos? gy sin 1) Bets = = : q COS? Go [28G)+O04 
NY 3 9L9 Q IN 2 
(02)=( ri e?b? 1 sin 20) Bo +e*b cos? yo(—1+b tan gp —e? sin? yy)bay= —q cos ¢ bday +O, 


The formulae (21) have not been examined up to the higher orders, but the above result 
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suggest that this special case corresponds 
closely to Hetmert’s formulae. In my view, 
the formulae (21) seem to be applicable for 
connection of two different trianguration sys- 
fems. 


9. As for J. F. Hayrorn’s formulae used 
in his famous determination of the figure of 
the earth, (1909) they may not be suitable for 
getting the accurate geographical coordinates 
of triangulation points, because small quan- 
tities unnecessary for his own purpose have 
been neglected. In deriving the formulae, 
Hayrorp has postulated three processes, 
namely (i) a simple rotation 42=—dA) around 
the minor axis, (ii) translation of the datum 
point dg, along the principal meridian, keep- 
ing s and a as constant, (iii) rotation of the 
vertical section da) around a normal at 
a point corresponding to the datum point s 
being kept constant. The assumption that 
the arc length of any vertical section s is 
constant appears to be plausible. However, 
strictly speaking, since what ought to be kept 
constant are spatial positions of the two end 
points, the assumption is not reasonable. 


10. As already mentioned, products, squares 
and so on of Jda/a, de*, dgo, ---, have been 
neglected in my formulae in the same manner 
as in the others’. Hence, in order to get ac- 
curate values of dg and 64 by means of (10) 
(11) (12), we must take this fact into con- 
sideration. The reliable estimation of the 
errors therefrom seems to be impossible 
without performing lengthy calculations, but 
the effects on the final results may be 
avoided by a method of successive approxima 
tions. For example, let dg and 64 be required 
with an accuracy of 0’, 0001=5x10-", then 
V5 x 1020/2 =1.6 x 10-°=3”’ may be tentatively 
taken as a limit of da/a, de?, go, --- in (10) 
12), within which the successive approxima- 
‘ion should be carried out. 
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11. There remains, however, an important 
problem, namely, which case is the nearest to 
the reality in applying the differential formulae 
for changes of the datum point and of the 
reference spheroid. It seems to me that there 
is no apriori basis for deciding which set of 
the formulae is the best, but it should be 
empirically judged from the synthetic point of 
view, taking into account triangulation sur- 
vey, astrometric observations and gravimetric 
measurements. Investigations in this line will 
be left in future. 

In this connection, it is interesting to show 
the corrections to be applied to the old azimuth 
as well as the old zenith distance in conformity 
with the changes of the datum point and of 
the reference spheroid. For, they offer some 
important relations, namely the astro-geodetic 
conditions in which the Laplace condition is 
included as a special case. 


12. Taking an arbitrarily chosen point Pas 
an origin, consider two rectangular coordinates 
systems P-UVW and P-U’-V’W’, where U, 
U’; V, V’; W, W’ are directed eastwards, 
northwards and upwards, each with reference 
to the old spheroid O and to the new one O’ 


respectively (Fig. 2). Letting (gy, 4) and 
(Zenith) 
(North) W’ 


Big. 2: 


(y’=y+69, ’=A+0A) be the geographical 
coordinates of P referred to the old spheroid O 
and the new one O’ respectively, then 


(Of —sin (2’— Ap) cos (2A’— A) 0 x 
Vv |=[ —sin¢’ cos(4’—4) —singsin(A’—49) cos 9” es (22) 
Ww’ cos g sin (A’ —Ap) sin ¢’ ZR 


cos gy’ cos (A’— Ag) 
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xX —sin(A—Ay) —sin g cos(A—Ay) cos gcos(A —Ap) Gf 
Y }={ cos(€A—’) —singsin(A—d’o) cos gy sin (A—Ao) V (23) 
Z 0 cos @ sin @ W 


In (22) X’, Y’, Z are expressed by X, Y, Z by the aid of (1), then the relations between U’, 
V’, W’ and U, V, W are obtained by means of (23). Thus, finally 


4g cos g sin (A—A)) Ag sin g sin (A—Ap) 
. +(62+42) sin g —(0A+42) cos g 
U il : : : : U 
+Aa[sin gp sin g —Aalsin gp sin g 
--cos vp cosy Cos(A—Ay)] —COS gp cos gy cos(A—A;)] 
—Ag cos g sin (A—Ag) 
. —(0A4+4,2) sin g —d¢g+4g¢ cos (A—Ap) 
V\= : ¢ 1 ; V 
—Aal[sin go sing — Aa COS go sin (A—Ag) 
+C0S @p cos gy Cos (A—Ay) 
—Ag sin g sin (A—Ap) 
; +(04+42) cos ¢ d0g—Ag cos (A—Ap) 
WwW’ : : 1 W 
+ da[sin @) cos g +4a COS gp Sin (A—2A 4) 
—COS Gp Sin g cos (A—Ay) J 
(24) 
Especially at the origin, it reduces to 
Lie 1 (Ap +42) sin gGot4a (d¢y)+44A) cos go a 
Vo’ | =| —(6a) +4) sin gyda 1 Og ,—de 1 (24’) 
Wy’ (0Ayp +42) COS Po —Og )t4e 1 Wo 


13. Let the direction cosines of any line passing through the point P be (/, m, n), 
(/=1+6l, m’=m+6m, n’=n+6n) with reference to the old coordinate system P-UVW and 
to the new one P-U’V’W’ respectively, then by means of (24), we get 


lV’ =1+6l=14+m{4¢ cos ¢ sin (A—Ay) + (0A +42) sin — 

+4a[sin gp sin g+CoOs Yo COS Y COS (A—Ay)|} 

+n{dg sin g sin (A—Ay)—(02+ 4A) cos g —da[sin gp cos gy —COs ¢Yy Sin ¢g Cos (A—2y)|} 
m =m+dm=m-+n{—d¢ +4¢ cos (A—’y) — Aa COS Gp Sin (A—Ay)} 

+1{ —4¢ cos ¢ sin (A—Ap) —(62+- 44) sin yg —dafsin go sin g+cos gp cos ¢ cos (A —A,)]} 
n’ =n+0n=n-+l{ —4¢ sin g sin (A—2y)+(04+42) cos 

4a[sin go COS y—COSs gy Sin g cos (A—A,)]} 

+m{dy—Ag cos (A—2y)+da@ cos gy sin (A—Ay)} 


(25) 


Writing /, m,n in terms of azimuth .9/ and the zenith distance 2 


> 


J=sin Z sin , m=sin 3 sin 7, n=cos 2 


Ae 


we have 


a 


/ = - o-.: ae : oe 
V’—l=dl=cos sin (6d 2 +sin ¥ cos 97 d_Y, 
m —m=d0m=cos &CcOS.WO & 


aN 


—sin sin 976.97, n’—n=On=—sin & 6 & 


> 
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and consequently 


sin 20 = cos .97d]—sin .o%dm , Sin 0 = =—0n,. (26) 


The required connections between 6.9% and é = are now obtained as follows. 


0.97 = Ag cos ¢ sin (A—4y)-+(64+42) sin 9+ 4a-B 
—cot 4 cos .%{—4¢ sin g sin (A—2,) + (64+42) cos y+-da- A 


+cot 2° sin 7% {06g—d4¢ cos (A—2y)+ 4a: D} (27) 
0 2 =sin {dg sing sin g sin (A—A)—(64+ 4d) cos gy—da- A} 
+cos .9%{—d¢+4¢ cos (A—Ay)—da- D} . (28) 


14. Two astro-geodetic conditions useful for comparison of the geodetic data with the 
astronomical ones can be derived from the formulae (27) and (28). In fact, if we put 


V=22, 
OA=1°—19=y SEC G , 6g=9°—g%=E , 0.97 =a"—as , 


Gl Se a =aet : A= 12, g=e!, ale 


we can write the formula (27) in a form, 


0.9% —7 tan g—cot = (E sin 9% —y cos _.7)=d¢ cos ¢ sin (A—Ay) +44 sin gp+4aB 
+cot 2 cos .97{4¢@ sin g sin (A—A,)—44 cos p—da- A} 
—cot = sin &%{d4¢ cos (A—A,)—4a-D} . (29) 
It is evident that the left-hand side represents nothing but the so-called Laplace’s residual, 
which has hitherto been regarded as null if the geodetic as well as astronomical values are 
free from errors* This, however, is not right in general, except only in special case (i) in §9, 
where 4g=4A=4a=0. For example, in the case (iii) where 4g=0¢ , dy=—0A,, J4A=0_ Yo, 
the Laplace-residual R is not zero but becomes 


R=6_Y-—7 tan g—cot 4 (E sin .%—y7 cos .)=£&) cos ¢g sin (A—Ay) —Oy tan P+0_WoB 
+cot 2 cos. {£, sin vy sin (A—Ay) +9 -8_W: A}—cot F sin {Ey cos (A—Aq)—8_%- D} (30) 


Thus, the Laplace condition does not hold in general. In this case, however, the Laplace 
residual is reduced considerably, if we put 0.97%)=7) tan g) or more elaborately 


6 Wo= tan yotcot F oo’ (Eo sin Woo’ — 0 COS Won’) - 


z imilar argument may be possible with ; 
ee oae s y P ae measurements in United States. Washington, pp. 
respect tod =, and some such condition may 73-92 


also be introducible. Haare ties Read 
1880 Hodhere Geodesie I, Kapitel 6. 
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An additional note: 


After already writing this paper, I had an 
opportunity to read LepEeRsTGER’s recent paper 
“Zur Ableitung der projectiven Lotabweich- 


ungsgleichungen.” (ZfV. 80, 84-89, 1955) |] 
understand through it that Veninc-MEINgEsz had 
derived the same formulae as (20) which cor- 
respond to the case (i) of my general for- 
mulae (10)~(12), perhaps by the similar 
mathematical treatment or by the “ projected 
method”. Consequently, my paper should have 
been revised in some passages. I have decid- 
ed, however, to publish my paper in its original 
form, because the Veninc-MEINESz’s paper as 
well as LEpDERSTEGER’s one “Die translative 
und die projektive Methode der astronomischer 
Geodasie.” have not yet come to my hand. 
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Observations of Near-by Microearthquakes 
with Ultra Sensitive Seismometers 


By 
Toshi Asapa 
Geophysical Institute, Faculty of Science, Tokyo University 


Summary 


In order to see whether or not there is a lower limit in the magnitude of earth- 
quakes, in other words whether or not “minimum ” earthquakes exist and to investigate the 
rate of occurrence of microearthquakes, the present writer has designed and constructed 
Ultra Sensitive Seismometers of short period, which have a magnification of 1.3107 at 
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20 cps. 


microearthquakes were recorded in 48 hours in total. 
thousand microearthquakes are occurring in and near Kanto district every year. 
microearthquakes are not aftershocks of a particular large earthquake. 
Occurring in the area within 100 km from the Tsukuba observatory. 

Their earthquake-motions have a prevailing frequency of some tens 


10!° ergs to 108 ergs. 


They were temporarily installed at Mt. Tsukuba, Kanto district. 


About 400 
This shows that about one hundred 
These 

Most of them are 
Their / are about 


of a cycle per second, and so they can be recorded only by seismometers extraordinarily 


sensitive for 10 cps to 100 cps. 


Tens of sensible earthquakes are being felt at Mt. Tsukuba every year. 


Based upon 


this fact and the number of microearthquakes recorded at the place during 24 hours, it 
can be concluded that the equation for magnitude frequency relation, 


log N(M)=a-—bM , 
holds for earthquakes of M from less than zero up to about 5 occurring in a certain 


seismic region. 
with that of earthquakes of M 4 or 5. 


Contents 
§ 1. Introduction 
§ 2. Seismometers 
§ 3. Seismograms 
§ 4. Frequency distributions of the maxi- 


mum amplitude A, and the magnitude, 
M 

§ 5. Energy of microearthquakes 

§ 6. Geographical distribution of the foci 
of microearthquakes 

§7. Résumé 


$1. Introduction 


In Japan, the routine seismological observa- 
jon has been and is being made with com- 
yaratively less sensitive seismometers. As a 
natter of fact, these instruments are appro- 
yriate for the purpose, because many earth- 


The space distribution of foci of microearthquakes agrees approximately 


quakes take place in and near Japan which 
are large enough to be recorded by them 
with the frequency much higher than in any 
other areas of the world. For instance, 
several tens of near-by earthquakes are felt 
by people in Tokyo every year. On the other 
hand, the seismological stations situated in 
other areas of the world which are seismically 
less active are usually equipped with seismo- 
meters which are sensitive for earthquake- 
motions of a few seconds or longer in period. 

The predominant period of seismic waves 
due to near-by earthquakes are shorter than 
0.5 sec., and also those due to very small 
earthquakes are still shorter. The seismo- 
meter which is sufficiently sensitive in a fre- 
quency band higher than 10 cps will be the 
only instrument that can profitably register 
microearthquakes having hypocentral distances 
less than several tens of a kilometer, 
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Most kinds of electromagnetic seismometers 
now in routine operation do not appear to be 
sufficiently sensitive in such a frequency band 
of seismic waves. For instance, the response 
curve of the Galitzin seismometer has its 
miximum at a period of 14 seconds. The 
sensitivity maxinum of the Benioff short 
period seismometer lies at a period around 1 
second, and its period characteristic curve 
falls in proportion with the frequency of 
earthquake-motions in the period range shorter 
than 0.5 seconds. 

In order to be capable of recording micro- 
earthquakes having a small £, the sensitivity 
of the seismometer must be higher than 
10° mm/kine. (or mm/cm/sec) In another ex- 
pression, a magnification of 1.310" at a 
frequency of 20 cps is needed for recording 
near-by microearthquakes. Only the seismo- 
meters having a sensivity of about 10° mm/ 
kine or higher deserve to be called “ultra 
sensitive” for near-by earthquakes. Not 
even such an Ultra Sensitive seismometer is 
sensitive enough to record microearthquakes 
with hypocentral distances more than about 
100km. So far as the writer is aware of, 
no observation of near-by microearthquakes 
by Ultra Sensitive seismometers has been 
made thus far. The purpose of the experi- 
ments with Ultra Sensitive seismometers to 
be described in the following pages is mani- 
fold. 

i) One of the most important aims is to 
explore what will be the magnitude of the 
smallest recordable earthquake and to see if 
there is a minimum earthquake. What will 
be the energy of the minimum earthquake, 
if it exists at all?” Aftershocks occur 
generally in shallow surface layers, a few 
kilometers to 10 or 20 km in depth. Small 
aftershocks having energies 10” ergs or more 
can easily be recorded by the seismometers 
with a sensitivity of 10‘mm/kine for a fre- 
quency band from 10 to 100 cps, because in 
this case the hypocentral distances are always 
smaller than a few tens of a kilometer. 

It has not yet been known if microearth- 
quakes with such a small energy also occur 
in seismic regions where larger earthquakes 
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of magnitudes up to about 6 occur “ station- 
arily.” 

The upper limit of energy of an earthquake 
is considered to be about 10° ergs (Tsuso!).” 
The existence of the lower limit in the energy 
of natural earthquakes is not ascertained. I 
the limit does exist, it might be too small tc 
be detected even with our Ultra Sensitive in- 
struments. In other words, groundunrests or 
noise in the instruments are so large that 
smallest natural earthquakes may easily be 
obscured by them. 

ii) A magnitude-frequency relation giving 
log N as a linear function M holds for earth- 
quakes of M=6 or more (GuTENBERG, RICHTER,* 
Tsusor,” and Kawasumr?). The magnitude 
frequency relation for minor earthquakes has 
been studied by Suzuki and Asapa.”? In 
both cases, the relation can be expressed by 
the equation,” 


N(M )dM =const. x 10-°*dM , 


where, N(M) is the number of earthquakes 
between M and M+dM in magnitude, and Ll 
a constant. No physically plausible explana. 
tion for this frequency distribution of M has 
been given. 

Is the magnitude-frequency relation fot 
microearthquakes of MM zero or less also ex: 
pressed by the same equation ? 

If the rate of occurrence of microearth 
quakes in a certrin seismic region is known 
can the frequency of larger earthquakes of 
M 4 or 5 in the same region be estimatec 
according to the same equation ? 

These problems are important concerning 
the magnitude-frequency relation. 

ili) If the earthquakes exist, the energy o: 
which is so small that they can be recordec 
only by an Ultra Sensitive seismometer in 
stalled close enough to their hypocenters 
space distribution of their foci and tims 
series of their occurrence are the next pro 
blems to study. 

iv) The predominant periods of earthquake 
motions or their spectrum depend on_ th 
magnitude (GuTENBERG and RicuTEr,® Asapa, 
Axi,» Honpa™ and Ito, and KANAI, e 
al.‘), and our experiences tell us that th 
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predominant periods of the microearthquakes 
are remarkably shorter than those of larger 
earthquakes. This phenomenon should be in- 
vestigated more precisely. 

Some preliminary technical knowledges are 
needed before the observation by the Ultra 
Sensitive seismometers are to be started. 

For instance; 

1) Informations regarding ground unrests 
of periods between 1 sec and 1/100 sec are 
necessary.? Microseisms, however, do not 
disturb the observations of microearthquakes, 
because microearthqakes do not involve the 
waves of periods longer than one second. 

2) The preliminary informations regarding 
the spectrum of the waves of microearth- 
quakes are also indispensable. The studies 
hitherto carried out of the spectrum of very 
small earthquakes indicate that the pre- 
dominant frequencies in their waves are more 
than 10 cps. 

3) The seismometers must be designed so 
as to have such a characteristics as stated 
above. Studies of ground unrests, or the 
short period disturbances, and of techniques 
of the seismometer designing which is Ultra 
Sensitive in a band between 1 cps and 30 cps 
have been much improved along with the 
recent progress in the explosion seismological 
research of the earth’s crust. The frequency 
response of the seismometers to be used for 
the microearthquakes observations is required 
to cover a frequency range of lower cycles 
than the instruments for the seismic explora- 
tion and also a range of higher cycles than 
the seismometers for explosion seismological 
observation of artificial earthquakes (N. Dan"). 

The short period disturbances have their 
energy cencentrated in a range between 1 cps 
and 10 cps. Therefore, a velocity seismo- 
meter having a flat frequency response bet- 
ween 10 cps and 100 cps is appropriate for the 
observation of microearthquakes. 


§ 2. Seismometers 


The purpose of the present study is to in- 
vestigate the microearthquakes occurring in a 
seismic region where the shocks of M up to 


4 or 5 are occurring “ stationarily.” For this 
purpose, it is highly important that the 
Ultra Sensitive seismometers be installed at a 
station where the ‘‘ passive seismicity’’ is not 
low. In and near Kanto district many shocks 
are occurring approximately “stationarily,” 
and any place in this district can be chosen for 
setting up a temporary seismological station 
for this purpose. However, the place where 
the Ultra Sensitive seismometer is to be 
operated must be extremely quiet. Consider- 
ing the above mentioned conditions, the pre- 
sent writer determined to set the Ultra 
Sensitive seismometers at the seismological 
observatory at Mt. Tsukuba, which belongs to 
the Earthquake Research Insitute, the Uni- 
versity of Tokyo. Mt. Tsukuba is geological- 
ly composed of granite. The seismometer 
pier is firmly cemented to the bed rock. The 
velocity amplitude of the short period ground 
unrests at the place is’ about 10-® kine or 
less, in other expression, a few Angstrom at 
10 cps, so that the order of magnitude of 
them is less than one thousandth of the 
ground-unrests in ordinary city areas in am- 
plitude. 


Seismometers 


The present writer designed and constructed 
a new electrodynamic seismometer of moving 
coil type with an amplifier for the observa- 
tion of microearthquakes. One of the simplest 
methods for obtaining an extremely highly 
sensitive velocity seismometer to adopt an 
electrodynamic transducer and a C-R ampli- 
fier. It must be emphasized that a displace- 
ment seismometer is not appropriate for 
observing very small earthquakes in which 
amplitudes are not large enough in compari- 
son with the ground unrests of short periods. 

The predominant period of microearthquakes 
becomes shorter, as the M becomes smaller. 
The predominant period of the ground un- 
rests is longer than about 0.1 sec. Because 
of these facts, we can get a better s/n by 
making the observation of microearthquakes 
with the velocity seismometers than displace- 
ment seismometers. Therefore, an electro- 
magnetic seismometer connected with the 
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galvanometer of proper period longer than 
0.1 sec records ground unrests rather than 
microearthquakes. The seismometer which 
the present writer designed according to the 


Fig. 1. Ultra Sensitive short period seismometer 


with the cover removed. The instrumental 


constants are given in Table 1. 
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already mentioned requirements is shown in 
Fig. 1. The instrumental constants are given 
in Table 1. The proper period of the pendu- 
lum is adjusted to be 0.1 sec. The impedance 
of the transducer is made rather low for 
direct connecton with the grids of the first 
vacuum tubes in amplifier, so that the trans- 
ducer and lead wire do not pick up electric 


Table 1. Constants of the seismometer shown 
ins Higws he: 

Mass of the pendulum 146 grams 
Moment of inertia BUNS Okey 
Equivalent length 6.57 cm 
Length between center of 

gravity and rotation axis 3.91 cm 
Mechanical magnification 1E25 


Sensitivity | {376 volts/em at 10 eps 
Natural period 3~10 cps 
Internal resistance | 5000 ohms 
Critical damping resistance| 60 x 103 ohms 

|  (To=10 cps) 


Va 


R:=100 ka Vib 
C0) fd R,=500 ko V.=6 oa 
C;=0.002 nfd R3=100 ka V3=6 SL7 
C,=20 pd Rie 5OkG Vi=12AU7 
alin ue R;=200 ka 
6= 20 pfd Re=100 ka = 
C;=20 pfd R;7=100 ko ais 
Re= 30a 
Rea 20 ka 


Pige2, 
The transducer is shown in Fig. 1. 


Circuit diagram of the low frequency 


amplifier for the Ultra Sensitive seismometer. 
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Fig. 3. Frequency response curve of the ampli- 
fier shown in Fig. 2. The sensitivity is 5 
micro-volts (peak to peak)/1mm, when a 
galvanometer with the sensitivity of 30 mm/ma 
is connected. 


10" kine/mm 


5x10" 


10 CPS 1OOCPS (OOO CPS 


Fig. 4. Overall frequency response curve of the 
Ultra Sensitive velocity seismometer consisting 
of the tranceducer shown in Fig. 1 and the 
amplifier shown in Fig. 2. 


Fig. 5. Highly sensitive electrodynamic seismometer having natural period 1 sec. 


disturbances. The coils are connected push- 
pull type for the same purpose. 

The circuit diagram of the amplifier is given 
in Fig. 2. Input and output transformers are 
not used. A galvanometer with the proper 
frequency of 300 cps is connected with the 
output terminals. The frequency response of 
the amplifier connected with the galvano- 
meter is shown in Fig. 3. The circuit noise 
generated in the vacuum tube and in other 
parts of the first stage is about 1 ~ volt»_» if it 
is reduced to the value at the input terminals. 
The overall frequency characteristic of the 
seismometer is shown in Fig. 4. Its maximum 
sensitivity is 1.2x 10° mm/kine (mm/cm/sec). 
Namely, the magnification of this seismo- 
meter can be raised up to 1.5x10" at a fre- 
frequency of 20 cps. 

The pendulum is designed to have no un- 
favourable peaks in its response curve in a 


frequency range between 10 cps and 200 cps. 
The noise due to the effect of Brownian 
motions due to molecules of the air surround- 
ing the pendulum is negligibly small in the 
present observation. The energy of vibration 
of the ground particles with a velocity ampli- 
tude as small as 10-7 cm/sec is yet sufficiently 
large for a passive transducer to convert it 
into the energy of electric current above the 
noise level of the amplifier. In the present 
observation, three seimometers were used. 
Two of them were of the type described 
above; one was set horizontally and the 
other, vertically. The third one was also a 
horizontal seismometer set up in the same 
direction as that of the 10 cps_ horizontal 
seismometer. This is an electromagnetic seis- 
mometer with a proper period of one second 
(Fig. 5), its output being amplified. The cir- 
cuit diagram is shown in Fig. 6 and the 
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Vs Ve 
(il pfd R,=500 ka Vi=6SJ7 
C,=0.01 pfd R2=300 ka V2=6SJ7 
C3=100 ufd R3=250 ka V3=6SJ7 
Ry= 50 koa V.=6J5 
V5=6J5 
B+ =200 Vs=6SJ7 


Fig. 6. Circuit diagram of the very low frequency amplifier for a seismometer having natural 


period 1 sec. 


quency reponse in Fig. 7. The overall 
frequency characteristic of this 1 second 
seismometer and amplifier is shown in Fig. 8. 
Its sensitivity is higher than 10° mm/kine. 
Operating these three seismometers, one 10 
cps vertical, one 10 cps horizontal and one 
lsec horizontal seismometers, legible records 
of near-by microearthquakes can be obtained 
with clearly recognisable commencements of 
both P and S phases. An example of a 
near-by microearthquake recorded is shown in 
Fig. 9. 

By making the temporary observation only 
for 24 hours by these seismometers, we shall 
be able to collect the records of many near-by 
earthquakes which are sufficiently numerous 
for the seismometrical research. However, 
since one of the writer’s purposes is to study 
the secular trend of microearthquake occur- 
rence in comparison with that of larger earth- 
quakes of M 4 or 5, the temporary observa- 
tion by the Ultra Senstitive seismometers 


was planned to be made for 2 or 4 hours 
every night, when disturbances or ground- 
unrests are small. The periods, during which 
the seismometers and oscillograph are operat- 
ed, and the number of the shocks registered 
in each period are given in Table 2. 
According to the results of the former 
studies on microearthquakes, a hundred or 
more near-by earthquakes are expected to be 
recorded every day at Mt. Tsukuba, and 
therefore a large quantity of the recording 
bromide papers will be needed. In the present 
observation at Mt. Tsukuba, multi-elements- 
oscillograph and photographic paper tape recor- 
der were used. The tape was 125mm wide, 
and the running speed was 20cm per minute. 


§ 3. Seismograms 


In case that the observation is made at a 
place not so quiet, a small earthquake can 
only be detected by comparing two or three 
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IcPs lo cps \Osec (sec tocps 100 cPS 
Fig. 7. Frequency response curve of the ampli- Fig. 8. Overall frequency response of the high 
fier shown in Fig. 6. The sensitivity is about sensitive velocity seismometer consisting of the 


seismometer having natural period 1 sec and 


2 micro-volts (peak to peak)/1 mm. 
the amplifier shown in Fig. 6. 


Fig. 9. Trace a: Seismogram by 1sec horizontal seismometer. 
Trace b: Seismogram by 10 cps horizontal seismometer. 
Trace c: Seismogram by 10cps vertical seismometer. 


seismometers installed at different positions at one point. Consequently the discrimination 
tens of a meter apart, since the records of between earthquake motions and disturbances 
the ground unrests at the positions tens of must be made by the difference in features 
meter apart have little correlation to each of the vibrations on the record. 

other, whereas the records of an earthquake The seismograms of near-by shocks of 
recorded at those positions are quite similar. hypocentral distance less than a hundred kilo- 
In the previous experiments’, the present meters are shown in Figs. 10 and Lie srigw, 
writer employed a tripartite station method to the seismograms of the earthquakes of hypo- 
discriminated small earthquakemotions from central distance about several hundred kilo- 
disturbances. But in the present experiment meters are given. In those figures, the trace 
such a tripartite station method could not be “a” is the seismograms recorded by the 1 
employed, and the seismometers were set up second seismometer connected with the ampli- 
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Table 2. List of time during which the seismo- 
meters were operated, and the number of 
microearthquakes registered. 


The time during which | 


Photographic the seismometers were | Number 

PvE operated of 
eis ees Sa ta ak 
i | | 
2 h m hm | 

3 | Aug. 1. 23 45~01 35 33 
+ 2. 02 05~04 10 30 
5 3. 00 05~01 43 © & 
6 3. 02 10~03 04 16 
7 4. 00 06~01 55 18 
8 5. 01 00~02 48 g 
) 5. 03 09~05 08 12 
10 5. 05 36~07 30 11 
11 6. 00 17~02 14 13 

12 { 6. 23 02~23 40, 6 | 

| 7. 18 55~20 01 4 J 
13 7. 20 30~22 20 8 
14 7. 22 36~00 29 8 
15 9. 20 45~21 46 5 
16 9. 22 08~23 58 6 
17 10. 23 00~00 50 16 
18 11. 19 40~21 34 10 
19 11. 23 40~01 36 | 17 
20 12. 13 15~15 16 51 
21 12. 22 50~00 44 13 
22 13. 20 55~22 51 12 
23 13. 23 40~01 38 18 
24 14. 15 45~17 45 14 
25 14. 20 55~22 51 | 1a 
260 15. 00 05~02 03 12 

fier given in Fig. 7; the trace “b”, the 


seismogram of the 10 cps seismometer of the 
horizontal component; and the trace “c”, the 
seismogram of the 10 cps vertical component 
seismometer connected with the amplifiers 
shown in Fig. 3. 

In examining the seismograms shown in 
Figs. 10, 11 and 12 the following points can 
be clearly noticed: 

1) In the seismograms ‘“‘b’’ and ‘‘c’”’ of 
Fig. 10, the waves of very short periods are 
predominant. On the other hand, the “P 
waves” of the trace “a” have small ampli- 
tude and are almost masked by the ground 
unrests. 


This phenomenon can be interpreted by the 
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lifference of the frequency responses between 
the one second seismometer, seismogram of 
which is shown in trace “a”, and 10 cps 
seismometers, seismogram of which is shown 
m trace “b” and “c” 

_ The fact described above shows that the 
Beriods of the P waves of very small near-by 
earthquarkes are short and they are recorded 
only by the seismometers with the response 
given in Fig. 4. The seismometer which is 
sensitive in a band between 1 cps and 30 cps 
is not appropriate for recording the P waves 
of such a small near-by earthquake. 

_ The oscillogram in the trace “a” can hardly 
be identified as a record of an earthquake. 

_ This fact leads to a conclusion that the 
predominant frequency of the P waves of 
microearthquakes is more than 50 or 60 cps. 
2) On the contrary, the seismograms of the 
rather distant earthquakes as shown in Fig. 
12, give larger trace amplitudes in the trace 
fa” than in the trace “b” and “c”. The 


Fig. 


o 


20 to 30 cps, and those of the S waves are 
about 13 cps to 30 cps. But the waves of 
ten odd cps oscillate only for a few times. 

In the seismograms of very near earth- 
quakes, the waves of short periods which 
have never been recorded in the usual routine 
observation are predominant, but in case of 
more distant earthquakes, the predominant 
periods are not very short. They are about 
9.1 second or longer. 

This phenomenon will be interpreted from 
two different points of view. One is that the 
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seismograms given in Fig. 12 show that the 
predominant periods of earthquake-motions of 
the shocks of hypocentral distance about 
several hundred kilometers or more are longer 
than 0.1 second. 

The seismometer which gives the trace 
“a” is sensitive in the range of longer periods 
than the seismometer shown in Fig. 1 and 
Fig. 2. The seisitivity of this seismometer 
falls about 6 db at a frequency of 40 cps and 
is much higher in shorter periods than the 
usual electrodynamic seismometers for the 
teleseismic recording. Nevertheless, the sen- 
sitivity in short periods is quite insufficient 
for recording very small near-by earthquakes. 

The seismogram in Fig. 13 is an example 
recorded on the photographic paper running 
faster. The hypocentral distance of this 
earthquake is about 60 km, and this is twice 
that of the shocks shown in Fig. 10. The 
seismogram in Fig. 13 shows that the pre- 
dominant frequency of the P waves is about 


Scm 


13. 


waves of shorter periods decay more rapidly 
than the waves of longer periods, and as the 
hypocentral distance becomes longer the 
waves of longer periods become more pre- 
dominant. The other one is; among the 
earthquakes registered by a certain seismo- 
meter, the shocks with smaller hypocentral 
distances have smaller magnitude than the 
more distant shocks with the trace amplitudes 
of the same order, and the predominant 
periods become longer with larger magnitudes, 
therefore the waves from more distant earth- 


94 Toshi ASADA 


quakes have the longer predominant periods. 
In the latter alternatives, the effect of the 
viscoelesticity of the earth’s crust is estimat- 
ed to be small. In the present writer’s 
opinion, the second reason seems to be more 
important. 


§4. Frequency distributions of the maxi- 
mum amplitude A, and the magnitude 


M. 


When an instrumental magnitude scale was 
defined by C. F. RicuTEerR™ in 1935, the 
quantity MM was naturally proposed for more 
detailed study in the statistical seismology. 
The principal object in the statistical seis- 
mology is the investigation of time and 
space distribution of earthquake occurrences. 
As the number of earthquakes occurring in a 
certain area depends wholly on their magni- 
tudes, the frequency distribution of magnitude 
is one of the most fundamental relations in 
the statistical seismology. The study con- 
cerning the number of earthquakes will be 
meaningless if the magnitude frequency rela- 
tion is neglected in course of that study. 

The problems on the magnitude-frequency 
relation have been studied by GutTenpere and 
RicHTER®, Kawasumrl®, Tsusor, IsHimoto and 
ipa,” Suzuki and Asapa?, The equation of 
the relation given by Gurensera and RIcHTER is, 


log N=a+b(8—M) , (4.1) 


where, the value of “a” depends on the 
period of the observation as well as the seismic 
activity of a certain closed region being studied, 
and the coefficient “5” has also a geophysical- 
ly important significance. “6b” is one of the 
quantities which characterise a seismic region, 
when “6” is determined from the data on 
shocks occurring there. 

In studying the relation, it is possible to 
adopt Ao, the amplitude at the hypocenter 
as a variable instead of M. Maruzawa2 
and Suzuki discussed the relation between 
frequency distribution of Ay and that of A, 
the amplitude at a station. They proved 
mathematically that if the functional form 
of the frequency distribution of Ay is Ay-™, 


the frequency distribution of A, f(A)dA, is 
expressed by an equation, /(A)dA= const. 
A-™dA, under the following assumption. The: 
assumption is that the decay of seismic. 
waves depends only on the properties of the | 
medium and not on the magnitude of the 
earthquake. However, the decay of waves 
actually depends on M, but the degree of 
dependence is negligible in studying near-by | 
earthquakes of magnitude up to about 5. 

By making use of this relation between. 
A, and A in studying the magnitude fre-_ 
quency relation, the records of a seismometer 
operated at a single station can provide 
sufficient data for studying the amplitude fre- 
quency relation of near-by earthquakes. As 
will be explained later, the equation, f(A)dA 
=const. A-"dA, is almost equivalent to the 
equation, log N=a+d(8—M). 

According to the statistical study by Isut- | 
moto and Itpa™, Suzuki, and Asapa,?” the 
equation concerning the frequency distribu- 
tion of A, 


NA)dA=const.A-"dA , (4.2) 


always holds good. The value of m is about 
1.7 to 1.9. Since M is defined as 


M=const.+log Ay , 


we can compare the equation (4.2) with the 
result obtained by GuTeNserRG and RIcHTER. | 
Transforming A in the equation, (4.2) into 
M, we have, 


NM )dM=const. 10-@-D¥dM 


Since m is from 1.7 to 1.9, this may be 
written as follows, 

N(M)dM =const. x 10-°7™™ dM , 
or, 


NM )dM=const. x 10-°9* dM . 


This is the same as the equation (4.1), be- 
cause the values of “5” determined by 
GUTENBERG and RicHTER, Tsusor, and other 
authors are from 0.7 to 1.0. 

(4.1) and (4.2) are only empirical equations 
and were not derived theoretically. The fre- 
quency distribution of magnitude could be 
expressed by some other equations (Ko- 
muRA)*)?), Moreover, the equations (4.1) and 
(4.2) are not proper as a function for pro- 


ability density, because N becomes infinitely 
large as Ay approaches zero. According to 
he observations made thus far, the number 
f small near-by earthquakes increases as ex- 
ressed in (4.2) when the sensitivity of seismo- 
meters is made higher in a range from 10 
cps to 100 cps, and appears to keep on in- 
creasing indefinitely if the observation could 
be made with higher sensitivity than the 
present level. Little information has been 
accumulated as to the lower limit of the range 
of A in which the equation (4.2) is valid. 

In carrying out theoretical research con- 
cerning the statistical quantities in seismo- 
logy, it is convenient to adopt the relation, 
(4.2) as a magnitude-frequency relation. 

The upper limit in the equations, when the 
statistics are made on the shocks occurring 
in a certain closed part of the crust, has not 
yet been studied either. This problem will 
likely germinate a clue to forecasting the 
eventual occurrence of major earthquakes in 
the region. 

In studying these problems, it is desirable 
to determine the magnitude of every micro- 
earthquake. But it is technically difficult to 
carry out the routine observation by the 
Ultra Sensitive seismometers at several points 
simultaneously. Consequently, the accurate 
determination of the magnitude of individual 
earthquake is also difficult. 

As the frequency distribution of A is almost 
equivalent to the frequency distribution of 
magnitude, our purpose to investigate in 
detail the magnitude-frequency relation in the 
microearthquake occurrences will be accom- 
plished through studying the frequency dis- 
tribution of the maximum trace amplitudes 
instead. 

In the present observation, the trace am- 
plitude less than several mm are masked by 
the ground-noises, and on the other hand, the 
trace amplitude larger than sixty or seventy 
mm is too large for registration. The statis- 
tics are made on the trace amplitudes, A, 
between 5 mm and 60 mm in the present 
study. The frequency distribution of A is 
obtained from the data on the maximum 
trace amplitudes of the seismograms recorded 
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by the short period horizontal seismometer 
with the frequency response given in laiez, 
The result is shown in Table 3 and in Fig. 14. 

The graph given in Fig. 14 shows that this 
frequency distribution is well expressed by 
the equation, MA)dA=kA-"dA. The theo- 
retical values for the frequency calculated ac- 
cording to the equation, (4.2), are listed in 
Table 4 compared with the observed values. 
The difference between two sets of values is 
not statistically significant. The velocity 


Table 3. Frequency distribution of trace 
amplitudes 

Trace cies Frequency 
0.5~1.0 | 60 
1.1~1.5 46 
1.6~2.0 51 
2.1~2.5 29 
2.6~3.0 26 
3.1~3.5 9 
3.6~4.0 14 
4.1~5.5 7 
4.6~5.0 12 
5.1~5.5 3 
5.6~6.0 6 
6.1~6.5 4 
6.6~7.0 2 
(eee) 4 
7.6~8.0 2 
8.1~8.5 1 


Trace Amplitude > 


lem Dee SrA MS 10 
Ground Amplitude > 
1.7 x 10-5 cm/sec 
Fig. 14° Frequency distribution of velocity 
amplitudes. 


1.7x10-4cms 
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magnification of the ground-motions is about 
8.6x10-® kines per 10 mm of the trace am- 
plitude. In other expression, the amplitude 
magnification of the ground particle is 6.8 
10-7 cm per 10 mm of the trace amplitude at 
20 cycles per second. 

The results given in Tables 3 and 4 demon- 
strate that the equations (4.1) and (4.2) hold 
good for the magnitude frequency relation in 
the occurrence of very small earthquakes as 
well as of larger earthquakes. The range of 


Table 4. Calculated values of the frequency 
distribution of trace amplitudes. 


‘ Frequency distribution 
Trace amplitude 

ee) \Observed values ee 
1.1~2.0 97 111.48 
ESS ee 55 43.46 
3.1~4.0 23 23.08 
4.1~5.0 19 14.39 
5.1~6.0 9 9.81 
6.1~7.0 6 Took 
(el 18.0 6 yea | 


| \ 


A in Tables 3 and 4 is from 5 mm to 60mm, 
or from 4.3x10-* cm/sec to 5.2x10-5 cm/sec. 
The above mentioned results show the vali- 
dity of the equation in this range only. 

It is also an important problem whether 
these are valid also within a wider range 
from micro-earthquakes of M less than 0 to 
earthquakes of M 4 or 5. If the frequency 
distribution of the amplitude of earthquakes 
of magnitude from less than zero to 4 to 5 
is expressed by the equation, (4.2), we shall 
be able to estimate the number of shocks 
of magnitude 4 or 5 occurring in a certain 
region and during a certain period of time, 
through the observation of micro-earthquakes 
during much shorter period of time in the 
same region, 

Denoting the number of shocks recorded by 
a seismometer with the magnification v; by 
N;, and denoting the number of shocks re- 
corded in the same period of time by other 
seismometer with the magnification v, at the 
Same position by Ns, the following relation 


is obtained as will be proved in the late 


( =) m-1 eel 
Vy Nz : 


where m is the exponent in the equation) 
(4.2). | 

The information on the relation between 
the number of near-by shocks recorded by a : 
Ultra Sensitive seismometer and the number 
of those registered by the seismometer with 
lower sensitivity operated at the same posi+ 
tion is indispensable for our purpose. 

MiyamurA” and his colleagues carried out 
preliminary experiments of tele-recording wit 
ultra short radio waves. Their experiment 
were to have the earthquake-motions at Mt. 
Tsukuba recorded at the Earthquake Researc 
Institute in Tokyo. The frequency respons 
of the seismometer employed by them was 
almost the same as the one given in Fig. 8. 
Its sensitivity was 7x10? mm/kine or th 
magnification was 90,000 at 20 cps. Thi 
value is about 1/160 of the maximum sensiti- 
vity of the seismometer shown in Fig. 4. 

It is not appropriate, however, to adop 
this value, 1/160, as vi/v.. The seismograms. 
recorded by the Ultra Sensitive seismometer: 
can be identified as earthquake-motions only? 
when their double trace amplitudes are larger: 
than 5 mm, for the seismograms of smaller’ 
amplitudes are usually masked by the dis-- 
turbances. However, the seismograms record-- 
ed by the low sensitive seismometer having: 
the magnification of 90,000 at 20 cps are not: 
disturbed by ground-noise, and those with 
trace amplitude larger than about 1 mm can: 
be distinguished from the noise. Consequent- 
ly, the maximum amplitude of the minimum: 
earthquake-motions which can be _ identified’ 
clearly by the Ultra Sensitive seismometer is. 
not 1/160, but about 1/30 of that of the 
minimum earthquake-motions recorded legib- 
ly by the seismometer used by Mryamura 
and his colleagues. 

About 200 earthquakes were registered 
during the period of every 24 hours by our 
Ultra Sensitive seismometers placed at Mt. 
Tsukuba. The seismometer with the sensi- 
tivity of 7x10? mm/kine, recorded about 10 


(4.3) 


arthquakes per 24 hours, 
Putting these values into the equation, (4.3), 
jiwe get m=1.9, as 


( 30 1 re 200 
1 10 


This value of m is conformable with the 
iresults of studies to date.2” 

|| Further, the relation between the number 
\jof small near-by earthquakes registered only 
by the Ultra Sensitive seismometer and the 
number of near-by earthquakes felt by people 
\jat Mt. Tsukuba should be discussed. The 
\jseismograms recorded by the accelerometers 
jjat Mt. Tsukuba show that the earthquake- 
imotions of 12.5 cps are predominant.2” The 
i lower limit of earthquake-accelerations sensible 
i by human feeling is 1.7 gals at Mt. Tsukuba. 
| Accordingly, this limit corresponding to 2.2 
10- kines, if expressed in terms of velocity. 
Now we put v,=5x8.7x10-7 kines, and v.= 
2.2x10-2 kines, and since m is 1.9, we get 
1/3300 as the value of N,/N2, according to 
the equaion, (4.3). 

This result shows that one sensible earth- 
quake will occur satistically every 16 days, 
the highly sensitive seismometer records 200 
near-by earthquakes every 24 hours. In other 
words, about 20 to 30 near-by earthquakes 
will be felt by people at Mt. Tsukuba every 
year. The sensible earthquakes observed at 
Mt. Tsukuba are listed in Table 5. Forty 
five earthquakes were left there in 1955. 
The number of sensible earthquakes calculat- 
ed by the equation (4.3) agrees with the 
actual value. The list of earthquakes of 
magnitudes above 4.5 occurring in the vici- 
nity of Kanto district in 1955 is given in 
Table 4. Fig. 15A shows the distribution of 
their epicenters. The geographical distribu- 
tion of the epicenters of earthquakes occurring 
during the period from 1914 to 1940 is also 
given in Fig. 15B. 

The above described results lead to the 
following conclusions. The frequency distri- 
bution of A may be expressed by the equation 


NA)dA=kA-“dA , 


and the amplitude range where the equation 
holds good is from 2%x10-° kines to the am- 
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Fig. 15A. Epicenters of earthquakes of magnitude 
above 4.5 occurring in the vicinity of Kanto 
district during 1955. 


plitude of near-by earthquake of M about 5. 
If A denotes the displacement amplitude of 
near-by earthquakes, the equation is valid 
within a range of from 10-’cm to 107? or 
Om actor 

Consequently the equation, 


log N=a+b(8—M) , 


holds good also in a wide range of M. This 
magnitude range from —1 to 5 will be dis- 
cussed in detail in the later chapter. 

The Wood-Anderson seismometer is not 
operated at the observatory of Mt. Tsukuba, 
but now it is possible for us to estimate how 
many near-by earthquakes would have been 
registered by the seismometer of this type 
during a certain span of time, if it had been 
operated there. As the magnification of this 
seismometer is 2,800, and the predominant 
frequency for near-by earthquakes of M 3 or 
4 is from 10 cps to 15 cps here, the velocity 
amplitude of the earthquake with the maxi- 
mum trace amplitude of 2 mm is about 6.8 
10-3 kines. According to the equation (4.3), 
putting 1.=5x8.6x 10-7, v2=6.8x 10"°, m=1.9, 
and N,=200, we have the result that one or 
two near-by earthquakes would be registered 
every day by the Wood-Anderson seismometer, 


98 Toshi ASADA 


Table 5. Date of occurrence and intensity of earthquakes felt at Mt. Tsukuba 1955 
(intensities are expressed in Japanese C.M.O. scale). 


OMNA TO FwWNY 


ie) 
Ea 
=F 


a | | | | | 


if it were operated at Mt. Tsukuba. 
Summarizing the above statements, we can 
reasonably say that the number of micro- 
earthquakes can easily be estimated if the 
number of sensible earthquakes during a 
certain span of time is known, because the 
equations (4.1) and (4.2) hold good not only 
in a magnitude range between 4 and 6, but 
also in a wider range between —1 and 6. On 


| | | | 


the contrary, if the observation of micro- 
earthquakes by the Ultra Sensitive short — 
period seismometer is made during a period — 
of tens of an hour in total, the number of the 
sensible near-by earthquakes can be estimat- 
ed by putting the number of microearthquakes 
and the value of sensitivity of the seismo- 
meter into the equation (4.3). The value of 
m in the equation N(A)dA=const. A-™qd Aig 
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Table 6. List of earthquakes of magnitude above 4.5 occurring in the vicinity 
of Kanto district in 1955. 


Time of occurrence (J.S.T.) Hypocenter 
Date 4 ow = | Magnitude 
h m s N E (H) km 
Jan ip eae 16 35 46 36.2 139.9 50 D 
12 16 18 32 37.3 141.5 | 90 fe 
16 04 5 | 15 35.2 141.5 20 
17 11 7 ed 49 | 35.5 140.2 80 @ 
19 13 58 | ay | 36.4 140.7 50 D 
a | oO | 32 32 | 36.6 141.0 | 30 D 
Feb. oh 04 | 56 46 | 34.3 141.0 | co 
4 «| 03 01 53 36.5 141.0 | 40 | 
March 2 07 ean | meets 35.5 138.9 | 20 D 
g 07 DO 27 35.5 138.90) 4 10 D 
4 18 24 15 35.3 139.5 |  80~100 | 
16 01 09 29 36.4 10.7 50 D 
18 17 59 21 35.4 140.9 0~10 E 
20 02 16 15 36.2 139.8 50 D 
29 13 12 | 23.2) 35.7 140.4 60~70. | E 
April 5 04 | 13 07 Cd 35.6 141.7 30~40 
9 to 33 59 7.2 141.7 40 
is} 02 | 10 50 35.5 141.1 40 D 
15 | 13 52 46 36.3 141.2 40 D 
Ae iA eet 926 39 36.4 140.8 40~60 
Ge | 09 03 ot | 36.6 140.9 90 
May 12 00 48 52 Ci Seat) - 141.6 40 C 
12 17 36 | OT. | 35.9 | 140.9 20 E 
13 10 59 11 35.7 140.3 | 40 E 
21 Ta 53 54 oye! 141 oO 40 E 
26 15 | 48 06 37.1 141.9 40 D 
26 21 45 43 36.2 141.3 40 
27 05 18 25 36.2 139.8 50 
June 15 01 34 57 36.3 142.1 50 D 
15 02 21 57 36.3 142.1 40 € 
16 01 43 08 37.4 141.8 40 
23 17 39 AT | 37.4 141.4 60 C 
July 9 19 | 04 ve 35.8 140.1 80 E 
14 19 17 5 | 36.4 141.2 40 D 
17 17 03 12 35.7 140.1 80 E 
24 20 02 11 35.7 140.7 7 «| G 
Aug. eer 13 09 27 36.4 141.1 30 D 
Sept. os 10 D4 | 48 35.6 140.3 70 E 
3 17 50 30 35.6 140.4 70~80 E 
13 15 23° | 20 35.9 140.3 80 E 
20 | 07 =| 33 | 04 35.7 140.8 40 D-i 
Oct. Coe ae i 46 06" (ureka5i6or| ©) 140:5 40 D 
ice, 1s 56 09 | 35.6 | 140.9 10d E 
19 11 04 23 B54 i 140.3 40~50 E 
20 | 22 10 40 | epige me 14057} 50 D 
| 16 00 00 36.5 140.7 | 40~60 E 
oe | 12 eid 08 | 36.7 141.0 70 2 
Nov. 3 ar 304 58 11 36.3 141.8 40~60 D 
Dec. 2 07 57 26 35.6 140.5 70~80 
8 01 05 iW 36.5 140.7 50 
03 Al 07 35.75 141.25 40 D 
S 18 29 01 36.6 140.9 80 E 
15 17 34 29 36.1 139.9 50 2 
1Susor 13 24 4 01 35.7 142.2 40 ues 
18 14 | a3 | 10 | 35.2 139.4 | 50 ne 


100 


also necessary. In Mt. Tsukuba region, m 
was estimated to be 1.9 in the present in- 
vestigation. According to our studies to date, 
the value of m is from 1.7 to 1.9 in any 
region in Japan. 

The time series of earthquake-occurrences 
is always assumed to be stationary in the 
course of the calculation mentioned above. 
This assumption is, of course, not strictly 
correct, but it may be safely assumed; when 
the accuracy of calculation is low, and the 
period of observation is sufficiently long, for 
instance, more than one year for sensible 
earthquakes, and more than two or three 
weeks for microearthqaukes. 


The constant A in the equation (4.1), 
denotes naturally the displacement amplitude, 
but in this chapter A denotes the velocity 
amplitude which has the dimensions of kine 


or cm/sec. The equation should be written, 
N(Av)dAv=const. Ay~™ dAo , (4.4) 
instead of 


N(A)dA =const. A-™dA , 


where A» is the maximum velocity amplitude 
of an earthquake. 

The results in Table 3 and Fig. 14 confirm 
the validity of this equation. The results of 
other studies!?!® with velocity-seismometers 
have supported also the validity of the equa- 
tion (4.4). In other words, the frequency 
distribution of the maximum displacement 
amplitude and that of the maximum velocity 
amplitude are expressed by the equations in 
the same functional form. 

This fact may be explained as follows.1 
In seismic waves dealt with in the present 
study, the wave train which gives the max- 
mum displacement amplitude is assumed to 
be also that which gives the maximum in the 
velocity amplitude. The assumption seems 
adequate from our study in which the ac- 
celeration seismograms were compared with 
the displacement seismograms of the same 
earthquakes. 

Now we have, 


Ay=(22/T)A , (4.5) 
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where TJ is the period of the wave train of 
maximum amplitude both in displacement and 
in velocity. We may call “J” the pre- 
dominant period of the earthquakes-motions. 
Another relation between A, and A, 
Av=C-A?, (4.6) 
is assumed, where C is a quantity independent. 
of Aand A», and p is a constant. Then, the: 
equation, | 
Avy=CoAo? , (4.7)) 


holds good also, where the quantities with) 
suffix o are those at hypocenter. Cy is ai 
constant. 
As already referred to, the frequency dis-- 
tribution of A, when a number of near-by’ 
earthquakes are registered at a certain re-- 
ceiving position, may be expressed by the: 
equation, 
N(A)dA=const. A-™ dA. (4.8) 


On the other hand, the studies on the magni- 
tude frequency relation showed that the 
equation, 

IN(Ap)dAp=const. Ay ™ dAy 5 (4.9) 


is adequate. As we have the relation, A»,= 
CyAy”, we can transform the variable Ay in 
(4.8) into Av). 
We have, 
N(Ap)dAp=k- Co“ V? Avg ™!?d Ay ; 
MAv)dAry=h: Cot? Ag 20 dA 


vO 
Care 
pay ened LE eg 8 AS. : 


p 


| 


or 
N Av)dAvn=KAv,™ dAv, , 
pm 
= —— | il 5 
p (4.10) 


jen lake 
p 


m 


=const. 


If the equation (4.10) holds good concerning | 
earthquakes occurring in certain region, we > 
shall have the relation, 


N(Av)dAv=const. A,-™’dAy, (4.11) 


as already verified mathematically and em-_ 
pirically too. This is the frequency distribu- 
tion of their maximum velocity amplitudes at 
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Fig. 15B. Epicenters of earthquakes occurring in the vicinity of Kanto district during a period 
from 1914 to 1940. The relation between M; and M is as follows: 
M=4.85+0.5 Mz . 


A: Tsukuba B: Hongo, Tokyo 


a receiving position. mometers at Mt. Tsukuba shows m’ being 
According to the results of the observation 1.88 as already mentioned. 

by velocity seismometers, the value of m’ is As we have, m=1.74 and m’=1.88, we get 
usually slightly larger than m, the exponent p=0.84, putting the values of m and m’ into 
in the equation, (4.8), expressing the frequency (4.10). From (4.6) and (4.7), we have, 
distribution of the maximum displacement am- T=CA 

plitudes. The value of m obtained experi- } marin 

mentally by Isuimoto and Ima is 1.74. The 
result of the observation by the velocity seis- Pr=conste An’ 


P being 0.84, we get 
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The accuracy in determining the value of 
pis low. The second decimal is of little sig- 
nificance. 

The facts that the frequency distribution of 
A,» has the same form as that of A, and m’ 
is slightly larger than m can be explained 
clearly by assuming the predominant period 
T being dependent on Ap or M. 


The list of earthquakes of magnitude larger 
than about 4 in and near Kanto district in 
1955 is given in Table 6. Their epicenters 
are plotted in the map in Fig. 15A and 15B. 
Every year tens of earthquakes are felt by 
people in Mt. Tsukuba region, but great earth- 
quakes have never occurred there. Accord- 
ingly no aftershocks have been observed in 
the region. Thus the microearthquakes which 
amounted to 200 in number in 24 hours by 
the observation of the Ultra Sensitive short 
period seismometers are definitely not after- 
shocks of great earthquakes. 

As the seismicity of this region is shown 
in Table 5 and 6 and Fig. 15A and 15B, the 
seismicity of any other region can be com- 
pared with that of this region. By means 
of the equation (4.3), it will be possible to 
estimate how many microearthquakes will be 
recorded by an Ultra Sensitive seismometer at 
a place. 


§5. Energy of microearthquakes 


The frequency distribution of very small 
near-by earthquakes were described in the 
previous chapters. The estimation of the 
energy released in seismic waves in such 
small earthquakes will be discussed in this 
chapter. 

The hypocentral distances of the greater 
part of very small earthquakes which have 
been recorded at Mt. Tsukuba are approxi- 
mately between 20km and 100km. It is self- 
evident that the energy of the earthquake 
that has the shortest hypocentral distance is 
the smallest among those of earthquakes 
which gave maximum trace amplitudes of the 
comparable order. 

An example of the seismograms of some of 
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the smallest earthquakes is given in Fig. 10. 


This particular earthquake has shown the 
maximum trace amplitude of 8mm and the 
hypocentral distance of approximately 30 km. 

In estimating the energy of seismic waves 
of such a small earthquake, it is not appro- 
priate to employ the magnitude scale system 
of GuTenserc and Ricuter. The relation be- 
tween M and E in that system does not 
seem to be reliable when the values are small. 
Other methods should be adopted for this 
purpose. We have two methods for roughly 


estimating the seismic wave energy of a small 


near-by earthquake. 
of E by estimating the kinetic energy of the 


One is to get the value | 


spherical waves radiated from the hypocenter. | 


The other is to obtain the value of E of an 
earthquake by comparing it with the energy 
of the waves from distant explosions. In 


both methods, the uniformity of the energy © 


radiated in all directions is assumed. 
In 1954, a blast of 2.7 tons of dynamite 


was made at Nozori,*® Gunma Prefecture, and 


the waves from this explosion were observed 


at Mt. Tsukuba by a vertical component seis- _ 


mometer of proper frequency 3 cps, connected 
with the amplifier having the same response 
as shown in Fig. 6. The maximum trace 
amplitude in S phase was 50mm. The sen- 
sitivity of the seismometer (pick-up) was 4.5 
volts/em/sec or 300 volts/em at 10 cps, and 
the gain of the amplifier and the galvanome- 
ter was adjusted to be 5.04 volts,-,»/mm. 
The overall sensitivity of this seismometer is 
9x10° mm/kine. The approximate energy of 
the above mentioned microearthquake can be 
estimated by comparing these data of the 
Nozori blast with the seismogram given in 
Fig. 10. 

The total energy of the explosion of 2.7 
tons of dynamite is about 2.7 x 10‘ calories, or 
11.3x10" ergs. We roughly estimate the 
seismic energy radiated from it to be about 
one thousandth of the total energy of the 
blast. Thus the energy of this artificial earth- 
quake propagated in the form of seismic 
waves is about 10! ergs. The waves of 
this artificial earthquake gave the maximum 
velocity amplitude of 2.8x10-5 kines at the 
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epicentral distance of 150 km. 

The microearthquake given in Fig. 10 had 
the maximum amplitudes of 0.35x10-° kines, 
and its hypocentral distance was about 30 km. 

The ratio of the maximum amplitude of 
the earthquake in Fig. 10 to that of the arti- 
ficial earthquake is 1:8. E of the latter should 
have been 1.610" ergs, if the waves from 
this explosion had given the amplitude 1/8x 
of 2.8x10-° kines at Mt. Tsukuba, which is 
equal to that of the natural earthquake. 

The ratio in the hypocentral distances of 
two shocks is about 1:5. According to the 
study by N. Den,!©® the maximum velocity 
amplitude of the waves from distant explosion 
decays inversely proportional to 1.2 to 1.5th 
power of the epicentral distance. 

The artificial earthquake of 10 erg would 
have given a velocity amplitude of (2.8 x 10-5) 
x51-5 kines at Mt. Tsukuba by using a value 
1.5 if the blast had been made at a distance 
of 30km which is equal to the hypocentral 
‘distance of the natural earthquake. Compar- 
ing the difference in the maximum amplitudes 
in two earthquakes, we get 1.4x10!! ergs as 
the value of & of the natural earthquake. 

We have another simple method for estim- 
ating the energy of an earthquake. The value 
of E can be calculated if it is identified with 
the kinetic energy of the spherical waves 
passing though a sphere witha radius, 7. An 
approximate value of E will be given by this 
method. 

The kinetic energy of the spherical seismic 
waves is expressed by the formula, 


E=2zr? vto(Ap)? , 


where v is the velocity of the waves, ¢ is the 
time of the duration of vibrations, o is the 
density of the medium, and A, is the veloci- 
ty amplitude of the waves. 
Putting into the formula, 

r=30 km, 

v=3.0 km/sec, 

#=1.0 sec, 

O=ell 

Ay=0.35 x 10-5 kines, 


we have E=1.4x108 ergs. However this 
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value of F is not equal to the energy of this 
earthquake, because the seismic waves decay 
on the way of propagation by several reasons. 
No comprehensive studies have likely been 
made concerning the decay of the seismic 
waves propagated in the crust. 

The values of the decay coefficient obtained 
by various authors are those of waves with 
longer periods, and therefore these are not 
applicable here. The decay coefficient of the 
waves with periods tens of a cycle per second 
is needed for the present purpose. When 
the waves are assumed to decay according to 
1/r-e-’” where + is the hypocentral distance, 
the value of h is 0.02km~™! according to 
Kawasuml,”? and 0.005km-! according to 
Wapatr and Hrrono.* The former value of 
h is applicable to the waves of frequency 
about 3 cps, and the latter is obtained from 
the data on wave of period about several 
seconds. They used the data on the waves 
of 3 cps and those of periods of several sec- 
onds respectively in the course of their studies. 
As the effect of the term 1/r is already in- 
volved in the formula representing the kinetic 
energy of spherical waves, we must have the 
energy of the spherical waves multiplied by 
e”” in order to obtain the value of the energy 
emitted from the hypocenter. 

Now, putting h=0.02km-! and r=30km, 
we have e~”""=0.55. In this case the waves 
lose only 20 percent of their energy in the 
propagation path 30km long. If it is assumed 
that the decay is due to the internal friction 
(Sezawa),*» it increases proportional to the 
square of the frequency of the waves. The 
predominant frequeney of the S waves of this 
microearthquake is above 30 cps, and the 
value of h for the waves of this phase will 
be above 2km-! according to the result by 
Kawasumi. Since y is 30km, we get, e’”"= 
This value of e~”” is too small. 

The hypothesis that the decay of seimic 
waves is due to such an internal friction as 
studied mathematically by Sezawa has not 
been supported by experimental results up to 
the present. Axi! found that the spectrum 
of the waves of near-by earthquakes depends 
on magnitude more than on hypocentral dis- 
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tance; in other words, the decay of the waves 
due to internal friction in the medium, if it 
exists, is not so large as to have any appre- 
ciable effect upon their spectrum during their 
propagation within several hundreds of km. 
According to the statistical study by Asapa” 
concerning near-by earthquakes of M_ less 
than 6, the dependence of the decay of the 
waves upon magnitudes is not very significant. 
Consequently, such a degree of dependence 
of the decay upon periods of the waves is 
also small. Although the results of these 
studies do not wholly disapprove the effect of 
the internal friction in the crust, they showed 
that the main part of the decay of seismic 
waves must be due to other causes than in- 
ternal friction. It seems inappropriate to 
assume the absorption spectrum of the crust 
to be proportional to the square of the fre- 
quency of seismic waves. Consequently it is 
not reasonable to estimate the energy loss of 
microearthquake waves during its propagation 
from the data of larger earthquakes of longer 
predominant periods. 

Assuming that the amplitude of wave decays 
proportionally to the 1.5th power of the hypo- 
central distance, we can show the loss of the 
wave energy in travelling 30km to be 70 
percent of its initial value. Accordingly the 
energy of this earthquake is evaluated at 
about 5x 108 ergs. 

We obtained two considerally different val- 
ues of EF of this earthquake according to the 
two different methods. This discrepancy is 
due to the lack of our knowledge on the 
absorption of the wave energy in the shallow 
part of the crust. 10" ergs will be a probable 
value for the E of this earthquake. 


We can have the approximate value of its 
magnitude also, if we assume that the decay 
of the waves is proportional to 1.5 th power 
of the hypocentral distance. As the maximum 
displacement amplitude of the waves of fre- 
quency 30 cps, or 0.35 x 10° kines in velocity, 
is 2x10-8cm, the maximum trace amplitude 
of the waves should be 0.084 microns, if they 
are recorded by a ‘‘ standard torsion seismo- 
meter ’’ of magnification 2,800 at a hypocentral 
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distance 100km. Accordingly the value of 
its M is —1.1. 

GUTENBERG and RIcHTER have given data 
regarding a shock of magnitude zero in their 
second paper on magnitude. According to) 
them, log B is —2 at a hypocentral peer | 
30 km, here B denoting the maximum trace: 
amplitude by Wood-Anderson seismometer, 
and so the maximum ground amplitude of the 
shock of magnitude zero is about 3.6 x 10-7cm. | 

The discussion stated above results in a 
conclusion that most earthquakes registered — 
at Mt. Tsukuba by the Ultra Sensitive seis- 
mometers have E, from 10% to about 10%. 
ergs. Consequently the range of applicability 
of the frequency distribution of the magnitude 
or Ay, described in chapter 4 are quantitative- 
ly determined. The equations (4.1) and (4.2) 
hold in a magnitude range between —1 and 
5. Transforming the variable M in (4.1) into | 
E, we have, 


NE)dE=const. E-/-1dF , 


using E=const. x 10°. 

This equation concerning the frequency dis- 
tribution of F holds in an energy range of 
from 10% to 1074 ergs. The equations (4.1), 
(4.2) and (5.1) show that in the part of the 
crust and the upper part of the mantle where 
the shocks of M 4 or 5 occur about three 
times a month, more than two hundreds micro- 
earthquakes of magnitude zero or less are 
occurring in every 24 hours. In an area where 
one or two shocks are felt in every year, an 
Ultra Sensitive short period seismometer will 
record 3,000 or 4,000 micro-earthquakes dur- 
ing the same period. 

The result of the experiment made by the 
present writer is that the equations 4.1 and 
4.2 always represent the magnitude frequency 
relation in shocks of magnitude between —1 
and 5, when they are occurring in a certain 
region. This result leads to the conclusion 
that we can always estimate how many micro- 
earthquakes will be registered by a highly 
sensitive seismometer at a certain place, if 
we know the number of larger shocks record- 
ed at the same place during a certain period. 
We can also calculate the rate of the occur- 


(5.1) 
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rence of shocks of M 4 or 5 ina certain seis- 
mic region, if we know the number of micro- 
earthquakes occurring there during a short 
period, say one week. 

But it must be noticed that the magnitude 
of the largest earthquake which will occur in 
a certain seismic region depends closely on 
the dimensions of the region. For instance, 
an earthquake swarm occurring in a very 
small region involves no shock larger than a 
certain limit which is related to the dimen- 
sions of that region. Consequently it is mean- 
ingless to calculate merely mathematically 
the rate of occurrence of shocks larger than 
a certain limit. 

However, it can be said without such a re- 
striction that no earthquake of magnitude 5 
or larger will occur more frequently than once 
per several year (4~7 years) in a part of the 
crust in which no microearthquakes of M 0 
or —l occur ina period of any 24 hours. 
- This inference is supported by an actual ex- 
ample. The present writer found that no 
microearthquakes had occurred in the part of 
the crust shallower than 20 km in the vicini- 
ty of Mt. Tsukuba, and in this part the earth- 
quakes larger than 4 in M have not occurred 
since 1924. This fact suggests that it may 
be possible to make a progress in the statis- 
tical study for forcasting the occurrence of 
large earthquakes by observing microearth- 
quakes. 

Historical studies on the destructive earth- 
quakes in and near Japan have shown that they 
occur cyclically at a certain specified place. 
In other words, even in Japan where pepole 
‘have experienced many destructive earth- 
quakes, the areas can be found, which have 
never been attacked by destructive earth- 
quakes since the dawn of the history. For 
instance, in the vicinity of Mt. Tsukuba 
destructive earthquakes have never occurred, 
although tens of earthquakes are felt there 
by pepole every year. Since destructive earth- 
quakes occur always in the shallow part of 
the crust, this fact seems to have a physical- 
ly close relation to the fact already described 
that the earthquakes seated shallower than 
about 20km have never occurred there. 
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The destructive earthquakes occurring on 
the Japan Sea side of Honshu are smaller 
in their magnitudes than those on the Pacific. 
The energies of the former earthquakes are 
considered to be stored in a part of the crust 
with an area of approximately 20 x20 km? and 
10 to 20km in depth. 

It must be emphasized that it is one of the 
hopeful methods of foretelling the possibility 
of the occurrence of destructive earthquakes 
in a certain area to operate the Ultra Sensitive 
seismometer during one or two weeks and 
determine whether microearthquakes are oc- 
curring there or not. 


§6. Geographical distribution of the foci 
of microearthquakes 


It is an important problem whether the seis- 
micity of microearthquakes shows the same 
geographical distribution as the seismicity of 
earthquakes of WM above 4. As explained in 
the previous chapter, the magnitude fre- 
quency relation of the earthquakes of M from 
1 to 5 can be expressed by the equation, 
N(M)dM=const. 10-’""dM. This fact suggests 
that the microearthquakes occur in the same 
area as the earthquakes of M 4 to5. But this 
is not a sufficient reason for the final con- 
clusion that the geographical distribution of 
microearthquakes is similar to that of larger 
ones. 

It is desirable to locate the hypocenter of 
every microearthquake by the observation at 
several stations. However, it is actually im- 
possible, because the stations equipped with 
the Ultra Sensitive seismometers must be 
located less than 20km apart from one an- 
other for this purpose. Moreover, some of 
the stations must be established not only on 
the surface but also under the surface of the 
crust. 

It is a matter of course that the distribu- 
tion of the foci determined by the observa- 
tion with the network of the stations is dis- 
torted when the number of stations is insuf- 
ficient (SAVARENSKy)*”. 

We have a reasonable expedient for this 
purpose. It is to compare the P-—S frequency 
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distribution of microearthquakes with that of 
larger near-by earthquakes observed at the 
same position. This is the only practical 
method, because the whole information ob- 
tained at one station can efficiently be used 
for investigating the distribution of the foci. 

P-S interval is approximately proportional 
to the hypocentral distance, as is already 
well known. However, it is a misinterpreta- 
tion? to consider the frequency distribution 
of P-S intervals as the space distribution of 
foci themselves. 


In the frequency diagram of P-S intervals, 
the total number of the observed earthquakes 
is expressed as an area enclosed by the fre- 
quency curve and the abscissa axis. The 
maximum amplitude, A, of the respective 
earthquake being left out of consideration in 
counting the number of earthquakes, each 
one occupies an equal space in this area. But 
putting an earthquake of larger amplitude to 
a lower position in each P-S class interval in 
the frequency diagram so that the amplitude 
of the earthquake plotted at the lowest posi- 
tion is largest among those of the earthquakes 
of the same P-S interval, the frequency curve 
of P-S intervals may be regarded as repre- 
senting the boundary between the observable 
earthquakes and the ones too small to be ob- 
served. The earthquakes, amplitude of which 
are small and masked by disturbances or be- 
yond the sensitivity of the seismometer, are 
seated outside the area enclosed by the fre- 
quency curve and the abscissa axis of the fre- 
quency distribution diagram of P-S intervals. 

Investigating the diagram from this point 
of view, it can be easily seen that the curve 
representing the frequency distribution of P-S 
intervals or hypocentral distances, D, depends 
not only on the geographical distribution of 
foci but also on the magnitude frequency re- 
lation and the relation between the decay of 
seismic waves and D. Moreover, as the pre- 
dominant period of earthquake-motions de- 
pends on its maximum amplitude, it is natural 
that the frequency response of the seismome- 
ter also has an effect on the form of the 
curve, 
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A simple model is desirable for furthering 
quantitative studies on this problem from the 
above explained point of view. One of the 
simplest models is that earthquakes are oc- 
curring uniformly on the surface of the crust. 
Now let us study the frequency diagram of 
P-S intervals obtaines by observing these 


earthquakes at a position on the surface. 


Here, A denotes the maximum trace amplitude | 
. ’ 

of an earthquake, and A,, the maximum am- | 

plitude of the same one at its hypocenter. 


N(Ao, 4)dAy d4 denoting the number of the 
earthquakes of Ay between Ay and Ay+dAo 
and of 4, or hypocentral distance, between 4 
and 4+d4, we have, 

No(Ao, 4)dApd4d=274p(Ap)dArd4 , (6.1) 
where p(Apy) is the frequency distribution of 
Ap in unit area. 


Putting, 
p(Ay)=kAo™, (6.2) 

we get, 
No(Ao, 4)d Acd4=27AkRAo™dAcd4 (6.3) 


where m and & are constants. 

Since the seismic waves decay as they travel, 
the magnitude of the minimal earthquake 
recorded at a position becomes larger with 
increasing hypocentral distance. Aomin. obs. 
denoting the Ay of the minimum earthquake 
which is observable by a seismometer having 
a certain sensitivity at the hypocentral dis- 
tance, 4, Aomin. obs. iS expressed by a function 
of 4, or 

Ao min. obs.=a@(4) 2 
The trace amplitude of the small earthquakes 
which can barely be identified as earthquake 
being denoted by Amin., 


Amin: =a , 
where @ is a constant which depends on the 
magnification of the seismometer and the 


amplitude of groundunrests at the station. 
Here we get, 


a= VAomin. obs. /(4) > (6.4) 


V being the magnification of the seismometer 
and (4) is the relation between the decay of 
seismic waves and 4. 

As 


A= A) Vf(4) (6.5) 
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we have, 
a(4)=a/Vf(4) . (6.6) 


N(4)d4 denoting the number of the earth- 
quakes, 4 of which is between 4 and 4+d4, 
we have, 


b(A) 
Naas=| getdkAs™ dA, dd (6.7) 


where we may put b(4)= approximately. 
Thus we get, 


oak 


N(4)d4= -4-a(4)-™dd (6.8) 


or, 


N(4)d4= 


bah m—1 m—1 
perme etki and sfid maida 


(6.9) 


The error in NM4)d4 introduced by this ap- 
proximation is below a few percents. 

The equation (6.9) shows that the form of 
the frequency distribution of P-S intervals 
_ depends on the magnitude frequency relation 
and the relation between the decay of seismic 
waves and the hypocentral distance. The 
total number of earthquakes observed at a cer- 
tain position during a certain period is deter- 
mined by the magnification of the seismo- 
meter, the magnitude frequency relation and 
the constant, & or the seismicity. When the 
magnitude frequency relation is expressed by 
the equation M(A))dAp=RA,.-"dAo, the total 
number of the recorded earthquakes is pro- 
portional to V™-!. In other words, the value 
of ‘‘m’’ in the equation, (4.2) usually being 
co or —)'* in (4.1). being. 0.9, the total 
number of the observed earthquakes will be 
63 times, when the sensitivity of the seismo- 
meter is increased 100 times. 

For deducing the frequency distribution of 
the P-S intervals of earthquakes, the foci of 
which are on a plane witha depth of h, when 
they are observed at a point on the earth’s 
surface we have to transform the variable, 
4 in (6.9) into D, where 4?+h?=D*. This 
relation is obtained under a simple assump- 
tion as to the wave trajectory. Transform- 
ing 4 into D in the equation, 


2rkai-™ 
m 


i amet f(A) dA : 


N4)p4= 
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we have, 
N(D)dD 
. .2rraS er rae | 
or, 
ork 
N(D)dD= song - Vin-l. D-f(D)"-1dD . 
(6.10a) 
In case D’<h, 
N(D’)dD'=0 . (6.10b) 


The equation (6.10a) has the same form as 
(6.9). 

Now, let us take upa more practical model. 
Earthquakes are occurring in a layer having 
a definite thickness, boundaries of which are 
respectively AH, and H, in depth. This model 
is similar to the one adopted by Kawasumr2” 
except for the assumption on magnitude fre- 
quency relation. 

Denoting the frequency distribution of P-S 
intervals of the earthquakes occurring in the 
above described layer and being observed on 
the surface of the crust by N,(D)dD, we get, 


Nir(D) dD=0, (6.11a) 
where, D<H;; 


Na(DdD=NDyaD\ dh 


T, 
= ND\(D—H,)dD , 
where, M,<.D<.H2; and 


(6.11b) 


Na(D\dD= ND)a DV “ah 
1 


= ND) H.—M;)dD , (6.11c) 


where, H.<D. 

In this case, the dimension of k in the equa- 
tion, p(Ay))=kAy-™, differs from that of & in 
the equations 6.9 and 6.10. 

Considering that the surface area of a belt 
which is a part of a sphere of radius D sliced 
by two parallel planes of distance h is 2xDh, 
it will easily be seen that the result is correct. 

Some might think that the above mention- 
ed model for earthquake occurrences is too 
simple and abstract for representing actual 
situations. However, in observing the earth- 
quakes occurring in the vicinity of Kanto dis- 
trict, the model seems to stand approximate- 
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Fig. 16. Frequency distribution of P-—S intervals calculated according to the formula 6, 11. 


a: m=2, f(D)01/D3, 
c: m=1.8, f(D)o1/D?, 


ly for the actual distribution of foci. 

In Fig. 17 the curves for the frequency dis- 
tribution of the P-S intervals of the earth- 
quakes recorded at Hongo and Mitaka in 
Tokyo are shown. Putting H,=4x7km, and 
H,=10x7km, and assuming the decay of the 
seismic waves to the proportional to the square 
of the hypocentral distance in the equation 
6.11, we get the results shown in Fig. 16. 
The calculation was made for both cases, 
m=2 and m=1.8. The similarity of the 
curves in Fig. 16 and those in Fig. 17 sug- 
gests that the above stated theory on P-S 
interval frequency distribution is fairly reason- 
able. 

The depth of the deeper boundary, 70 km, 
seems to be too large, but this value is ap- 
propriate in this case, because in the vicinity 


b: m=1.8, f(D)o1/D3, 
d: m=2, f(D)o1/D. 


of Tokyo Bay and Kinugawa area in Kanto 
district the foci are distributed form 30 km 
down to bout 70km in depth. Most of “ shal- 
low’’ earthquakes in and near Kanto district 
are occurring under the Moho-discontinuity, 
or in the upper part of the mantle. 


The frequency distributions of P-S intervals 
of both microearthquakes and earthquakes of 
M above 4 recorded at Mt. Tsukuba are shown 
in Table 7 and Fig. 13. They show that 
both are approximately similar except that 
the curve for microearthquakes has its maxi- 
mum between 10 sec and 11 sec and the one 
for larger earthquakes, its maximum between 
9 sec and 10 sec. 

48 microearthquakes among those which be- 
long to the class-interval of 10 sec and 11 sec 


P-S 10 20 30 40 50 60 
Sec 


a 
Fig. 17. 


arithmetic moving-average-method. 


are the aftershocks of the earthquake at 13h 
09m, Aug. 12, 1955 (M=5.4). As most of 
microearthquakes belonging to the class- 
interval between 10 sec and 11 sec in the 
frequency diagrams in Fig. 18 have the P-S 
intervals of 10.3 or 10.4 sec except these 
aftershocks, and the accuracy in determining 
the commencement of the S phase being low, 
we may neglect the small difference in the 
position of maxima in both diagrams. Now, 
we may safely conclude that the geographical 
distribution of the microearthquakes of M be- 
low 0 and that of larger earthquakes of M 
4 or 5 agree statistically with each other. 
There is another small discrepancy between 
the two geographical distributions. Some of 
the microearthquakes have the P-S intervals 
from 3 sec to 4 sec, but none of the earth- 
quakes of M above 4 has the P-S interval 
is less than 4 seconds. In other words, none 
of them has occurred in the region within 
30 km from Mt. Tsukuba. Since the depth of 
the Moho-discontinuity is about 25km in Kan- 
to?® district, it is safely concluded that the 
‘crust is more ‘‘stable’’ than the upper part 
of the mantle in the vicinity of Mt. Tsukuba. 
Eight earthquakes of M above 4 and of the 
P-S interval from 4 to 5 seconds were record- 
ed at Mt. Tsukuba during a period from 1924 
to 1931, but since 1931 none of earthquakes of 
M above 4 has had the P-S intervals shorter 
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Frequency distribution of the P-S intervals of shocks recorded by a displacement 
sismograph: a, at Hongo, Tokyo; and b, at Mitaka. 


The curves are smoothed by the 


than 5 seconds. On the contrary, about 2 per- 
cent of the microearthquakes recorded there 
have the P-S intervals shorter than 5 sec, 
or occur in a part of the crust and the mantle 
within about 35km from Mt. Tsukuba. 

In observing the frequency distribution 
diagrams of the P-S intervals in Fig. 18, we 
may safely conclude that in the vicinity of 
Mt. Tsukuba no earthquakes occur in the 
part of the crust shallower than about 20 km. 
In this part no microearthquakes were ob- 
served during the period of the present 
writer’s experiment. Accordingly, as already 
stated in Chapter 4, it is inferred that 
the earthquakes of M larger than 4 will not 
occur there more than once in several years. 
Actually, in the crust near Mt. Tsukuba earth- 
quakes of M above 4 have not occurred in 
the past 20 years since the beginning of seis- 
mological observations at the station. This 
fact seems to support the reasoning described 
in chapter 4 and 5. 

The curves in Fig. 19 suggest that the ge- 
ographical distribution of the foci of the sen- 
sible earthquakes observed at Mt. Tsukuba 
did not change sensibly during a period be- 
tween 1924 and 1941. 


§7. Résumé 


The results of the present experiment and 
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Table 7. Frequency distribution of P-S intervals 
of shocks recorded at Mt. Tsukuba. 


No. of micro- 


Number of felt 


'1924~ 


1941 


Sia bk 


earthquakes 

ti BE a 

ene hours in total 1930 1937 1941 

in Aug. 1955 | 

3.1~ 4.0 3 Gero 0 
4.1~ 5.0] 3 8} Oo 0 
5.1~ 6.0 19 38 | 39 «16 
6.1~ 7.0 14 27.| 426 7 
7.1~ 8.0 a7 27| 29] 10 
8.1~ 9.0 28 35-1400 GTi 
9.1~10.0 | 27 60| 56| 22 
10,1~11.0 30-243 4) 923) 292 9 
11.1~12.0 | 26 Dp e20 6 
(7113.0 9 19} 16 8 
13.1~14.0| 9 9} 8| 7 
14.1~15.0 3 12) iit 3 
15.1~16.0 | 5 11 9 4 
16.1~17.0 2 9 walLtoala 
17.1218.0 3 1 ag a 5 
18.1~19.0 3 4 4 5 | 
19.1~20.0 0 1 4 6 
20.1~21.0 2 1| 10 7 
21.1--22:0 4 ine 2 
22.1~23.0 3 T | ae 5 
23.1~24.0 1 0 2 | 2 
24.1~25.0 0 1; 8 | 2 
25.1~26.0 1 0 | 4 2 
26.1~27.0 1 o| | 2 
27,1~28.0 2 et ee 
28.1~29.0 1 OSLO 
29.1~30.0 2 Oo} 3 2 
30.1~31.0 2 o| “oh we 

>31.0 24 | 


a 


their analysis have led us to the followin 
conclusions: 

1) The present writer made a tempora 
observation of near-by earthquakes at Mt 
Tsukuba by Ultra Sensitive short period seis 
mometers. Their sensitivity is about 10°mm 
kine, or the magnification is 1.3 x 10’ at 20 cps 
The seismometers are of velocity meter type, 
and their frequency response curves are flat 
between 10 cps and 200 cps. For obtainin 
good records of very small near-by earth+ 
quakes it is essential that the seismometers 
are sensitive for the earthquake-motions o 
frequency tens of a cycle per second. The 
seismometers must be placed on bed rock. 

2) About 200 microearthquakes were re- 
corded at Mt. Tsukuba every 24 hours. Th 
hypocentral distances of most of them are 
from 30km to 100km. They are not after-: 
shocks of a particular large earthquake. They 
seem to occur approximately stationarily., 
Their earthquake-motions have the prevailing 
periods of 1/20 sec or less. 

3) A» denoting the maximum trace ampli-- 
tude, the frequency distribution of A, can) 
be expressed by an equation 

NAv)dAv=kRAv ™dAd , 
where, m is 1.88. This relation is valid for: 
a range of A, from 10-*cm/sec to 10-2cm/sec: 
or more. 

4) The energies of earthquakes, which are! 
recorded by the seismometers described above: 
and hypocental distance of which are less: 
than about 100 km, are from 10 ergs to 10%! 
ergs. The magnitude of the smallest one is: 
estimated about —1. The result described in. 
3) is equivalent to the following relation; the. 
magnitude frequency relation can be express- 
ed by the equation for a magnitude range of 
from —1 to 5. 

NM)dM=const. x 10-°* dM , 


where, }=0.8. The relation between m and 
b has been discussed 


5) The number of earthquakes of magni- 
tude 4 or 5 possibly occurring in the vicinity 
of an observing station during one year can 
be calculated according to the equation given 
in 3) or 4), when the number of microearth- 
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_ period of tens of hour are known. 


shocks recorded at Mt. Tsukuba. The upper dia- 
gram is the frequency distribution of P-S intervals 
of microearthquakes, and the lower one is that 
of felt earthquakes observed during a period bet- 
ween 1924 and 1941. 


quakes recorded at the same station during a 
At Mt. 


Tsukuba the rate of occurrence of microearth- 


“quakes are about 200 for every 24 hours. 


~ less than once per each 24 hours. 


Accordingly the equation in 3) or 4) shows 
that about 20 or 30 earthquakes will be felt 
there by people every one year. Actually, at 
Mt. Tsukuba, 45 earthquakes of magnitude 
about 4 or 5 were felt in 1955. 

6) The results obtained in the present ex- 
periment lead to the conclusion that the earth- 
quakes of M 4 or 5 do not occur more than 
once in several years in a part of the crust 
where microearthquakes (E~10" ergs) occur 
For in- 
stance, in the vicinity of Mt. Tsukuba no 


Fig. 19. Frequency distribution of P-S intervals 
of felt earthquakes observed at Mt. Tsukuba. 
The curves are made from the data obtained 
during three different intervals of time: 1924— 
1930: 1931-1937: 1938-1941. 


microearthquakes have occurred in the part 
shallower than 20km. in depth during the 
term of the experiment; in this part, the 
the earthquakes of M larger than 4 have not 
occurred since 1924. It is noteworthy that 
according to the historical study, no destruc- 
tive earthquakes have occurred in the vicinity 
of Mt. Tsukuba. 

7) The facts described in 5) and 6) suggest 
that the observation of the microearthquakes 
with Ultra Sensitive seismometers provides 
valuable data for the statistical study on 
the occurrence of destructive earthquakes. 
More concretely speaking, we may say that 
we can find out a part of the crust where 
destructive earthquakes will never occur, by 


112 


making the observation of microearthquakes 
and determining whether they occur there or 
not. Even in Japan where destructive 
earthquakes occurs every 2 or 3 years, we 
can find the areas where no destructive earth- 
quakes have occurred since the beginning of 
history. 

8) The geographical distributions of the 
seismicity of microearthquakes of FE about 
10” ergs or larger agrees approximately with 
that of earthquages of M 4 or 5. As the 
microearthqakes are apt to occur in swarms, 
the period of the observation of microearth- 
quakes should be sufficiently long in order to 
avoid distorted results concerning the geo- 
graphical distributions of their seismicity, or 
the period of the temporary observation by 
Ultra Sensitive seismometers should be longer 
than a week. 

9) The results of the present experiment 
is insufficient for concluding that the micro- 
earthquakes of E smaller than 10” ergs still 
exist. Smaller microearthquakes would prob- 
ably be recorded by seismometers which are 
more sensitive than the ones used in the 
present experiment, if they could be operated 
at a more quiet place, although it is technical- 
ly extremely difficult. 


The present writer and his collaborators in- 
stalled the Ultra Sensitive seismometers at 
the Matsushiro Seismological Observatory in 
August in 1956. Microearthquakes of P-S 
intervals about 1 sec were recorded, and the 
smallest of them had the energy E~109 ergs. 
The experiment at Matsushiro will be re- 
ported in another paper. 
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